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PREFAZIONE 


Il linguaggio BASIC conosce attualmente un grande sviluppo data la diffusione 
di microcalcolatori personali e professionali. 


Questo linguaggio è certamente il più diffuso in questo campo: è utilizzato per la 
gestione, l’acquisizione dei dati, i calcoli scientifici e per le applicazioni grafiche ed 
il controllo di processo. 


Bene, la programmazione si impara essenzialmente con la pratica. 


Ecco una collezione completa di esercizi concreti e completamente svolti che 
potrete seguire passo passo: enunciazione del problema, analisi, flow chart e 
commenti, programma, esempio di esecuzione. 


Noi speriamo che dopo aver effettuato un numero sufficiente di questi esercizi e 
dopo averne letto la spiegazione, voi avrete effettivamente appreso ‘‘praticamente 
... Il BASIC". 


Ciascun esercizio proposto è trattato in modo completo: 
@ enunciazione del problema 
@ analisi 
@ flow chart con commenti 
@ programma 
@ esempio di esecuzione 


Questo metodo permette al lettore di verificare ad ogni passo quali progressi avrà 
compiuto nell’apprendimento dell’analisi di un problema. Mostra inoltre come 
scomporre il problema in “sotto-problemi‘ più semplici. Così è possibile risolvere 
separatamente ogni difficoltà e costruire un programma modulare di facile lettura e 
di facile adattabilità ad applicazioni concrete. Questi sotto-problemi possono in 
seguito servire da base per la costruzione di programmi più complessi. 


Questo libro costituisce un complemento necessario ad una buona assimilazione 
del linguaggio. Un’appendice presenta in maniera coincisa le principali istruzioni 
del BASIC. 


Alfine di rendere generalizzati i programmi l’autore ha utilizzato nella maggior 
parte dei casi, un BASIC molto diffuso (sotto-insieme dell’MBASIC della MICRO- 
SOFT). Per alcuni programmi sono proposte più versioni alfine di mostrare come 
sia possibile aggirare le lacune di certi BASIC. 
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CAPITOLO 0 


IL VOSTRO PRIMO PROGRAMMA IN BASIC 


0.0 INTRODUZIONE 


Questo capitolo mostra al lettore che ciascuno può imparare a programmare 
partendo da esempi semplici e che, tranne per una certa categoria di problemi, la 
programmazione non é cosa da specialisti. Per comprendere questo capitolo, non é 
dunque necessario conoscere il BASIC e, partendo da un semplice esempio, il 
lettore apprenderà differenti categorie di istruzioni BASIC e, come sia possibile 
migliorare progressivamente un programma. Il tema scelto é il calcolo del bilancio 
familiare che seppur di limitato interesse é comunque utile per applicazioni più 
interessanti quali il calcolo dell’imposta sul reddito delle persone fisiche (IRPEF) 
che costituisce l’oggetto di due esercizi nel Capitolo 6. 


Per migliorare la conoscenza del linguaggio, il lettore potrà da una parte leggere 
un libro teorico che presenti in maniera più precisa le differenti istruzioni del 
linguaggio, certamente in modo meno gradevole di questa, dall’altra fare degli 
esercizi che sono indispensabili nell’apprendimento della programmazione. 


0.1 CALCOLO DEL BILANCIO FAMILIARE 


Desideriamo ora calcolare il bilancio familiare di una famiglia partendo dal 
reddito imponibile e dal numero dei componenti della famiglia stessa. Chiamere- 
mo: 

B bilancio familiare da calcolare 
R reddito imponibile 
N numero dei componenti la famiglia 


Per ottenere il bilancio familiare B é sufficiente effettuare la divisione 
R/N. Proporremo quindi un programma in grado di svolgere le seguenti operazio- 
ni: 


leggere R e N 20 INFUT R,N 
calcolare Q= R_ 80 Q=R/N 
40 PRINT Q 
N 
stampare Q 50 END 


Da questo piccolo programma potremo trarre le seguenti constatazioni: 


— ogni linea comporta un numero di linea, 

— ogni linea comporta un’istruzione, 

— l’istruzione di lettura INPUT corrisponde all’impostazione di una 
informazione nel computer, 

— l’istruzione di calcolo B= R/N si scrive su di una sola linea e la barra 
obliqua é il simbolo della divisione. Questa barra é chiamata 
“SLASH*, 


— il programma termina con l’istruzione END. 


L’esecuzione del programma fornisce il seguente risultato: 


? 100000,3 
33333.,3 


Quando il calcolatore deve eseguire l’istruzione INPUT, stampa sul terminale un 
punto interrogativo che sta a significare l’attesa di dati da parte dell’operatore. In 
questo esempio di esecuzione, l’operatore ha battuto due numeri separati da una 
virgola: 100000 e 3. Poiché la lista delle variabili che segue INPUT é R, N, il primo 
valore é attribuito a R ed il secondo a N. Il computer effettua la divisione e la 
stampa del risultato. 


Questo risultato é corretto dal punto di vista matematico, ma la presentazione é 
poco elegante. Comanderemo ora al computer di stampare un commento modifi- 
cando leggermente il programma. 


Per far stampare un testo, é sufficiente piazzarlo tra virgolette nell’istruzione 
PRINT. Cerchiamo ora di migliorare il programma stampando un commento 
all’inizio dello stesso. 


nre pei ii. virgola impedisce 
10 FRINT "DARE I VALORI DI R ED N" : ; : 
20 INFUT R,N il passaggio 
30 Q=R/N i 
40 FRINT "EILANCIO FAMILIARE = "3E alla linea seguente 
50 END 


Esempio di esecuzione: 


DARE I VALORI DI DI KR ED N ? 100000,8 
BILANCIO FAMILIARE «= 33333,3 


Quando il computer ‘‘attende‘ i dati, anticipa l’utilizzatore 
ponendo in generale un punto di domanda. 


Con alcuni computer, in particolare con numerosi home-computer, si otterrà lo 
stesso risultato scrivendo: 


20 IINFPUT"DARE T VALORI DI KR ED N",KR,N 


In questa forma l’istruzione INPUT stampa un testo prima che il calcolatore 
legga i dati introdotti dall’operatore. 


0.2 SECONDO CALCOLO DEL BILANCIO FAMILIARE 


Consideriamo ora il caso dell’impiegato a stipendio fisso. I dati di base sono il 
reddito netto dal quale dedurremo forfettariamente un 10% per renderlo più reale, 
poi applicheremo le altre deduzioni quali: l’assicurazione sulla vita, interesse di un 
mutuo, etc., il calcolo del reddito imponibile é dunque: 


reddito imponibile = 0.80 x reddito netto +pensioni e redditi vitalizi 
a titolo gratuito - deduzioni 


Il reddito imponibile é quindi utilizzato per il calcolo del bilancio familiare. 


100 FRINT'IMPOSTARE IL REDDITO NETTO"? 

110 INFUT RI 

120 PRINT'IMFOSTARE IL MONTANTE DELLE PENSIONI E LE RENDITE VITALIZIE!"; 
130 INPUT F 

140 FRINT"'IMFOSTARE IL MONTANTE DELLE DEDUZIONI": 

150 INFUT D 

160 Re, B80XFE14+P-D 

170 PRINT"'IMFOSTARE IL NUMERO DEI COMPONENTI Là FAMIGLIA"; 

180 TINFUT N 


190 N 

z00 NT 

210 SINT"REDDITO IMFONIETLE= "R$" BILANZIO PAMILIARE= "E 
220 END 


Esempio di esecuzione: 


IMPOSTARE LA RENDITA NETTA 2 140990 

IMPOSTARE Il. MONTANTE DELLE .FENSIONI E DELLE RENDITE VITALIZIE ? 1000 
IMFOSTARE IL MONTANTE DELLE DEDUZIONI ? 7000 

IMPOSTARE IL NUMERO DEI COMPONENTI LA FAMILIA ? 4 

REDDITO IMFONIETIL 109200 

EILANCIO FAMILIARE= 27300 


0.3 CONCLUSIONI 
Si può ora modificare il programma per fare in modo che consideri: 


— il numero di adulti ed anziani che incidono ciascuno per una certa 
parte 

— il numero di bambini che incidono ciascuno per metà 

— il reddito netto ordinando al computer di eseguire il calcolo delle 
deduzioni 


Si può inoltre domandare al calcolatore di fare il calcolo del 
montante delle imposte. Questo sarà trattato nel Capitolo 6. 


Con questo semplice esempio si potrà facilmente proseguire il lavoro. Per altri 
problemi invece, saranno necessarie nozioni più precise concernenti per esempio i 
flow-chart che saranno presentati in seguito. 


Come tutti i linguaggi di programmazione, il BASIC é costituito da un “alfabe- 
to proposto nell’appendice. 


CAPITOLO 1 


FLOW-CHART 


1.0 INTRODUZIONE 


Il capitolo precedente ha mostrato come si possano facilmente apprendere i 
rudimenti del linguaggio BASIC, che ci ha permesso di scrivere un semplice 
programma. Tuttavia man mano che i problemi diventano più complessi, non é più 
possibile scrivere direttamente il programma. E’ necessario prima analizzare il 
problema quindi disegnare un ‘‘flow-chart“. In effetti l’esperienza ha mostrato che 
questo costituisce un aiuto prezioso per la programmazione, in particolare per chi é 
all’inizio di questa attività. 


Lo scopo di questo capitolo é mostrare come si costruisce un flow-chart. I 
capitoli seguenti permetteranno al lettore di applicare ciò che avrà appreso ed 
inoltre di completare le proprie conoscenze. 


In seguito, con l’esperienza, diviene possibile ridurre il tempo per disegnare il 
flow-chart, questo é comunque da sconsigliare al neofita. 


1.1 SCOPO DEL FLOW-CHART 


Il flow-chart costituisce una rappresentazione grafica del problema trattato. Lo 
stato attuale della tecnica non consente al calcolatore di comprendere questa 
rappresentazione, ciò pone, dunque, il problema della sua utilità. 


Il flow-chart permette di verificare che al momento dell’analisi del problema, 
certi casi non siano stati dimenticati e facilita il dialogo tra diverse persone che 


partecipano alla elaborazione del programma. Per il neofita il flow-chart di detta- 
glio costituisce una tappa preliminare che aiuta ad una buona programmazione. 


1.1.1 DIFFERENTI TIPI DI FLOW-CHART 
Nella pratica si distinguono tre tipi di flow-chart: 


— flow-chart di catena: soprattutto utilizzati in gestione, mettono in 
evidenza i legamenti tra archivi e programmi, 


— flow-chart di principio: molto utili per descrivere in modo macro- 
scopico grossi programmi che comportano un ‘‘concatenamento di 
algoritmi‘, é di interesse molto limitato per piccoli programmi, 

— flow-chart di dettaglio: costituiscono una rappresentazione precisa 
e completa dell’argomento trattato, tolgono ogni eventuale ambi- 
guità e permettono una facile programmazione. 


Tuttavia il flow-chart deve restare il più possibile indipendente dal linguaggio di 
programmazione utilizzato. 


1.1.2 NORME 

Le tecniche del flow-chart sono state l’oggetto di una norma AFNOR Z 67-010 
che é stata omologata nell’ Aprile 1966. Alla fine del capitolo troveremo i principali 
simboli utilizzati per il disegno del flow-chart. 

NOTE: 
1. In italiano la traduzione del termine flow-chart é DIAGRAMMA DI FLUSSO. 


2. Esistono altri metodi per la rappresentazione degli algoritmi, ma questi presup- 
pongono un’approfondita conoscenza delle tecniche di programmazione. 


Vediamo allora un semplice esempio di mini flow-chart, ciò permetterà di 
constatare che la soluzione non é una sola prima dello studio del problema. 


1.2 IL VALORE PIU’ ALTO DI DUE NUMERI A E B 


Supponiamo che X assuma il valore più alto dei due numeri A e B. Come 
procedere per minimizzare la stesura? 


1.2.1 Prima soluzione: si raffrontano A e B, se A = Ballora si pone il valore di A 
nella X altrimenti si pone il valore di B nella X. 


Questo metodo può essere rappresentato per mezzo del flow-chart di Fig. 1.1 che 
contiene un ‘“rombo* nel quale si trova il confronto A = B, e due “rettangoli‘* che 


corrispondono a delle ‘‘assegnazioni“‘. 


Il listato di Fig. 1.2 corrisponde alla sequenza BASIC relativa. 


100 IF A >» = E THEN 1530 
110 X os E 

120 GOTO 140 

130 X à 

140 seguito 


Figura 1.1 Figura 1.2 


Utilizzando un BASIC più evoluto si potrà scrivere 
100 IFÀA: BO THEN X = A ELSE X E 


1.2.2 Seconda soluzione: per evitare il salto della linea 120, si può modificare il 
flow-chart di Fig. 1.1 sostituendo una delle istruzioni di assegnazione, da qui il 
flow-chart di Fig.1.3 che potrà essere scritto in BASIC. 


110 IF À ® E THEN 130 


Figura 1.3 Figura 1.4 


Sotto questa forma si economizza nel listato un GOTO. Questo permette di 
ridurne un pò l’ingombro. 


Con un BASIC evoluto si avrà: 
100 Xs A 


110 1F [Ri AUTHEN Xo = E 


1.2.3 Terza soluzione: certi interpreti BASIC, sfortunatamente rari, comportano le 
funzioni MAX e MIN. F' sufficiente allora scrivere: 


100 X = MAX CALLE 


Attualmente queste funzioni MA X e MIN sono raramente disponibili sui perso- 
nal computer. 


NOTA: quando queste funzioni sono disponibili, ammettono un numero qua- 
lunque di parametri: si può quindi scrivere per esempio 


Yo= MAY 0X,9,Z,0) 
oppure ancora: 
Yo MINO OX + Z, xk, KXSINCA)) 
1.3 ESEMPIO DI FLOW-CHART: ELEMENTO MAGGIORE DI UNA 
TABELLA 


Supponiamo di cercare l’elemento più grande di una tabella A di 100 numeri. Il 
metodo proposto é il seguente: 


- porre X= A(1) 
- assegnare successivamente a I i valori 2,3,4 fino a 100: 


- confrontare X e A(I) 
- se X< A(I) trasferire il valore di A(I) nella X altrimenti continuare. 


Quando si é terminato, X conterrà il valore più alto. Questo metodo può essere 
rappresentato con il flow-chart seguente 


—— lettura 


LEGGI A 


*——- calcolo 


<«—__ per ripetere 
i confronti 
confronto 


Figura 1.5 


Questo schema pone in evidenza che: 


— le istruzioni di Input/Output sono rappresentate da un parallelo- 
gramma, 

— le istruzioni di calcolo sono rappresentate come un rettangolo 

— le istruzioni di confronto sono rappresentate da un rombo. 


Constatiamo inoltre un tipo di istruzione curiosa per l’‘‘apprendista‘“ 
I=1I+1 
In maniera più generale un’istruzione di calcolo si scriverà: 
variabile = < espressione > 


Questa istruzione significa calcolare il valore numerico dell’espressione e porre il 
risultato nella variabile situata a sinistra del segno = (uguale). Per questa ragione 
detta istruzione é chiamata di ‘‘assegnazione‘‘ ed il carattere= (uguale) rappresenta 
appunto il simbolo dell’assegnazione. 


Per contro in un rombo I= 100 significa confrontare I e 100 e vedere se essi hanno 
il medesimo valore. In nessun caso bisogna porre il valore 100 nella I. Nel rombo, il 
carattere = (uguale) é un simbolo di confronto. 


1.4 COME VERIFICARE UN FLOW-CHART 


Se si procede alla programmazione partendo da un errato flow-chart, il risultato 
non sarà certo quello auspicato. E’ dunque consigliabile, dove possibile, assicurarsi 
dell’esattezza del flow-chart prima di cominciare la programmazione. 


Per questo si può far ‘“girare‘ il flow-chart ‘a mano“, in questo caso l’essere 
umano simula la macchina che funzionerà nel rispetto di quanto previsto nel 
flow-chart potendo comunque intervenire ove fosse necessario. 


Riprendiamo il flow-chart di Fig. 1.5 e interessiamoci ad una tabella più piccola, 
per esempio di 5 numeri. 


Figura 1.6 


All’inizio poniamo X = A(1), X assumerà quindi il valore 3. Possiamo ora far 
girare il loop di calcolo. La tabella 1.7 mostra l’evoluzione del contenuto di X in 
funzione di I. 
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si confronta A(2) con X e X è più grande. 
poichè A(3) è più grande di X si pone A(3) nella X 
— lè più piccolo di 4 


6 è più grande di X quindi si copia il 6 nella X 


Figura 1.7 


Constatiamo che X si sposta verso il valore più alto della tabella. Si può dunque 
programmare questo flow-chart. 


NOTA: Questo metodo si può applicare solo con flow-chart molto semplici. Il 
flow-chart della Fig.1.5 può essere tradotto in BASIC in differenti modi, per 
esempio: 


100 DIM ACILGO) 

110 FOR I = 1 TO 100 

1200 READ ASI) 

130 NEXT XY 

140 XY = A(1) 

G0 FOR I = 2 TO 100 

160 IFOX 2 ss ACI) THEN 180 
170 Xos ACI) 

180 NEXT I 

190 PRINT "MASSIMO NUMERO DELLA TABELLA= "3 X 
200 DATA 0), 

210 DATA L0606, 

410 END 


Figura 1.8 


Questa forma non è la migliore, ma è di facile comprensione: 
— le linee da 110 a 130 sono destinate a far leggere tutta la tabella, 


— le linee da 140 a 180 inclusa, corrispondono alla ricerca del valore 
più alto, 


— le linee da 200 a 400 contengono normalmente i valori numerici dei 
100 elementi della tabella. 
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Critica a questo programma 


Questo programma non può funzionare correttamente se la tabella non contiene 
esattamente 100 elementi. Di solito si preferisce far leggere prima il numero n degli 
elementi che costituiscono la tabella ed in seguito prevedere un programma capace 
di adattarsi al numero n di elementi. Il programma della Fig.1.9 è, da questo punto 
di vista, senz'altro migliore. 


100 DIM ACLOO) 

105 READ N 

110 FOR I = 1 TON 

120 READ ACI) 

130 NEXT I 

140 X = ACL) 

150 FORT = 2 TON 

160 IF X 3 o0# ACI) THEN 180 
170 Xo= ACI) 

180 NEXT I 

190 PRINT "MASSIMO NUMERO DELLA TABELLA= ", 
200 DATA 5 

210 DATA 3,-2,34,5,0 

410 END 


bas 


Esempio di esecuzione: 
MASSTMO NUMERO DELLA TAEELLA= 34 


Figura 1.9 


Commenti: 


— l’istruzione 105 permette di leggere il numero “n‘ degli elementi della tabella, 


— la linea 200 contiene il valore numerico 5 corrispondente qui a 5 elementi, 


— la linea 210 contiene i 5 valori numerici, 


— sotto questa forma, questo programma è limitato dall’istruzione DIM A(100) 
a 100 elementi. Modificando questa dichiarazione, è possibile adattare il 
programma a tabelle di maggiori o minori dimensioni, 


— questa forma riduce le modifiche del programma consentendone inoltre una 
migliore leggibilità. 


In generale è consigliabile che il valore di chiusura del loop FOR, sia 


una variabile piuttosto che una costante. 
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Gli elementi di un flow-chart 


Simbolo generale di trattamento 


Chiamata ad un sottoprogramma 


Simbolo di test 


Rinvio 


Inizio e fine programma 


Input/Output 
simbolo generale 


Input da tastiera 


Output su carta (stampante) 


[1 
SD 
O 


=) GO) 
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Un altro metodo più coinciso consisterà nello scrivere: 


10 DIM AC100) 


15 N = 100 

20 MAT READ A 

sì ca Aa La attenzione la funzione MAX 

50 i = MAX(XACI)) non è disponibile 

50 NEXT IO su tutti i sistemi 

70 PRINT"!MASSIMO NUMERO DELLA TABELLA? "3X 

80 

90 

. istruzioni DATA 

160 

170 END 

Figura 1.10 

1.5 I TEST 


In un flow-chart, un test ha un’entrata ed in generale 2 0 3 uscite. La rappresenta- 
zione si ottiene nel modo seguente: 


, SI NO 
2 uscite 


3 uscite < > 


Il simbolo :: è utilizzato come simbolo di confronto. 
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Si può concepire a livello di flow-chart un test con più di 3 uscite poichè il 
flow-chart è una rappresentazione dell’algoritmo. La norma non prevede la rappre- 
sentazione per test che abbiano più di 3 uscite. Si può per esempio utilizzare la 
seguente rappresentazione: 


1.6 UNA TECNICA: LA BILANCIA 


Come rappresentare un flow-chart di dettaglio in modo che la parte sinistra sia 
percorsa ad ogni loop dispari e che la parte destra sia percorsa ad ogni loop pari, 
fino al termine del loop stesso? 


DISPARI 


Figura 1.11 
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Un metodo semplice consiste nell’utilizzare una variabile ausiliaria, il cui valore 
permetterà di effettuare questa “bilancia‘‘. Per esempio si attribuirà alla variabile B 
il valore 0 prima di entrare nel loop. Successivamente un test sulla variabile B farà 
andare a sinistra se B vale 0. Nella posizione di sinistra si aggiungerà un’istruzione 
B=1 in modo che al test successivo ci si sposti nella parte destra. In questa parte si 
aggiungerà l’istruzione B=0 che per avere al giro seguente uno scambio verso 
sinistra. Questo porta al flow-chart di Fig. 1.12. 


Figura 1.12 
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999 


1.000 


1500 


2 000 


IF B=0 THEN 1 500 


B=0 

PARTE B 
GO TO 2 000 
B=1 

PARTE A 

PARTE C 


IF 000 THEN 1000 


Figura 1.13 


A partire dal flow-chart di Fig.1.7, la programmazione in BASIC è facilmente 
realizzabile nella sequenza mostrata dalla Fig.1.13. I numeri di linea sono dati solo 


a titolo di esempio. 


NOTA 1: con il FORTRAN 77 si potrà scrivere: 


999 B=0 
1000 IF (B.EQ.0) THEN 
B=1 
partie A 
ELSE 
B=0 
PARTIE B 
ENDIF 
Partie C 
IF ( ) GO TO 1 000 
Figura 1.14 


NOTA 2: Con il CBASIC (1) si potrà scrivere: 


999 B=0 
1.000 IFB=0 THEN 
B=1: 
Partie A 
ELSE 
B=0 
Partie B 
Partie C 
IF THEN 1000 
Figura 1.15 


1.7 CONCLUSIONI 


Non abbiamo riporta- 
to i numeri di linea, 
che sono si obbliga- 
tori, ma, comunque 
non necessari alla 
comprensione. 


In questo capitolo abbiamo visto gli elementi di base che permettono di tracciare 
un flow-chart. L’appendice presenta un caso generale la cui conoscenza non è utile 
per il ‘‘neofita‘* ma che è molto interessante per coloro i quali s’interessano di 
problemi più complessi. Studiando gli esercizi dei capitoli seguenti, il lettore potrà 
perfezionare la propria conoscenza dei flow-chart e soprattutto il modo per realiz- 


zarli. 


(1) Il CBASIC costituisce il marchio depositato del BASIC esteso commercializzato dalla società 
americana SOFTWARE SYSTEMS e che funziona sotto il sistema operativo CP/M, disponibile sui 


sistemi compatibili INTEL 8080 Z80, INTEL 8080, INTEL 8085. 
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APPENDICE 


Generalizzazione di una bilancia a p posizioni 


Bisogna ora che il primo percorso sia fatto nel ramo I, il secondo nel ramo 2, il 
p-esimo nel ramo p, il p+l-esimo nel ramo l e così di seguito. In generale si dirà che il 
percorso dev'essere fatto nel ramo i modulo p. 


RAMO 1 RAMO 2 


Figura 1.16 


Non è più possibile rappresentare il metodo migliore in maniera coincisa per 
mezzo di un flow-chart, poichè sarà più opportuno utilizzare un’istruzione di 
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GOTO calcolato di Fig.1.17 anzichè una serie di test. Una prima soluzione è data 
dalla sequenza 1.17 


999 
1 000 
1 100 


1 300 


1 500 


1 600 


1 800 


B=1 
ON BGO TO I 100, 1 300, 1 500, 1 600 
B=2 

Ramo 1 
GO TO 1 800 
B=3 

Ramo 2 
GO TO 1 800 
B=4 

Ramo 3 
GO TO 1 800 
B=1 

Ramo 4 
Parte C 


IF... THEN 1 000 


Figura 1.17 


Si constata allora che è possibile accorciare la sequenza utilizzando l’istruzione B 
= B+ | posizionata nel “tronco‘ comune, tale istruzione fornisce la sequenza di 


Fig.1.18. 


999 
1000 
1 100 


1 300 


1 500 


1 660 


1 800 


B=1 
ON B GO TO 1 100, 1 300, 1 500, 1 600 


Ramo 1 
GO TO 1 800 


Ramo 2 
GO TO 1 800 


Ramo 3 
GO TO 1 800 


Ramo 4 
B=B+1 
IFB>4THENB=1 

Parte C 
IF ... THEN 1 000 


Figura 1.18 
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CAPITOLO 2 
ESERCIZI CONCERNENTI I NUMERI INTERI 


2.0 INTRODUZIONE 


Questo capitolo propone degli esercizi relativi ai numeri interi. I flow-chart non 
sono sempre facili da costruire e, da questo punto di vista, questi esercizi presen- 
tano un interesse didattico. Si consiglia al lettore, in caso di difficoltà, di non 
fermarsi troppo su certi esercizi, ma di passare agli altri capitoli per ritornarci in un 
secondo tempo. 

Le soluzioni proposte valgono per gli interpreti BASIC “standard‘*. Tuttavia 
sempre più frequentemente si trovano sistemi sui quali l’interprete BASIC accetta 
variabili intere, la cui rappresentazione si ottiene con una lettera seguita dal 
carattere %. 


Esempio A%, B% 


Per questi esercizi l’utilizzo di variabili intere che non superino il valore di + 
32767 (1), permette di ridurre sia il tempo di esecuzione che lo spazio in memoria. 


Comunque non tutti i sistemi accettano le variabili intere. Inoltre, insieme a 
queste variabili le funzioni standard SIN, COS, SQR, etc. sono raramente disponi- 
bili. 

Una difficoltà che si incontra in questo tipo di problema è la necessità di 
effettuare delle ‘divisioni intere‘* e di calcolarne il resto. Per esempio si desidera 
conoscere Q e R tali che: 


A=BxQ+R 


Per questo bisogna effettuare la divisione intera A/B per la quale il BASIC non 
dispone di operatori particolari. Si utilizzerà dunque la funzione INT e si scriverà: 


Q=INT(A/B) 
R=A-QxB 


(1) Questo limite è proprio della maggior parte dei personal computer. Tale limite è più elevato quando si 
utilizzano dei ‘‘mini‘‘* o dei grossi calcolatori. 
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Così si otterà il divisore Q e il resto R. Quando si ha la possibilità di utilizzare 
delle variabili intere, è sufficiente scrivere: 


Q%R=A%/B% 
R%A=A%-B%*Q% 


oppure se interessa solo il resto: 


— nel primo caso: 
R=A-B*INT(A/B) 
— nel secondo caso: 


R% = A% - B®R*(A%/B%) 


2.1 NUMERI INTERI COME a? + b° = c° 


Problema: Trovare tutti i numeri a e d interi compresi tra 1 e 100 in modo che a° + 
b° sia uguale ad un quadrato perfetto. 


Si chiede: 


— di analizzare il problema, 
— di stabilire un metodo e disegnare il flow-chart, 
— di scrivere il programma BASIC corrispondente. 


Soluzioni: Per prima cosa notiamo che sarà necessario eliminare una delle due 
soluzioni identiche salvo una permutazione dei valori. Per esempio: 


a=3 a=4 
b=4 e È =3 costituiscono due soluzioni identiche 
c= 5 c= 5 


per evitare questo si cercheranno dei numeri d che siano > a. Si potrà quindi far 
mutare una variabile intera i da 1 a 99 e la seconda variabile j da è + 1 a 100. 

Bisogna ora sapere se i + j è un quadrato perfetto, vengono presentati due 
metodi: 


I° metodo. Si procederà con prove successive su una variabile 4 alla quale si 
faranno assumere dei valori numerici partendo da j + 1. 


Se si trova i + j" = K° si è ottenuta una soluzione 
Se si trova il + j > K? proseguire le prove incrementando & di 1 


Se si trova i + j < K° la coppia i, j è errata. 
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Flow-chart 


Secondo metodo Primo metodo 
nio Rappresentiamo di seguito solo la 
parte che differisce dal secondo. 


K=INT(K) 


Figura 2.1 Figura 2.2 
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2° metodo. Si calcola: 


se K è intero allora si è trovata una soluzione, altrimenti la coppia i, j è errata. Per 
sapere se è intero, è sufficiente comparare X ed INT (k). Siamo così giunti ai 
flow-chart Figg.2.1 e 2.2 nei quali constatiamo: 


— la presenza di un loop di calcolo su J per fare variare J da I+1 a 100, 


— la presenza di un loop di livello superiore per fare variare I da 1a 99. 


Partendo dai flow-chart di Figg.2.1 e 2.2, si potranno facilmente scrivere i 
programmi riportati in Figg.2.3 e 2.4. 


5 N=100 lista del programma con il secondo metodo 
10 FORI = 1 TO N-1 
20 J=I+1 TON 
30 K=SQRCIXI+JXJ) 
40 IF K <> INT(K) THEN 60 
50 FRINT I}4}K 
60 NEXT J 
70 NEXT I 
80 END 
Figura 2.3 
5 N=100 lista del programma con il primo metodo 
10 FORI = 1 TO N 
20° FOR J = I+1 TO N 
30 S=IXI+JXJ ? K=dJ 
40 {=K+1 3 K2=KXK 
50 IF K2 « S THEN 40 
60 IF K2 > 5 THEN 80 
70 PRINT I}J}K 
80 NEXT J 
90 NEXT I 
100 END 
Figura 2.4 


Per facilitare la lettura di un programma, consigliamo al lettore di scrivere in 
“modo strutturato” scalando di circa due caratteri le istruzioni riguardanti i 


loop FOR ... NEXT. 
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I listing delle Figg.2.5 e 2.6 mostrano l’utilità dell’istruzione PRINT USING che 
permette di migliorare la presentazione dei risultati. 


100 N=100 

110 FORI = 1 TON 

120 FOR J = I+1 TON 

130 S=IMI+JMJ 

140 {= 

150 &a=K+1 

160 K2=KX*K 

170 IF K2 « S THEN 150 
180 IF K2 > S THEN 200 
190 FRINT USING 220,I,4,K 
200 NEXT 4 


210 NEXT I 
220 SERE 414 #44 
230 END 


Figura 2.5 
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Il listing di Fig.2.5 corrisponde ad un programma la cui uscita è molto ‘‘lineare‘* 
quindi un pò ingombrante. Per ridurre questa lunghezza si può preparare un’‘‘edi- 
zione‘ su tre colonne. Per questo si utilizza una variabile ausiliaria B che fà eseguire 
un passaggio alla linea seguente dopo la stampa successiva di tre soluzioni. Abbia- 
mo dunque modificato leggermente il flow-chart precedente nel seguente modo: 


PASSAGGIO 
ALLA LINEA 
SEGUENT 


Questa soluzione porta ai listati di Figg.2.6 e 2.7 
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100 N=100 
105 B-1 

110 FOR I = 1 TO N 

12 FOR J = I+1 TO N 


130 S=IXI+.]XJ 

140 {=J 

150 K=K+1 

160 K2=KXK 

170 IF K2 « S THEN 150 
180 IF K2 5» S THEN 200 
190 FRINT USING 220,I,4,K, 
195 B=E+1 

197 IF E «=3 THEN 200 
198 FRINT 199 

199 B=1 

200 NEXT 4 


210 NEXT I 
220 EEE +4 #44 I 
230 END 


Figura 2.6 


Esempio di esecuzione: 


3 4 5 19 12 13 I 6 8 10 I 
7 24 25 I 8 15 17 I 9 12 15 I 
9 40 41 110 24 26 111 60 61 1 
2 16 20 112 35 37 I 13 84 25 I 
14 48 50 1 15 20 25 I 15 136 391 
16 30 34 I 16 63 65 I 1f 24 30 I 
18 80 82 1 20 21 22 120) 48 52 I 
20 99 101 121 28 35 121 72 75 I 
24 2 40 I 24 45 51 I 24 70 741 
25 60 65 1 27 36 45 I 28 45 53 1 
2 96 100 130 40 50 130 72 78 I 
2 60 68 I 33 44 55 I 33 e) 65 I 
35 84 9% 136 48 60 I 36 77 85 I 
39 2 65 1 39 80 89 140 42 58 I 
40 75 85 1 40 96 104 142 56 701 
45 60 75 1 48 55 73 148 b4 80 I 
42 90 102 I 51 68 25 154 72 MI 
56 90 106 157 76 55 160 63 87 I 
60 80 100 160 9% 109 163 24 105 I 
65 72 97 I 66 ea 110 I 69 2 115 I 
2 96 120 175 100 125 If 24 116 I 


Figura 2.7 


2.2 NUMERI DI ARMSTRONG 


Si chiamano numeri di Armstrong, quei numeri che sono uguali alla somma dei 
cubi delle cifre che li compongono. Per esempio 153 è un numero di Armstrong 


poichè: 
153=1+5+3 
Scrivere un programma per ottenere tutti i numeri di Armstrong compresi tra l e 
2000. 
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Analisi di un problema: Per analizzare ogni numero, si dovrà scomporlo in 
ognuna delle cifre che lo compongono ed in seguito calcolare la somma dei cubi 


delle cifre. 
INIZIO 


R è uno dei valori cercati 


accumulo con gli altri cubi 
già calcolati 


quando K = 0 si sono ottenuti 
tutti i valori che compongono 


il numero I 
SI 


Figura 2.8 
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Per ottenere la cifra delle unità, è sufficiente calcolare il resto della divisione 
intera di questo numero per 10. Se I è il numero, calcoleremo come segue: 


Q= INT (1/10) 
R=1-10xQ R dà la cifra delle unità 


Per ottenere ora la cifra delle decine, sarà sufficiente ricominciare la stessa 
operazione partendo da Q. 


Q= INT (Q/10 ) 
R=Q-10-Q1 
Un modo comune di procedere può essere quello di continuare finchè il quozien- 
te della divisione è uguale a zero. Se si limita N = 2000, non si supereranno le 
quattro cifre. 
Piuttosto che calcolare Q1 poi Q2 e così di seguito, si può operare nel seguente 
modo: 


1. PorreK=I 
2. Calcolare Q = intero (K/10) 
R=K_-10xQ 


porre S = S + R3 

porre K = Q per la seguente iterazione 

se K > 0 ritornare in 2 altrimenti andare in 3 
3. Provare se S= I 


Questo ci porta al flow-chart di Fig. 2.8. 


Il listato di Fig. 2.9 corrisponde al flow-chart di Fig. 2.8. L'esempio di esecuzione 
mostra che i numeri di Armstrong non sono poi tanti. 


10 N=2000 

20 FRINT"NUMERI DI ARMSTRONG TRA 1 E 2000 " 
30. PRINT 

40 FOR I=1 TON 

50 5=0 

60 K=I 

70 Q=INT(K/10) 

80 R=K-10xQ 

90 S=5+RXRXR 

100 
110 
120 
125 
130 


140 Figura 2.9 


Esempio di esecuzione: 


NUMERT DI ARMSTRONG TRA 1 E 2000 
1 


153 
370 
371 
407 Figura 2.10 
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2.3 SCOMPOSIZIONE DI UNA FRAZIONE IN FRAZIONI EGIZIANE 


Una frazione quando ha per numeratore 1 è chiamata frazione EGIZIANA (1), 
esempio 1/3, 1/10. 


Si chiamano frazioni proprie quelle che hanno il numeratore più piccolo del 
denominatore. 


Problema: scomporre una frazione propria in una somma di frazioni Egiziane. 
Per far questo proponiamo di utilizzare l'algoritmo di FIBONACCI (2). 


Sia a/b la frazione da scomporre. Per determinare la prima frazione delia 
scomposizione, si prende la più grande frazione Egiziana il cui valore è inferiore ad 
a/10. Si sottrae questa frazione da a/b e si continua con questa procedura fintanto 
che il resto sia uguale a zero. 


Esempio: 
a=?2 b=3 a=- 2 
b 3 
la più grande frazione da prendere è } 
Qa_1-2_1-1 
b 2 3 2 6 
dunque 
3 2 6 


la scomposizione non è sempre così rapida. Ad esempio per 8/11 si trova 


Bode pd 


Il 2 5 55 110 


Esercizio: Costruire un programma che scomponga una frazione in frazioni 
Egiziane secondo l’algoritmo di Fibonacci. Porre particolare attenzione alla rappre- 
sentazione dei risultati. Si discuterà inoltre circa i limiti del programma e le 
precauzioni da prendere. 


dl pra frazioni erano utilizzate nell’antichità dagli Egiziani che non riuscivano ad utilizzare altri tipi 
i frazioni. 


(2) FIBONACCI: Leonardo da PISA conosciuto sotto il nome di Fibonacci nato verso il 1175 a Pisa, 
pubblicò questo algoritmo nel 1202. 
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Soluzione: Questo problema è in apparenza molto semplice, è necessario: 
- cercare la più grande frazione Egiziana inferiore ad a/b 
- calcolare la frazione restante. 

Se si parte da C = intero di B/A, si avrà, 


dd 
C B 


ma C sarà molto vicino al denominatore cercato. E’ sufficiente dunque porre C= 
C+ 1 fino a che 1/c < A/B e si è allora trovata la frazione cercata. 
La frazione restante si ottiene da: 


drei 
B C 


che è dato come 


AC — B= numeratore 
BC = denominatore 


Partendo da questa analisi si può disegnare un primo flow-chart. Ma si può già 
fare un’importante considerazione: i calcoli devono essere effettuati unicamente 
con numeri interi che possono diventare molto grandi. E’ noto che la rappresenta- 
zione dei numeri in BASIC dà una precisione limitata. E’ necessario dunque 
prevedere un test per assicurarsi che non sia stata superata la precisione del 
calcolatore utilizzato. 

Questo test dovrà essere fatto a priori dopo ogni moltiplicazione AxC e BxC, ma 
quando B > A, è sufficiente fare un test sulla seconda moltiplicazione. Questo porta 
al flow-chart delle Figg. 2.11 e 2.12. Il test di arresto è B< P. 


NOTA: Come nei problemi precedenti, si possono utilizzare le variabili intere se il 
computer utilizzato lo permette. Tuttavia per i microcalcolatori questo porta, 
generalmente, ad una capacità inferiore. E’ dunque preferibile lavorare con delle 
variabili normali. 


Il programma è stato realizzato in due parti: 
- un programma principale che effettua l’I/O e qualche test, 
- un sottoprogramma che a ciascuna ‘iterazione‘‘ cerca la più grande 
frazione Egiziana possibile, e la frazione restante per la tappa 
successiva. 


Spezzando il programma in due parti, si aumenta un pò il numero delle istruzio- 
ni, ma la sua scrittura è resa più semplice. 


Nel flow-chart di Fig. 2.11, si fà appello alla variabile L che assume solo 2 valori: 
0 e 1. Inizialmente questa variabile è posta a zero e permette, grazie al test, di non far 
precedere le prime frazioni trovate dal segno +. Poi, assume il valore 1 e ogni 
frazione seguente è preceduta dal segno +. Si ottiene così una stampa molto 
elegante‘ sebbene tutto sia sistemato su una sola linea. 
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INIZIO 
LEGGERE 
PRECISIONE 
LEGGERE 
AEB 


STAMPA LA 
FRAZIONE EGIZIANA 
OTTENUTA 1/C 


CODICE 
ERRORE 


Figura 2.11 


INIZIO 


C= INTERO A/B 
A1= A/B 


RITORNO 


RITORNO 


Figura 2.12 


Commenti al programma: 
- Per un calcolatore la precisione dei dati è costante. Si può dunque: 


@ sia indicare al programma il numero più grande possibile nel programma per 
mezzo di un’istruzione di assegnazione oppure delle istruzioni READ e DA- 
TA. Quest'ultima soluzione è quella adottata, 


@ sia domandare all’utilizzatore di indicare il numero più grande ammissibile 
durante l’esecuzione. Questa soluzione è meno pratica da realizzare. 


- Per fermare il programma è sufficiente introdurre due numeri A e Bin modo che A 
sia maggiore di B. 


Suggerimenti: Partendo da questo programma se ne può scrivere un’altra versione 
di tipo ‘““conversazionale‘* che permetterà all’utilizzatore stesso di provare un 
diversa scomposizione, senza che sia obbligato ad effettuare i calcoli. Il programma 
s’incarica di indicare la frazione restante. 
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100 FRINT "SCOMPOSIZIONE IN FRAZIONI EGIZIANE" 


105 READ FP 

110 PRINT 

120 PRINT 

130 FRINT "IMFOSTARE NUMERATORE, DENOMINATORE "3 
140 INFUT A,E 

150 IF A > = E THEN 800 

160 L= 0 

170 PRINT ? FRINT 

180 PRINT "FRAZIONE "3A;"/"B;" = o "} 


190 IF A = 0 THEN 110 
200 GOSUE 500 

210 IF L = 0 THEN 230 
220 PRINT" + "3; 

230 FRINT "1/"3}C; 

235 L= 1 

240 IF E « F THEN 190 
300 FRINT " SUFERO DI CAFACITA”" 
310 GOTO 110 

500 IF A > 1 THEN 600 
510 C = B 

520 A= 0 

530 RETURN 

600 C = INT (E / A) 
605 AL =A/ EB 

610 IF 1/C < Al THEN 640 
620C=C+1 

630 GOTO 610 

6AL ASAXC - B 

650 E= BxC 

670 RETURN 


700 DATA 32767 «=———_—_—___ massimo numero intero ammissibile 
800 END 
Figura 2.13 

Esempio di esecuzione: 
SCOMPOSIZIONE IN FRAZIONI EGIZIANE 
IMPOSTARE NUMERATORE, DENOMINATORE 22,3 
FRAZIONE 2/3 = 1/2 + 1/6 
TMPOSTARE NUMERATORE, DENOMINATORE  ?8,11 
FRAZIONE 8/11 = 1/2 + 1/5 + /37 + 1/4070 
TMFOSTARE NUMERATORE, DENOMINATORE 23,7 
FRAZIONE 3/7 = 1/3 + 1/11 + 1/231 
TMFOSTARE NUMERATORE, DENOMINATORE  ?7,13 


FRAZIONE 7/12 = 1/2 + 1/26 
TMPOSTARE NUMERATORE, DENOMINATORE ‘213,29 


FRAZIONE 13/29 e 1/9 + 1/9 + 1/262 SUFERO DI CAFACITA” 
TMFOSTARE NUMERATORE, DENOMINATORE  ?P11,2 

FRAZIONE 11/29 = 1/3 + 1/22 + 1/1914 

TMFOSTARE NUMEFATORE, DENOMINATORE 23,2 


Figura 2.14 
2.4 I NUMERI PRIMI 


Un primo metodo per trovare i numeri primi consiste nel cercare, partendo dal 
numero 3, i numeri dispari che accettano come divisori solo se stessi o la cifra 1. 
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Studieremo in seguito un metodo più raffinato. 


2.4.1 Primo metodo: Scrivere un semplice programma che stampi N numeri primi. 


Per l’esecuzione si prenderà N tra 10 e 60. Studiare poi come migliorare la presenta- 
zione. 


2.4.2 Soluzione: L’ossatura generale del programma corrisponde al flow-chart di 
Fig. 2.15. 


INIZIO 
NUMERI PRIMI 
1 2 E 3 


STAMPARE I 


Figura 2.15 
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— stampare i numeri 1, 2 e 3 

— poi cercare, in seguito, N-3 altri numeri primi facendo avanzare I di 
due passi per volta; di due passi perchè i numeri primi sono tutti 
dispari. 


Per determinare se I è un numero primo, sarà necessario eseguire test successivi 
utilizzando numeri dispari fino ad ottenere uno dei seguenti casi: 


— si ottiene come resto 0, questo significa che I non è un numero primo 
— si ottiene come resto un numero diverso da zero ed un quoziente 
inferiore o uguale al divisore, in questo caso I è un numero primo. 


Per rispondere alla domanda “I è un numero primo ?“, è necessario eseguire una 
parte del flow-chart di Fig. 2.16. Si possono ora tracciare i dettagli del flow-chart e 
soprattutto si può scrivere il programma. 


Q = INTERO (1/K) 
R=1-QxK 


| NON È PRIMO 


K=K+2 


| È PRIMO 


Figura 2.16 
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100 N = 60 
110 F 


‘RE 1° PRINT" di 


INT "I FEIMI "3N3" NUMERI PRIMI SONO? " 
PRINT 1,2,2 
2 
TON - 3 
(IZ K) 
x K 


THEN GOTO 150 
SO THEN  GOTO 230 

Kos Ko+ 2 

GOTO 170 


250 PRINT I, 
240 NEXT 4 
250 END 


Figura 2.17 


Esempio di esecuzione: 


I PRIMI 60 NUMERI FRIMI SONO? 
1 Q 3 


bi) 7 11 13 17 

19 29 chi 37 

41 47 53 59 

61 ta! 73 79 

e3 97 101 103 
107 113 127 131 
137 149 151 157 
163 173 179 181 
191 197 199 211 
223 229 233 239 
DAL 257 263 269 
271 

Figura 2.18 


2.4.3 Secondo metodo: Cominciando dal numero 5 incluso, tutti i numeri primi sono 
dati come 6 n + 1, con n intero. Inoltre è sufficiente cercare i divisori tra i numeri 
primi già trovati. Scrivere un programma che tenga conto di queste due note. 


2.4.4 Soluzione: Per ricercare gli eventuali divisori fra i numeri primi già trovati è 
necessario poterli memorizzare, questo rende necessario l’utilizzo di una tabella la 
cui dimensione limiterà il numero massimo dei numeri primi che si potranno 
cercare. Inoltre i numeri dati come 6 n + 1 sono indivisibili per 2 e per 3, sarà quindi 
sufficiente cercare dei divisori a partire dal numero 5. 


Scomponiamo il lavoro in due parti: 


— un programma principale inizializza la tabella T, aggiorna la variabile A e 
chiama un sottoprogramma. 
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— il sottoprogramma prova se la variabile A rappresenta un numero primo, se sila 
pone nella tabella T. Quando la tabella T è riempita il contenuto viene stampato. 


Si giunge così ai flow-chart di Figg. 2.19 e 2.20. 


Esempio di esecuzione: 


2.5 SCOMPOSIZIONE IN FATTORI PRIMI 


Scomporre un numero in fattori primi consiste nel cercare tutti i divisori primi di 
questo numero. 


2.5.1 Metodo elementare: Si cercheranno i divisori a partire dal numero 2. Dal 
momento in cui si trova un divisore lo si stampa. Quando un divisore non è più 
conveniente, si passa al divisore successivo. 


T(1)=1 1(2)=2 
T(3)=3 T(4)=5 
N=100 1=3 
A=5 


Figura 2.19 
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100 
110 
150 
160 
170 
180 
190 
200 
210 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 


INIZIO 


R=A-UX INTERO(Î) 


Figura 2.20 


DIM T(200) 
T(1) = 13102) = 23703) = 3704) = 5 
Bo = SII = 3 
N = 200 

GOSUE 500 
B= B+ 2 

GOSUE 500 

E = BE + 4 

GOTO 170 

FOR J = 4 T01I 

U = T(J) 

IFU XU x BE THEN 560 
{= B = INT (E / VU) x U 
IF R = 0 THEN RETURN 


IsI+1 


IF I N THEN RETURN 
PRINT "IL SEGUENTE LISTATO CONTIENE "iN" NUMERI PRIMI?" 
FRINT 
K_= 0 

FORI = 1 TON 

PRINT TAEC 7 * K)IsTC(I); 
Ka K+ 1 

IFK 3 10 THEN 670 

Kos 0% PRINT 

NEXT I 

END 


Figura 2.21 
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Esempio di esecuzione: 


IL LISTATO CONTIENE 200 NUMERI PRIMI 


1 2 3 ij 7 11 13 17 19 23 
29 31 37 41 43 47 53 59 61 67 
71 73 79 83 29 97 101 103 107 109 
113 127 131 137 139 149 151 157 163 167 
173 179 181 191 193 197 199 eil 223 227 


29 233 239 241 251 257 263 269 271 277 
281 283 293 307 311 313 317 331 337 347 
349 353 359 367 373 379 383 389 397 401 
409 419 421 431 433 439 443 449 457 461 
463 467 479 487 491 499 503 509 521 523 
541 547 557 563 569 571 577 587 593 599 
601 607 613 617 619 631 641 643 647 653 
659 661 673 677 683 691 701 709 719 727 
733 739 743 751 757 761 769 773 787 797 
809 811 821 823 827 829 839 853 857 859 
863 877 881 883 887 907 911 919 929 937 
941 947 953 967 971 977 983 991 997 1009 
1013 1019 1021 1031 1033 1039 1049 1051 1061 1063 
1069 1087 1091 1093 1097 1103 1109 1117 1123 1129 
1151 1153 1163 1171 1181 1187 1193 1201 1213 1217 


Figura 2.22 


Se, ad un certo punto, si trova un quoziente più piccolo del divisore: 

— se il dividendo è il numero dato, questo è il numero primo, 

— se il dividendo è inferiore, questo allora è un numero primo che divide il 
numero. 

Costruire un programma che effettua questa scomposizione, che poi richieda un 


altro numero finchè non gli venga dato un numero negativo oppure zero, il che 
costituisce un criterio di arresto. 


2.5.2 Soluzione: La struttura generale del programma è data dal flow-chart di Fig. 


2.23. Bisogna ora dettagliare la parte ‘“scomposizione*“. 


Per effettuare questa scomposizione si può utilizzare il seguente algoritmo: 
1. Ricopiare N in Ni 


2. Per I che assume successivamente i valori 2, 3, 4, 5, etc. 
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VISUALIZZAZIONE 
DI UN TITOLO 
STAMPARE 
‘IMPOSTARE N’ 
LEGGERE N 


SCOMPOSIZIONE 


Figura 2.23 


2.a Provare se I è un divisore di N, 
sia Q il quoziente, 
se I è un divisore, stampare I, porre N=Q e ritornare in 2.a, 
se I non è un divisore andare in 2.b. 


2.b Se il quoziente Q è più grande di I, incrementare I ritornando in 2, 
se il quoziente Q è « I allora: 


se N=1 è terminata la scomposizione 

se N=N allora N è primo, 

se NZN: allora N è un divisore e bisogna stamparlo alla fine della 
scomposizione. 


Questo metodo porta al flow-chart di Fig. 2.24 partendo dal quale si è scritto, 
senza difficoltà, il programma in Fig. 2.25. 
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42 


sottoprogramma 
di scomposizione 


Q= INTERO (N/1) 
R=N— QX1 


STAMPARE | 


SI 
STAMPARE N 
È PRIMO STAMPARE 


FINE 
SCOMPOSIZIONE 


Figura 2.24 


100 REM SCOMPOSIZIONE IN FATTORI FRIMI 


110 REM 

120 PRINT "SCOMPOSIZIONE IN FATTORI PRIMI" 
130 FRINT ? PRINT 

140 INFUT "IMFOSTARE IL NUMERO DA SCOMFORRE "N 
150 IF N < = 0 THEN STOF 

160 N1 = N 

170I=1 

180 I=I+1 

200 Q so INT (N / I) 

210 REN- QX I 

220 IFR > 0 THEN 290 


QQ0N= Q 
DAN FFINT I" "} 
250 GOTO 200 


290 IFQO> 
300 IF N= 
310 IFON « 2 Ni THEN 340 


320 PRINT TAEC 5)p"E” PRIMO." 
220 GNOTO 350 
240 PRINT N 


350 PRINT 
360 GOTO 130 
3Z8, END 


Figura 2.25 


Esempio di esecuzione: 
SCOMPOSIZIONE IN FATTORI PRIMI 
IMPOSTARE IL NUMERO DA SCOMPORRE 12 


2:23 


TMPOSTARE IL NUMERO DA SCOMPORRE 38 
2° 19 


IMPOSTARE IL NUMERO DA SCOMPORRE 262144 
2 


» . , , 5 
LE LA RR PERI? 


IMPOSTARE IL NUMERO DA SCOMFORRE 65784 


2223 27A1 


IMFOSTARE IL NUMERO DA SCOMPFORRE 1217 
E’ PRIMO, 


IMFOSTARE IL NUMERO DA SCOMPORRE 35427 
377 2AI 


IMPOSTARE IL NUMERO DA SCOMPORRE 0 
EREAK IN 150 


Figura 2.26 


2.5.3 Metodo per una elegante presentazione: La presentazione dei risultati ottenuti 
con il metodo precedente non è certo molto gradevole. Migliorarla per ottenere il 
seguente risultato: 


SCOMFOSIZIONE IN FATTORI FRIMI 


IMFOSTARE IL NUMERO DA SCOMFORRE 65784 
E” DIVISIEILE PER 2 3 VOLTE 
E‘ DIVISIBILE PER 3 1 VOLTE 
E‘ DIVISIBILE FER 2741 1 VOLTA 
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IMFOSTARE IL NUMERO DA SCOMFORRE 1217 


E‘ PRIMO. 
IMFOSTARE IL NUMERO DA SCOMPORRE 35427 
E/ DIVISIBILE FER 3 1 VOLTE 
E‘ DIVISIBILE PER 7 2 VOLTE 
E‘ DIVISIBILE PER 241 1 VOLTA 
IMFOSTARE IL NUMERO DA SCOMFORRE 262144 
E” DIVISIBILE FER 2 18 VOLTE 
IMFOSTARE IL NUMERO DA SCOMFORRE 19 
E‘ FRIMO, 
IMFOSTARE IL SCOMFORRE 14 
E” DIVIS 1 VOLTE 
E‘ DIVI! 1 VOLTA 
IMFOSTARE IL NUMERO DA SCOMFORRE 0 
EREAK IN 160 
Figura 2.27 


Si chiede di costruire il programma corrispondente. 


2.5.4 Soluzione: Si può partire dal flow-chart di Fig. 2.24 che dovrà essere modifica- 
to per stampare solo quando si termina con un divisore. La parte incorniciata del 
flow-chart precedente è dunque sostituita dal flow-chart che segue: 


STAMPARE È DIVISIBILE 
PERI J VOLTE 


Figura 2.28 
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E° ora sufficiente modificare il programma precedente per giungere alla nuova 
stesura. Oltre alle modifiche date per la parte di flow-chart, non bisognerà dimenti- 
care di modificare le istruzioni di stampa. Questo porta al listato di Fig. 2.29. 


240 
250 
260 
270 
280 
290 
200 
310 
220 


N 


INPUT 


SCOMPOSIZIONE 


SCOMPOSIZIONE 
$_ PRINT 


IN FATTORI FEIMI 


IN FATTORI FEIMI" 


IMPOSTARE IL. NUMERO DA SCOMPORRE "N 
Ls N 
IF Nè O THEN STOP 
he 
I+1 
0 
TINTOCN Z IM 
= Noe QI 
IFR£ 2 0 THEN 260 
#0 
Ja J +1 
cOTO 200 
IF 4 = 0 THEN 290 
PRINT" E DIVISIBILE FER "5135" "3 ThE: 10);Jj TAEC 
GOTO 180 
IFG > T THEN 180 
IF N = 1 THEN 250 
IFON NI THEN 340 
PRINT TAEG 5)p"E/ PRIMO," 
GOTO 2% 
PRINT" E” DIVISIBILE FER "INI TEC 105" 1"f TA[C 18);" 
PRINT 
GOTO 14N 
END 
Figura 2.29 


2.6 CONVERSIONE BASE 10 IN UN’ALTRA BASE 


La rappresentazione di un numero in differenti sistemi di numerazione (base 10, 
base 8, base 2, etc.) costituisce per “l’uomo della strada‘, un esercizio matematico 
senza concreta utilità. Invece, per gli “‘informatici‘*, questo lavoro presenta un’utili- 
tà reale soprattutto quando si programma in linguaggio assembler. Per questa 
ragione sarà presentato questo esercizio. Il principio della conversione è il seguente: 


— effettuare delle divisioni successive per la nuova base fino a che nonsi ottenga 
un quoziente inferiore al valore della nuova base; 


i 
Esempio: conversione di 83 scritto in base 10 nella base 8: 


8 3 


0 3 


10 


2 


Risultato 1 2 3 
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— la cifra delle unità corrisponde al primo resto. La cifra seguente corrisponde 
al resto seguente e così via. L’ultima cifra è l’ultimo quoziente. 


Così 83 dà in base 8 123, 
83 dà in base 7 146. 


2.6.1 Conversione in base minore di dieci: Scrivere un programma che stampi una 
tabella di conversione tra due numeri D e F dati dall’utilizzatore. La conversione si 
effettua da base 10 verso un’altra base inferiore a 10. 


Soluzione: La costruzione di questo programma ha qualcosa di familiare con i 
programmi precedenti: 


— utilizzo della ‘‘divisione intera“, 
— calcolo del resto, 
— inserimento in una tabella. 


Per meglio distinguere la struttura generale dell’algoritmo propriamente detto, si 
scomporrà il programma in due parti: il programma principale che effettua l’I/O e 
un sottoprogramma corrispondente all’algoritmo di cambiamento di base. 


La costruzione del flow-chart di Fig. 2.30 non presenta difficoltà. Mentre il 
flow-chart di Fig. 2.31 necessita di alcune spiegazioni. 


— Quando si comincia ad effettuare la conversione, non si sa a priori quante 
cifre significative avrà il nuovo numero: sarà dunque obbligatorio inserire le 
cifre nell’ordine in cui si ottengono. 


— Il metodo proposto fornisce prima le cifre relative alle unità, poi quelle delle 
decine (o più esettamente quelle corrispondenti a un esponente 1) e così di 
seguito. Per inserirle nella tabella A si procede partendo da J=° 1. 


A ogni cifra trovata A(J) =R 
J = J+ 1 (perla successiva cifra da inserire) 


Per l’ultima cifra si pone A(J)= Q 
Si può ora costruire il flow-chart di Fig. 2.31 del sottoprogramma di conversione 
e in seguito scrivere il programma completo. Per la stampa del numero desiderato, è 


necessario procedere nell’ordine inverso rispetto a quello utilizzato per ottenerlo; 
per esempio se il numero convertito è 127, si è ottenuto: 


in questo caso L = 3 
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Per stampare con ordine si scriverà: 


FOR I=L TO 1 STEP 1 


PRINT A(I); per restare sulla stessa linea 
NEXT I 
PRINT per passare alla linea seguente 


INIZIO 
LEGGERE 
B, D, F 


INIZIO 


S/P 
CONVERSIONE 
STAMPARE | 
E NUMERI CONVERTITI 


Figura 2.30 Figura 2.31 
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95 DIM AC15) LA NUOVA BASE ?16 


100 FRINT "LA NUOVA BASE "} PRIMO E ULTIMO NUMERO DA 
110 INFUT E CONVERTIRE ?16380,16390 
120 PRINT "PRIMO E ULTIMO NUMERO DA CONVERTIRE "3; 16380 3151512 

130 INFUT F,L 16381 3151513 

140 FORI =F TOL 16382 3151514 

150 PRINT 16383 3151515 

160 GOSUE 1500 16384 4000 

170 REM STAMFA LA TABELLA 16385 4001 

180 FRINT I; TAEC 6); 16386 4002 

190 FOR D = J TO 1 STEP - 1 16387 4003 

200 FRINT ACD); 16388 4004 

210 NEXT D 16389 4005 

220 NEXT I 16390 4006 

230 STOF EREAK IN 230 

1500 I1 = 1 1 

1510 J= 1 


1520 Q = INT C(I1 / E) 

1530 R= I1-Qx BE 

1535 I1 = Q 

1540 AC) = 1 

1545 J=J+1 

1550 IF Q®> = B THEN 1520 
1560 AC) = Q 

1570 RETURN 


1580 END 

Figura 2.32 
Esempio di esecuzioni. 
LA NUOVA BASE ?Z LA NUOVA BASE ?8 
PRIMO E ULTIMO NUMERO DA CONVERTIRE ?260,280 FRIMO E ULTIMO NUMERO DA 
260 100000100 CONVERTIRE 270,90 
261 100000101 70 106 
262 100000110 71 107 
263 100000111 72 110 
264 100001000 73 111 
265 100001001 74 112 
266 100001010 75 113 
267 100001011 76 114 
268 100001100 77 115 
269 100001101 78 116 
270 100001110 79 117 
271 100001111 80 120 
272 100010000 81 121 
273 100010001 82 122 
274 100010010 83 123 
275 100010011 84 124 
276 100010100 85 25 
277 100010101 86 126 
278 100010110 87 127 
279 100010111 88 130 
280 100011000 89 131 
EREAK IN 230 90 132 
J EREAK IN 230 

J 
Figura 2.33 


2.6.2 Conversione in base maggiore di 10: Estendere il programma al fine di permet- 
tere la stampa di una tabella di conversione per una base maggiore di 10. In questo 
caso si rappresenterà la cifra 10 con la lettera A, 11 conla lettera B e così di seguito. 
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Soluzione: Si utilizzeranno le funzioni relative alle stringhe di caratteri; si può per 
esempio dare una stringa B$ tale che: 


B$ = “0123456789ABCDEFGHIJKL“ 


A partire dall’A(L) precedente, è sufficiente estrarre il carattere in posizione 
A(L)+1 nella stringa (la cifra 0 corrisponde al primo carattere di B$). 
Utilizzando, in funzione del sistema disponibile, SUBSTR o MID$, si scriverà: 


00200 


PRINT SUBSTR (B$ , A (L)+1 ,1); 


per stampare il carattere appropriato. 


Questo porta al programma della Fig. 2.34. 


10 
50 
90 E 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
1480 
1500 
1510 
1520 
1530 
1535 
1540 
1545 
1550 
1560 
1570 
1580 
1 


REM FR 
DIM ACI 
$ = "01 
FRINT 
INFUT 
FRINT 
INFUT 
FOR I 
FRINT 
GOSUE: 
REM S 
PRINT 
FOR D 
PRINT 
NEXT D 
NEXT I 
STOF 
REM 
I1=1I 


Ii 


END 


OGRAMMA DI CONVERSIONE 

5) 

234567 89AECDEFGHIJKLMN" 

"LA NUOVA BASE "} 

E 

"FRIMO E ULTIMO NUMERO DA CONVERTIRE "3 


1500 

TAMFA LA TABELLA 

I} TAEC 6)} 

= J TO 1 STEP - 1 
MID$ (E$,A(D) + 1,1); 


SUEROUTINE DI CONVERSIONE 


-Qx E 


= B THEN 1520 


Figura 2.34 
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LA NUOVA EASE ?8 
FRIMO E ULTIMO NUMERO DA CONVERTIRE 275,90 


5 113 
76 114 

77 115 

78 116 

LA NUOVA BASE ?2 9, 117 
PRIMO E ULTIMO NUMERO 80 1£0 
DA CONVERTIRE ?260,280 - ct 
260 100000100 83 123 
261 100000101 84 124 
262 100000110 85 125 
100000111 86 126 
100001000 87 127 
100001001 8 130 
100001010 89 131 
10000L011 90 132 


100001100 
100001101 
100001110 
100001111 
100010000 


LA NUOVA EASE 216 
FRIMO E ULTIMO NUMERO DA CONVERTIRE ?1023,1035 


> ei 
100010001 Des 3EF 
1024 400 
100010010 ea 
: 1025 401 
100010011 1026 p 
ARA 26 402 
100010100 [e ; 
1027 403 
100010101 È 
PRAIA 1028 404 
100010110 ; È: 
de 1029 405 
100010111 4 
100011000 03 306 
, 1031 407 
1032 408 
1033 409 
1034 40A 
1035 40E 
EREAK IN 230 
1 
Figura 2.35 


Casi particolari: 


Alcuni calcolatori non dispongono delle funzioni SUBSTR o MID$ ma esiste 
comunque un’altra possibilità: dopo aver lanciato il programma dichiarare la 
lunghezza massima della variabile stringa di caratteri B$ si può estrarre una 
sottostringa scrivendo l’ espressione: 


B$ (I,J) 
nella quale I rappresenta la posizione del primo carattere della sottostringa e J la 
posizione dell’ultimo. E’ il caso, per esempio, del calcolatore LOGABAX LX 500. 


Si scriverà dunque in questo caso: 


AI=A(L)+1 
PRINT B$ (AI, A1) 


per stampare solo un carattere, da cui il programma di Fig. 2.36. 


50 


10 REM CONVERSIONE DI EASE 
20 REM AUTORE J.F., LOMOITIER 
300 DIM AC15),E$(30) 


40 E$ = "01234567 89AECDEFGHIJKLMNOFOR" 

100 INFUT "IMFOSTARE LA NUOVA BASE ",E 

110  FRINT 

120. INFUT "IMFOSTARE L'INIZIO E LA FINE DELLA TABELLA ",D,F 
125 FRINT 


130 FRINT "EASE 10", TAEC 10), "BASE ",E 
140 FORI = DTOF 

150 PRINT 

160 GOSUE 1500 

170 REM STAMFA 

180 FRINT I, TAEC 11), 


190 FOR L = 4 TO 1 STEP - 1 
195 A1 = ACL) + 1 TMENETARE LA NUOVA EASE 16 
200 FRINT E$(A1,A1), TWFNSTARE L'INTZIO E LA FINE 
210 NEXT L DELLA TABELLA 12,20 

220 NEXT I 

230 STOF BASE 10 BASE 16 

1500 IL = I 

1510 J= 1 12 oc 

1520 Q = INT (IL / E) 13 OD 

1530 R = IL - QX E 14 DE 

1535 Il = Q 15 0F 

1540 AJ) = R 16 10 

1545 J=J+1 17 ii 

1550 IF Q< = E THEN 1520 18 12 

1560 AJ) = Q 19 13 

1570 RETURN 20 14 

1580 END 


Figura 2.36 


2.7 CONCLUSIONE 


Questi esercizi di differenti difficoltà mostrano l’utilità della costruzione di un 
flow-chart per gradi, partendo dalla struttura generale del programma per poi 
procedere dettagliatamente fino a sentirsi pronti per la scrittura dello stesso. 


Per risolvere un problema, la soluzione in generale non è una sola sia sul piano 
del metodo che su quello della programmazione. I programmi presentati in questo 
libro non sono tra i migliori dal punto di vista dell’efficienza, ma sono in generale i 
più facili da comprendere e molto prossimi al flow-chart. Man mano che si 
acquisisce esperienza, si può ridurre il tempo destinato alla costruzione del flow- 
chart, contentandosi di un flow-chart di principio dal quale ricavare la scrittura del 
programma. 
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CAPITOLO 3 


SEMPLICI ESERCIZI DI GEOMETRIA 


3.0 INTRODUZIONE 


La geometria classica non si presta ad applicazioni numeriche, al contrario della 
«geometria analitica‘ elementare. Gli esercizi presentati in questo libro sono di 
facile comprensione e mostrano le capacità del microcalcolatore. Sono stati scelti 
per il loro interesse pratico e per la loro semplicità: il flow-chart si costruisce 
facilmente e la programmazione è ugualmente facile. I calcoli, necessitando di una 
certa precisione, risultano complessi se eseguiti a mano mentre sono eseguiti molto 
rapidamente con l’ausilio del calcolatore. 


Partendo da questi esempi il lettore accorto ne potrà realizzare altri più complessi 
o più adatti alle proprie esigenze. 


3.1 PERIMETRO E AREA DI UN TRIANGOLO 


Per calcolare la superficie di un triangolo qualunque, si può misurare la 
lunghezza di ogni lato e ottenere così l’area del triangolo applicando la formula di 
Erone. 


A=V s(s- a)(s-b)(s-c) 
con 


s= atb+e€ essendo a, be c le lunghezze dei lati 


Problema: Scrivere un programma che a partire da a, b, c calcoli il perimetro e 
l’area del triangolo. 


Soluzione: Partendo da a, b, c il calcolo non presenta difficoltà particolari: si calcola 
P=a+b+cche viene stampato, in seguito si applica la formula di Erone avendo 
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precedentemente calcolato il semiperimetro. Nel programma di Fig. 3.1 P rappre- 
senta prima il perimetro poi il semiperimetro. 


" "IMFOSTARE LE TRE LUNGHEZZE DEI LATI DEL TRIANGOLO"; 
JT A,E,C 


"PERIMETRO ="3F 


70 FRINT "S&S 
80 END 


Figura 3.1 


Commenti: 


Linea 10. Il punto e virgola, come la virgola posta alla fine di una linea, impedisce 
il passaggio alla linea successiva; questo permette dunque l’impostazione di dati 
sulla stessa linea. 


Linea 40. PRINT “PERIMETRO” =;P. Il punto e virgola fa si che la stampa del 
valore numerico della variabile P sia posta vicino al segno =. 


Linea 45. La sola istruzione PRINT provoca un salto di linea. Il suo utilizzo 
permette dunque la “buona“ presentazione dei risultati. 


Linea 50. Quando il valore del perimetro è stampato, non è più necessario questo 
valore, si può dunque utilizzare la stessa variabile per “inserire‘ il valore del 
semiperimetro, valore che sarà necessario in seguito. 


Critica del programma: 


Se vengono introdotte lunghezze di lati corrispondenti a un triangolo impossibi- 
le, per esempio 10, 20, 40, il calcolatore non avrà alcuna possibilità di rilevare 
questo errore e di segnalarlo. La macchina tenterà di estrarre la radice quadrata di 
un numero negativo, tutto ciò porta, generalmente, a spiacevoli conseguenze. 


10 PRINT "IMPOSTARE LE TRE LUNGHEZZE DEI LATI DEL TRIANGOLO" } 
20 INFUT A,LE,C 

30 F +0 

40 PRINT "FERIMETRO ="3F 
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45 PRINT 
50 P=, 
FE 

56 
57 
58 STOP 

60 S=SQR(PX(F-A)X(F-E)X(P-C)) 
70 FRINT "SUFERFICIE ="35 

80 END 


Figura 3.2 


Per rimediare al problema sarà necessario aggiungere un controllo di validità dei 
dati: la lunghezza del lato maggiore non deve superare la somma delle lunghezze 
degli altri due. Si può dunque aggiungere un test per verificare questa condizione 
oppure provare che 

(S_—- a)(s_ b)(s— ce) è positivo 


Questa prova porta al programma di Fig. 3.2 


3.2 DETERMINAZIONE DI UN CERCHIO PASSANTE PER TRE PUNTI 

Partendo dalle coordinate cartesiane di tre punti Mi, M: e M: determinare il 
cerchio che passa per questi tre punti, cioè trovare le coordinate del centro e la 
lunghezza del raggio. 


Analisi matematica: 


Siano (x1, y1), (x, y2), (x3, y3) le rispettive coordinate di M1, M2 e M+. La pendenza 
dell’angolo della linea che congiunge Mi: a M? è data dalla: 


Va 


XX 
quindi la pendenza di una perpendicolare è data da — 


l'equazione della mediana del segmento Mi M? è 


Vit ya XX ( 23%) 
= —— * — I 
di 2 Va7 Vi 
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dalla stessa equazione della mediana del segmento M: M: è 


y= pit yi xx (- =® ) 
2 V3- VI 


Queste due equazioni si possono scrivere sotto la forma: 


v= K:x+H: 

y= K3x + Hs 
Con K:- X2- XI k3=- X3- XI 
V2- V! V3- VI 


Hael!t}+ x:- x; Hj;= Vty:4 x- x, 
2 2(y2-y) 2 2(y2-y) 


Risolvendo il sistema lineare si ottengono le coordinate (Xo, Yo) del centro I: 


Xo= H:- H: 
K:- K; 


d;= K; H:- K:H; 
o= BID? A211: 
K3- K? 


Partendo dalle coordinate (Xo, Yo) del centro I si ottiene la lunghezza R del rag- 
gio: 


=V(—- xd +41 vò? 


Flow chart: 


La costruzione è realizzata senza difficoltà poichè è sufficiente seguire l’anda- 
mento del calcolo. 
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INIZIO 


CALCOLO DI K2 


CALCOLO DI 
Hz — Ha 
S Ki 


CALCOLO DEL RAGGIO 
R=v(x, — xo)? + (Yi — Yo)? 


VISUALIZZAZIONE 


Figura 3.3 


100 PRINT "DETERMINAZIONE DI UN CEKCHIO PASSANTE PER TRE PUNTI 


ORDINATE DEI TRE FUNTI DEVONO LSSERE FIAZZATE IM 
UNA TSTRUZIONE ‘DATA? FRIMA DELL'ESECUZIONE., 


160 IF D=0 THEN 230 

250 GRAY LAY (KQAXQEXIANIA)ZOY2=-Y1)) 
DL SX CY1+Y 34 XXI XNIRXI)ZOYB-Y13) 
H2-H3)/D 


210 RESQRCOXI-X0)YK2+(Y1-Y0)Mx2) 

220 PRINT "X0= "}X0j" YO= "3Y0}" R= "3R 

225 STOF 

230 FRINT "FUNTI ALLINEATI NESSUNA SOLUZIONE 
240 DATA 2,-1,0,1,2,3 

250 END 


Figura 3.4 


Esempio di esecuzione: 


DETERMINAZIONE DI UN CERCHIO FASSANTE FER TRE FUNTI 


Xo= 2 Yo= 1 Re 2 


Figura 3.5 


3.3 PERIMETRO DI UN POLIGONO 


Spesso i campi assumono forme geometriche corrispondenti a poligoni qualsiasi 
(per esempio quadrilateri). Si desidera conoscere la lunghezza del perimetro al fine 
di calcolare il prezzo di un ipotetico appezzamento. 


Problema: Si sono potute rilevare le coordinate cartesiane di ciascun vertice del 
campo o di una figura simile avente forma poligonale. Costruire un programma che 
calcola la lunghezza dell’appezzamento. 


Soluzione: Il lavoro consiste nel calcolare la lunghezza di ogni lato e nel fare la loro 
somma. Siano X(I), Y(I) le coordinate del vertice I. Il lato con l’estremo I eIl+1 
avrà come lunghezza: 


VOAa+bD—-YM)Y+(X(A+1)— X (MD). 


E’ sufficiente leggere prima il numero dei vertici (che è uguale al numero dei lati) 
poi leggere l’insieme X(I), Y(I) e procedere con il calcolo. 


Non bisognerà, nel frattempo dimenticare, che l’ultimo lato ha per estremità i 
vertici N e 1. Da ciò il flow-chart di Fig. 3.6. 
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CALCOLO DELLA LUNGHEZZA (1) 
DEL LATO DI VERTICE (I, 1+1) 


P-P+L(I) 
I=1+1 


CALCOLO DELLA LUNGHEZZA 
L(N) DEL LATO DEL VERTICE (N,1) 
P= P+L(N) 


VISUALIZZAZIONE 
DEI DATI 


Figura 3.6 


Il programma di Fig. 3.7 si divide in più parti. 
— un programma principale che non riempie alcuna delle funzioni figuranti nel 


flow-chart, 


— tre sottoprogrammi che effettuano rispettivamente le seguenti funzioni: 


— lettura dei dati, 


— calcolo della lunghezza dei lati del perimetro, 


— stampa dei dati 


e dei risultati. 
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Questo procedimento di esecuzione non è giustificato data la semplicità del 
problema, ma è un valido esempio di come procedere per la soluzione di problemi 
più complessi. 


100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
390 
400 
410 
420 
430 
440 
490 
500 
510 
520 
530 
540 
550 
560 
570 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 


P 


L 
FP 


L 
P 


REM CALCOLO DEL PERIMETRO DI UN FOLIGONO 
REM 

DIM XC(100),YC(100),LC(100) 

FRINT "FERIMETRO DI UN FOLIGONO" 


PRINT 

GOSUE 400 

GOSUE 500 

GOSUE 600 

DATA 5 

DATA 1,3,4,6,8,6,11,5,11,0 
STOF 

REM LETTURA DEI DATI 

READ N 


FORI = 1 TON 

READ XC(I),YC(I) 

NEXT I 

RETURN 

REM CALCOLO DELLA LUNGHEZZA DEL FERIMETRO 
= 0 

FORI =1 TON 


CI) = SOR COXCI) — XI + 1)) 4 2 + (YCI + 1) - YC(I)) 4 2) 
= P+ LAI) 
NEXT I 
(N) = SOR COXCON) — X(1)) A 2 + (YCON) — YC(1)) 4 2) 
= P + LCN) 
RETURN 


REM STAMFA DEI RISULTATI 

FRINT "NUMERO"} TABEC 10)}"X"} TABEC 22); "Y"} TAEC 30); "LUNGHEZZA" 
FRINT 

FORI = 1 TON 


FERIMETRO DI UN FOLIGONO 


NUMERO 


SAIBIQNE- 


EREAK IN 2 
1 


PRINT I} TAEC 10)}XCI)} TAEC 22)5YC(I)} TAEC 30);LCI) 
NEXT I 
PRINT 
PRINT 
PRINT TABC 22)}"FERIMETRO = "3 TAEC 30);F 
RETURN 
END 
Figura 3.7 
Xx Y LUNGHEZZA 
1 3 4,24264069 
4 6 4 
8 6 3.16227766 
11 5 5 
11 0 10,4403065 
FERIMETRO = 37,8452249 
00 
Figura 3.8 


3.4 TRACCIAMENTO DI UNA CURVA 


L’utilizzo di una stampante o di una macchina da scrivere come strumento per 
tracciare curve può, a volte, essere utile; in particolare quando non si ha la 
possibilità di avere un plotter, il problema è quindi trovare e mettere a punto un 
metodo che permetta di tracciare delle curve in condizioni accettabili. 
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Problema: Costruire un programma di tracciamento di curve secondo le seguenti 
istruzioni: 


1. Determinare il principio più semplice di lavoro per tracciare una curva y = f(x) 
per x che varia fra due valori dati A e B. 
— Come operare per ridurre l’effetto degli errori di arrotondamento ? 


2. Costruire il flow-chart di un sottoprogramma che effettua il tracciamento; e 
scrivere poi il sottoprogramma stesso: 


3. Scrivere in seguito il programma principale che chiama il sottoprogramma e 
provare a tracciare differenti funzioni come: 


=1/2 . 
— e ''’* cos 2x per x che varia da 0 a 10 


—_ e per x che varia da -2a + 2. 


Soluzione: Tenendo conto del fatto che con una macchina da scrivere è impossibile 
(salvo casi particolari) effettuare il ritorno della carta, e che è comodo far variare x 
incrementandolo, il metodo più semplice è quindi quello di scegliere l’asse y 
orizzontale orientato verso destra e l’asse x verticale orientato verso il basso. 


Figura 3.9 


Si pongono ora i seguenti problemi: 


— posizione degli assi, 
— come determinare per il valore dato di y, il numero di spazi da effettuare 
prima di stampare un punto? 


Queste due domande pongono il problema dell’arrotondamento. In effetti, una 
macchina da scrivere classica può effettuare solo un numero intero di spazi. Per un 
avanzamento teorico di y = 8,60 spazi, l'avanzamento reale sarebbe di 9 spazi. Per 
contro se y ha un valore di 8,40 gli spazi effettuati sarebbero, preferibilmente, 8. 


Poichè la funzione TAB normalmente opera per spazi interi, sarà necessario 


calcolare: 
Z = INT(y + 0,5) per ottenere l'arrotondamento 
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Prima di disegnare il flow-chart, bisogna risolvere il problema della posizione 
degli assi e del coefficiente di scala: come passare da y teorico a y reale sul terminale? 


Se la scala di variazione di y è: D = Y max - Y min, 


con L che rappresenta il numero massimo di caratteri per linea, si avrà un coeffi- 
ciente di scala K dato da: K = intero ((L — 1)/D). 
Il programma principale potrà quindi comunicare al sottoprogramma: 


=== Y max, Y min, Jo: 
— i valori di A e B, limiti di variazione di x 
— l’incremento di H per x 


Si suppone che Y min sia negativo così che l’asse x possa essere stampato nella 
colonna dalla quale ha origine. I punti della curva saranno rappresentati da punti 
ortografici. 


100 REM FROGRAMMA FER IL TRACCIAMENTO DELLE CURVE 

110 REM LA FUNZIONE FNA RAFFRESENTA LA CURVA DA TRACCIARE 
1200 DEF FN ACOX) = EXPO CO = XxX x 0,5) 

130 A= - 3 


170 GOSUE 500 


180 STOF 

500 D= Y2 - Yi 

510 K INT (CL — 1) / DI) 

520 Z INT (K * ABS C(Y1) + 0,5) 


FN ACX) — Yi 

INT CK x Z1 + 0,5) 

560 IF Z1 = Z THEN FRINT TAEC Z);"." 

570 IF ZI £ Z THEN FRINT O TAEC Z1);"."} TAEC Z);" 
580. IF Z1 > Z THEN FRINT TAEC Z);"I"} TAEC Z1);" 


530 FOR X = A TO E STEF H 


600 RETURN 


Figura 3.10 
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Figura 3.11 


3.5 CONCLUSIONI 


Questi esercizi danno l’impressione che la programmazione di formule matema- 
tiche non presenti difficoltà. In realtà è solo sufficiente scrivere le istruzioni 
opportune, il flow-chart rimane semplice (flow-chart lineare). L’esercizio di trac- 
ciamento di curve suppone già un flow-chart più complesso e una riflessione per la 
parte ‘‘presentazione“‘. 

Le suddette impressioni saranno confermate soprattutto quando incontreremo 
programmi che esigono manipolazioni complesse o che comportano numerosi test. 

Questo per esempio è il caso dell’esercizio ‘regine su una scacchiera‘ descritto 
nel Capitolo 10. 
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CAPITOLO 4 


ESERCIZI ORIENTATI ALLA GESTIONE 


4.0 INTRODUZIONE 


In questo capitolo ci si occuperà di esercizi semplici che presentano un interesse 
sul piano pratico e sul piano didattico. In effetti per le applicazioni gestionali, si ha 
frequentemente bisogno di effettuare delle ricerche, fusioni di tabelle o generazione 
di archivi. 

Gli esercizi proposti possono essere ripresi e incorporati in programmi di mag- 
giore complessità al fine di realizzare delle applicazioni complete. Per esempio la 
sequenza di ricerca data al paragrafo 4.1 può essere utilizzata per generalizzare il 


programma di funzione in 4.2 e in un certo modo può servire a generalizzare il 
programma della rubrica telefonica dato in 4.3. 


4.1 RICERCA SECONDO IL METODO DI SHELL 

Esistono numerosi metodi di ricerca di dati nella memoria del calcolatore. Il 
metodo più semplice è certamente quello della ““bolla‘ (in Inglese “bubble sort‘) 
che si trova in numerosi manuali. Qui viene proposto il metodo di Shell che, senza 
essere troppo complesso, fornisce le migliori prestazioni per quanto riguarda la 
velocità di esecuzione riducendo il numero dei confronti necessari. 


Si considera una tabella di N numeri da ricercare. Il metodo di Shell consiste 
nell’effettuare quanto segue: 


1. Determinate K in modo che: 


PE E e 


< seg è k 
- inizializzare D con il valore 2° - 1 


2. Effettuare una prima tappa di ricerca facendo variare l’indice I da 1 a N—D 
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2.1 Provare se A(I) < A(I+D) 
se si passare alla tappa successiva (in 3), 
se no scambiare A(I) e A(I+D) 
Porre k=I e passare alla 2.1 


2.2 Provare se A(K—D) < A(K) 
se si passare alla tappa successiva (in 3) 
se no scambiare, poi porre K = K—D 
e passare di nuovo in 2. 


3. Continuare i test incrementando I. Quando I ha assunto il valore N e quindi D > 
0, ritornare in 2 avendo calcolato a priori D= INT 


D_1I 
a GR 


Quando D ha assunto il valore 0, la ricerca è terminata. 
Problema: 


— disegnare il flow-chart di una subroutine di ricerca, 
— scrivere un programma che legga una tabella non ordinata e che chiami una 
subroutine di ricerca. 


Soluzione: Il programma può essere scritto facilmente utilizzando la precedente 
descrizione della tecnica di Shell. 


Scambiare Y e Z significa attribuire a Y il valore di Z e a Z il valore di Y. Si 
scriverà dunque: 


500 Y= Z 
SIOZ=Y 


In questo caso il valore di Y, essendo stato modificato dalla prima istruzione, non 
darà i risultati sperati nella seconda. Bisognerà inserire prima il contenuto di Y in 
una variabile ausiliaria X, si scriverà quindi: 


490 X = Y 
500 Y= Z 
SIOZ=X 


Questo tipo di scambio si trova nel programma (linee 590-610) e figura anche nel 
secondo programma della rubrica telefonica (4.5.2). 
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100 DIM ACIZ) 

110 N = 12 

120 PRINT "TABELLA TNIZIALE" 
180 PRINT 

140 FORT = 1 TON 

150 READ feI) 

160 FRINT ACI)" "; 

170 NEXT I 

180 5O£ 

190 PRINT 

195 PRINT 

PRINT "TABELLA ORDINATA" 
PRINT 

FORT = 1TON 

PRINT ACIO" "3 


NEXT I 
TOP 
= 1 
2% D 
IF D NOTHEN 510 
TNT CD 1) 2) 
TF_D e 0 THEN 700 
FORT = 1 TON D 
T 
dt D 
TEO fc) ASL) THEN 440 


RETIIEN 
DATA 3,-1,4,10,0,9,5,-10,-5,20,22,7 


END 


Figura 4.2 


Esempio di esecuzione: 


TABELLA TISIZIONLI 
2 1 4 10 93 9 5 10 5 PE 
TSBELLA ORDINATA 


10 Ò 1 3 4 S 7 8 9 10 22 25 


Figura 4.3 
4.2 FUSIONE DI DUE TABELLE 


Si richiede la fusione di due tabelle A e B contenenti numeri in ordine crescente, 
in una terza tabella C anch’essa ordinata secondo lo stesso criterio. 


Esempio: A 3 4 6 18 


Per ottenere: 
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Metodo da utilizzare: Utilizzare tre indici I, J, K per ciascuna tabella, questi indici 
saranno inizialmente ‘‘settati‘ a 1. 


Se ai < b; inserire a;in Cx, 
incrementare i e k. 


Se a:> bj inserire bj in Cx, 
incrementare j e k, 


quando una delle tabelle A e B è stata completamente trasferita in C, bisogna 
trasferire il resto dell’altra tabella in C. 


Problema: Disegnare il flow-chart, 
scrivere una subroutine in BASIC. 


Domande: a) Come fare se le tabelle A e B non sono ordinate? 
b) Come si può adattare il programma alla fusione di due files sequen- 
ziali ordinati? 


INIZIO 


RICOPIARE RICOPIARE 
A(1) INC B(J) IN C 
1=1+1 = JH1 


AFEII CÈIL 
COMPLEMENTO 
TRASFERITO IN C, TRASFERITO IN © 


TRASFERIRE TRASFERIRE 
IL RESTO 8INC IL RESTO AINC 


Figura 4.4 
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Soluzione: Il metodo proposto porta al flow-chart 4.4 che necessita di qualche 
spiegazione prima di poter cominciare la scrittura del programma. Il passaggio al 
flow-chart di Fig.4.5 è molto semplice a condizione che non si dimentichi di 
utilizzare tre indici. 


I indice per A, 
J indice per B, 
K indice per C. 


Nel programmare per evitare un GO TO, è preferibile inizialmente porre K a 0 e 
l’istruzione K = K + 1 situata nel mezzo del flow-chart all’inizio dell’algoritmo. 


Così anche i due piccoli loop possono essere scritti per mezzo delle istruzioni 
FOR ... NEXT a condizione di utilizzare una variabile ausiliaria. 


INIZIO 


C(K) = A(1) 
I=1+1 


NO 
RETURN RETURN 


Figuta 4.5 
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100 DIM AC100),E(100),C(100) 

110 REM LETTURA TABELLA A 

120 READ M 

130 PRINT "TABELLA A?" 

140 FORI = 1 TOM 

150 READ ACI)? PRINT ACI);: NEXT I 
160 PRINT 

170 FRINT 

180 REM LETTURA TABELLA E 

190 FRINT "TABELLA B:? " 

200 READ N 

210 FORI = 1 TON 

220 READ ECI): PRINT EC(I)}? NEXT I 
230 PRINT 

240 PRINT 

250 GOSUE 300 

260 FRINT "TABELLA Ci: " 

270 FORI= 1 TOM4+ N 

280 FRINT C(I);? NEXT I 

290. STOF 

295 REM SUEROUTINE FER LA FUSIONE DI Z E E 
300 I = 134 = 1îK = 1 

310 IF ACI) è» = EJ) THEN 350 

320 CK) = ACID:I = I+1 

330 IF I > M THEN 390 

340 K = K + 1? GOTO 310 

350 CK) = B(J)i4J=J+1 

360 IF Jz = N THEN 340 

365 REM COFIA IL RESTO DI A INC 
370 K = K + 1?CCK) = ACI)O:I = I+1 
380 IF I< = M THEN 370ELSE RETURN 
385 REM COFIA IL RESTO DI E INC 
390 K = K + 1?CCK) = B(J)iJ = J+ 1 


400. IF J < = N THEN 390ELSE RETURN 
410 DATA 5 
420 DATA 4,7,9,12,45 
430 DATA 4 
440 DATA -1,5,6,60 
450 END 
1 
Figura 4.6 


Esempio di esecuzione: 


TABELLA A: 

479 12 45 

TABELLA B: 

-1 5 6 60 

TABELLA Ci 

-14567 9 12 45 60 
EREAK IN 290 


Figura 4.7 


ESTENSIONI DI QUESTO PROGRAMMA 


Prima estensione: Come adattarlo a due tabelle A e B non ordinate. La prima idea 
che consiste nel fare una tabella unica non ordinata e di effettuare in seguito un 
ordinamento. Questo metodo conduce a un ordinamento più lungo di quello del 
secondo metodo che consiste nell’ordinare prima A e B ed effettuare poi la fusione 


(1). 
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INIZIO 
APERTURA 
FILE 
LEGGERE 
ARTICOLO 
FILE A 


LEGGERE 
ARTICOLO 
FILE B 


SCRIVERE ARTI- 
COLO A INC 


RICOPIARE 
ESTO DI BIN C 


RICOPIARE 
TESTO DI A IN Cl 


CHIUSURA 
FILE 


Figura 4.8 
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Per realizzare questi ordinamenti preliminari, si potrà utilizzare la sequenza 
precedente (linee dalla 500 alla 700 che bisognerà ripetere una volta per la tabella A 
e una volta per la tabella B) poichè la maggior parte degli interpreti e dei compilato- 
ri BASIC non accettano i sottoprogrammi con trasmissione di parametri (2). 


Seconda estensione: Fusione di due files sequenziali. 

Il flow-chart di Fig.4.4 costituisce un eccellente punto di partenza. Bisognerà 
includere le istruzioni di lettura e scrittura, pensare che raramente si conosce a 
priori il numero “degli articoli‘ dei files, dunque si prevederanno test periodici per 
la conclusione dei files. 


Questa estensione conduce al flow-chart di Fig.4.8. La programmazione dipende 
dal sistema utilizzato. In effetti le istruzioni per i files non sono normalizzate in 
BASIC. 


4.3 GIORNO DELLA SETTIMANA 


Partendo dalla data: GIORNO, MESE, ANNO, indicare di quale giorno della 
settimana si tratta. Questo problema è stato risolto con numerosi metodi, noi ve ne 
proponiamo uno: 

— Calcolare il valore del termine di correzione N, normalmente N = 0 ma se il 

mese è Gennaio o Febbraio N avrà per valore: 
2 se l’anno non è bisestile 
1 se l’anno è bisestile 
— Calcolare in seguito il ‘‘codice giorno“ con: 


C: = intero (365,25 x A2) + intero (30,56 x M)+J + N, 


AI rappresenta le due prime cifre dell’anno e, 
A2 le ultime due cifre, 
per esempio 1979 dà A1=19 e A2=79 


— Calcolare poi il numero S del giorno della settimana con 


S=C+3— 7xintero (==) 
7 

corrisponde a Lunedì, 

S = 2 corrisponde a Martedì, 

etc. 

S = 7 corrisponde a Domenica. 


(1) Per maggiori dettagli consultare pubblicazioni specializzate in algoritmi di ordinamento. 
(2) Questo costituisce un notevole vantaggio del FORTRAN sul BASIC. 
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Un anno è bisestile se: 
A2# 0 e A2 divisibile per 4, oppure A2 = 0 e AI divisibile per 4. 


Esempio: 1900 non è bisestile poichè 19 non è divisibile per 4. 
1984 è bisestile poichè 84 è divisibile per 4. 


4.3.1 Problema: Scrivere un programma che a partire dalla data G, M, A stampi il 
giorno della settimana corrispondente. 


4.3.2 Soluzione: Per poter utilizzare il programma in modo semplice, si costruirà un 
programma che domandi una nuova data fino a che il giorno introdotto sia 
negativo o nullo (0). Procedendo dal metodo proposto, la costruzione del flow- 
chart non presenterà difficoltà. 


E’ necessario conoscere: 
— come calcolare Al e A2, 
— se l’anno è bisestile. 


Si noti che Al è uguale al quoziente della ‘“divisione intera di A : 100‘ quindi: 
A1= INT (A/100) 


A2 rappresenta il resto della divisione intera precedente, quindi: 
A2=A— 100x Al 


Per sapere se A2 è divisibile per 4, è sufficiente calcolare il resto R per: 
R=A2—-4xINT(A2/4) 


Questo tipo di istruzione appare sovente negli esercizi del Capitolo 2. Si può ora 
scrivere una parte importante del programma fino al calcolo di S incluso. Bisogna 
studiare come si realizzerà la visualizzazione dei dati, vengono presentati differenti 
casi: 

Primo caso: Il sistema accetta delle tabelle di stringhe di caratteri. In questo caso è 
sufficiente utilizzare le seguenti istruzioni: 


DIM G$ (7) 


G$ (1) = “LUNEDI” 
G$ = “DOMENICA“ 


e stampare il giorno con: 


PRINT G$ (S) 
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INIZIO 
LETTURA 
DEI DATI 


CALCOLO DI A1 E A2 


CALCOLO DI C 
CALCOLO DI S 


VISUALIZZAZIONE 
GIORNO DELLA 
SETTIMANA 


Figura 4.9 
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Secondo caso: Il sistema non accetta delle tabelle di stringhe. Si potrà allora 
impostare una stringa di caratteri abbastanza lunga per contenere tutti i giorni della 
settimana. 


Il giorno più lungo è MERCOLEDI” con 9 caratteri. Una stringa di 9 x 7 = 63 
caratteri potrà dunque essere inizializzata convenientemente è sarà sufficiente per 
stampare la sottostringa per mezzo di un'istruzione del tipo: 


PRINT SUBSTR (G$,8 * S — 8,9 ) 


oppure 


PRINT MID $ (G$,8* S— 8,9) 


90 DIM G$(56) 

95 G$ = "LUNEDI MARTEDI MERCOLEDIGIOVEDI VENERDISABATO DOMENICA" 
100 REM CALCOLO DEL GIORNO DELLA SETTIMANA 

110 REM S RAFFRESENTA IL NUMERO DEL GIORNO (1 FER LUNEDI E 7 FER 

115 REM DOMENICA) 

120 INFUT "IMFOSTARE LA DATA "}G,M,A 


130 Ai = INT CA / 100) 

140 AZ = A - 100 x Al 

150 N= 0 

160 IF Mè= 2 THEN 300 

165 N = 2 

170. IF A2 = 0 THEN 220 

180 R = A2 - 4 x INT C(A2 / 4) 

190 IFR<& 0 THEN 300 

200N= 1 

210 GOTO 300 

220 R = Al - 4 * INT (Al / 4) 

230 IF R= 0 THENN = 1 

300 C = INT (365,25 * A2) + INT (30,56 Xx M) +N+G 
310 S= 3 +C - 7 Xx INT (CC + 2) / 7) 


320 FRINT MID$ (G$,9 x M - 8,9) 
330 PRINT 
340. GOTO 120 


Figura 4.10 


Nota: Nel corso dell’esercizio 4.4 vedremo l’utilizzazione di una funzione utente per 
ridurre il numero di istruzioni del programma. 


Il MICROSOFT BASIC non utilizza la funzione SUBSTR ma la funzione MID$. 
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IMFOSTARE LA DATA ?12,4,1979 
GIOVEDI 

IMFOSTARE LA DATA ?13,4,1979 
VENERDI 

IMPOSTARE LA DATA ?15,8,1979 
MERCOLEDI 

IMFOSTARE LA DATA ? 

EREAK IN 200 


Figura 4.11 


4.3.3 Analisi di un programma 


Il programma di Fig.4.12 permette ugualmente di ottenere il giorno della setti- 
mana utilizzando un altro metodo. 


100 REM FROGRAMMA FER CALCOLARE IL GIORNO DELLA SETTIMANA 
110 DIM G$C(7) 


120 G$(1) = "LUNEDI" 
130 G$(2) = "MARTEDI" 
140 G$(3) = "MERCOLEDI" 
150 G$(4) = "GIOVEDI" 
160 G$(5) = "VENERDI" 
170 G$(6) = "SABATO" 
180 G$(7) = "DOMENICA" 


190 FRINT "IMFOSTARE LA DATA "; 
200 INFUT G,M,A 

210 GOSUE 500 

220 FRINT G$(Z) 

230 GOTO 190 


500 IF A « = 1752 THEN 620 
SION = INT (0,6 + 1 / M) 
520L=A- N 

530 F = M+ 12X N 

540 D=L / 100 

550 Al = INT (CD) 

560 Z1 = INT (D / 4) 

570 Z3 = INT (5xL / 4) 

580 Z4 = INT (13 x (FP + 1) / 5) 


590 Z = Z4 + z3- A1 + Z41+G+5 

600 Z= Z- (7 x INT (Z / 7)) + 1 

610 RETURN 

620 FRINT "L'ANNO DEV'ESSERE FOSTERIORE A 1752 " 
630 STOF 

640 END 


Figura 4.12 


IMFOSTARE LA DATA 5,2,1979 
LUNEDI” 

IMFOSTARE LA DATA 13,5,1979 
DOMENICA 

IMFOSTARE LA DATA 14,7,1979 
SABATO 

IMFOSTARE LA DATA 


Figura 4.13 
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Domande: 


l. Cosa significa l’istruzione della linea 510 ? 


2. Come viene interpretata l’istruzione di linea 530 ? 
3. E’ possibile, senza modificare il metodo e senza aumentare il numero delle 


operazioni, ridurre il numero delle istruzioni del programma ? 


Soluzione: 
1. Istruzione di linea 510 
— se M vale 1 0 2, N assume il valore intero (0.6 + 1) o intero (0.6 + 1/2), che in 


ogni caso dà 1, 
— se M vale 3, 4, etc., N assume il valore 0. 


Questo è dunque da rapportare al programma precedente dove si teneva conto 
delle conseguenze provocate dall’anno bisestile. 


2. Istruzione della linea 530 

— P corrisponde al numero del mese se questo è Marzo, Aprile, etc. fino a 
Dicembre. Per Gennaio e Febbraio P assumerà rispetivamente il valore 13 e 
14. 
D-A1 nell’istruzione della linea 590 ha un valore V che è V=0 se si trova in un 
anno secolare. 


0<V< 1 se si trova in un altro anno. 


3. Poichè le variabili A1, Z1, Z3 e Z4 sono utilizzate solo nel calcolo di Z nella 
linea 590, si possono sopprimere le istruzioni 560, 570 e 580 e scrivere in 590: 


590 Z= INT (13-(P+1)/5)+INT(S-L/4)—INT(0)+INT(D/4)+J+5 


4.4 DIFFERENZA TRA DUE DATE 


Per conoscere la differenza tra due date si può calcolare “il codice giorno‘ di ogni 
data e fare la differenza. Si ottiene così il numero dei giorni che le separano. Questa 
informazione è indispensabile quando si debba calcolare l'ammontare di interessi 
per un periodo qualunque. 


Problema: Riferendosi al programma precedente costruire un nuovo programma 
che calcoli la differenza tra due date. 
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LEGGERE 

1° DATO 
SIP CALCOLO 
CODICE PERC 


LEGGERE 
2° DATO 
SIP CALCOLO 
CODICE C 


PER MEMORIZZARE IL RISULTATO IN VISTA DI UNA 
_ © VISUALIZZAZIONE ULTERIORE (QUI NON NECESSARIA) 


Ce C3 RAPPRESENTA LO SCARTO CERCATO 


VISUALIZZAZIONE 
DEL RISULTATO 


Figura 4.14 


Soluzione: E’ possibile verificare l’importanza del sottoprogramma, si dovrà in 
effetti calcolare per due volte successive il codice giorno e sarà preferibile costruire 
un sottoprogramma per questo calcolo. Si può dunque realizzare il flow-chart di 
Fig.4.14. 

Riprendendo la parte del programma di Fig.4.10 che calcola il ‘codice giorno“, 
si costruisce rapidamente il programma 4.15. 


4.5 RUBRICA TELEFONICA 


Il BASIC presenta alcuni vantaggi in confronto ad altri linguaggi, come l’‘anti- 
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100 REM CALCOLO DELLA DIFFERENZA TRA DUE DATE 
105 REM AUTORE 4J.F. LAMOITIER 

110 DEF FN RCA,D) = A - Dx* INT CA / D) 

120 INFUT "IMPOSTARE LA PRIMA DATA "3}G,M,A 

130 GOSUB 500 

140 C1 = C 

150 INPUT "IMPOSTARE LA SECONDA DATA "3}G,M,A 
160 GOSUB 500 


170 PRINT ; 

180 PRINT "TRA LE DUE DATE VI SONO "}C - C1}" GIORNI" 
200 STOP 

500 A2 = FN RC(A,100) 

5I10N= 0 

520 IF M > 2 THEN GOTO 570 

525 IFM«x = 2 THENN = 2 

530 IF A2 = 0 THEN GOTO 550 

535 IF A2 < > 0 THENR = FN RC(A2,4) 
540. IF R< > 0 THEN GOTO 570 

545 IF R = 0 THEN GOTO S60 

550 IF FN RCA1,4) <« > 0 THEN 570 
560 N = 1 


570 C = INT (365,25 x A2) + INT (30.56 Xx M) +N + G 
580 RETURN 
590 END 


Figura 4.15 


Esempio d’esecuzione: 


IMFOSTARE LA PRIMA DATA 23,2,1979 
IMFOSTARE LA SECONDA DATA 30,6,1980 
TRA LE DUE DATE VI SONO 493 GIORNI 


Figura 4.16 


co‘* FORTRAN, tra gli altri quello di trattare molto semplicemente le stringhe di 
caratteri. I due esercizi seguenti mostrano come si possono facilmente manipolare 
le stringhe di caratteri in relazione ad applicazioni pratiche. 


4.5.1 Esercizio 1: Scrivere un programma che legga delle linee DATA contenenti 
cognome, nome, codice, numero telefonico e che stampi in modo preciso dette 
linee. Si daranno a queste stringhe di caratteri rispettivamente i nomi N$, P$, B$ e 
T. Per questo esercizio si supporrà che la lista dei dati sia in ordine alfabetico. 


x x 


Soluzione: Il flow-chart di principio è semplice poichè è sufficiente leggere e 

stampare senza effettuare elaborazioni. La sola difficoltà è conoscere quando la 

lista dei dati è terminata. Due metodi possono essere convenientemente utilizzati: 
— utilizzare un’ultima registrazione particolare che ‘‘precede‘* la fine della lista, 
— utilizzare l’istruzione IF...END propria di alcuni BASIC. 


Qui verrà utilizzato il primo metodo che è alla portata di tutti i sistemi. 


Alla fine della lista, si aggiunge un nome particolare che è “ZZZ e quindi 
facilmente riscontrabile. In questo modo è possibile conoscere il termine della lista. 
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INIZIO 
VISUALIZZAZIONE 
INIZIALE 


LEGGERE 
ARTICOLO 


STAMPARE IL 
NUMERO DELLE 
PERSONE | 


VISUALIZZAZIONE, 
N$, P$, B$, T 


Figura 4.17 


Tutto questo conduce al flow-chart di Fig.4.17. 

In questo flow-chart è stata aggiunta una variabile I che conta‘ il numero reale 
di persone contenute nella lista. Così, per liste importanti, contando il numero delle 
linee di uscita (il numero delle persone) è possibile prevedere dei salti di pagina al 
fine di ottenere una buona presentazione. 


Questo programma, che non effettua alcuna ricerca, permette tuttavia di ottene- 
re un elenco ordinato se i nomi sono stati introdotti correttamente, per questo 
motivo le linee sono numerate di 10 in 10 così che altre linee possano essere inserite 
dove necessario. 


Mediante qualche piccola modifica (a livello di istruzione di lettura) si può 
lavorare come per un file sequenziale, che non limita la dimensione dell’elenco. 
Bisognerà nel frattempo assicurarsi del continuo aggiornamento del file. 


Nota: Alcuni BASIC permettono di ‘testare‘* la fine del file con l’istruzione 


IF....END; ciò evita di dover impostare una registrazione fittizia (come ZZZ). Il 
flow-chart è presente in Fig.4.19. 
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FUOCRAIMA PER LÀ RUBRICA TELEFONICA 


Val i 0 "EUERTRÀ TELCPORICA” 


f 
Ù 
I 
r 

140 Pr 

ì Ù 
I 
I NODGNORE": SPE 6) pMOME"S SECC 7) pNUMERO"S SPEC 99) "INTERNO" 
Ù 


200 FFAD 04,4, ET 

PAOOCTE DA = SZ THEN 250 

PR0O PRINT CEI SPCC 6)5N$5 SFICO 10){E$} SFCC 14);T 
2230) Ta 1 

24 N SOTTO 20N 

260 PETIT 

240 PRINT "NUMERO DELLE PERSONE = "SI 
PINO DATA FOSSI, ANDEEA,"9",210 
CONOODATA  FOSST,GTANNI,"3",2340 

PON DATA VFRDI, LUIGI, "5,400 

2000 DATA VINCI, CLAUDIO, "4",360 

808 DATA Z27,Z,3,4 


RUERICA TELEFONICA 


COGNOME NOME NUMERO INTERNO 
ANDREA 3 310 
GIANNI 3 340 
LUIGI 5 400 
CLAUDIO 4 360 


ONE = 4 
Figura 4.18 


NUMERO DELLE F 


VISUALIZZAZIONE 
DEL NUMERO 

DELLE PERSONE 
DELLA LISTA 


lire N$, P$, B$ et T 
VISUALIZZAZIONE 


Figura 4.19 
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4.5.2 Esercizio 2: Si tratta ora di costruire un programma più complesso che 
propone un ‘‘menu“ all’utilizzatore: 

ricerca di un cognome, 

ricerca di un cognome e nome, 

ricerca del solo nome, 

- ricerca del numero telefonico, 

- estrazione delle persone che hanno un interno, 


e che effettui in seguito il lavoro richiesto dall’utilizzatore. Per la ricerca si può 
riprendere la sequenza di ricerca già esaminata, a condizione che venga adattata al 
problema da risolvere. 


Problema: Costruire un programma che legga i dati che si suppone siano immagaz- 
zinati nelle istruzioni DATA del programma, e che proponga poi il seguente menu: 


In funzione della risposta data dall’utilizzatore, il programma effettua il lavoro 
selezionato poi propone di nuovo il menu. Il programma termina se l’utilizzatore 
introduce un numero «0 o = 6. 


Nota: Il BASIC non accetta, in generale, delle subroutine con trasmissione di 
parametri, questo è invece possibile con numerosi linguaggi come ALGOL, FOR- 
TRAN, etc. Si incontreranno difficoltà nèl costruire la subroutine di ricerca che 
possa funzionare correttamente in tutti i casi. Si potrà quindi considerare di 
trasferire “la chiave di ricerca‘ in una tabella ausiliaria. Bisognerà far attenzione a 
non dimenticare le permutazioni. 


Soluzione: Questo problema non presenta difficoltà particolari a condizione di 


lavorare con metodo. Si costruirà prima un flow-chart generale non dettagliato 
(Fig.4.20). 


Per utilizzare il sottoprogramma di ricerca già esaminato, sarà necessario ogni 
volta preparare i dati che possono essere usati dallo stesso. E’ stata scelta la 
seguente convenzione: 


I dati sono €C$ (I) cognome 
N$ (I) nome 
D$ (I) codice 
T$ (I) numero telefonico 


Una tabella ausiliaria B$ (I) riceve prima gli elementi da ricercare e in seguito 
viene effettuata la ricerca. Prima di effettuare la ricerca sul cognome si esegue la 
sequenza: 


FORI=1TON 
B$ (I) = C$ (1) 
NEXT I 
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LETTURA DELLE MATRICI 
COGNOME, NOME, CODICI 
E NUMERI DI TELEFONO 


PROPORRE IL MENU 
LEGGERE LA 
SCELTA | 
1 = INTERO (1) 


VIRGOLETTE 


RICERCA 


RICERCA SU RICERCA ESTRA- 
; RICERCA 
SUL 
Ki NOME COGNOM SUL ZIONE 
ME E COGNOME N° TEL. CODICE 


Figura 4.20 


Questo presuppone dunque che si debba effettuare una ricerca alfabetica, il che 
non è una difficoltà in BASIC con i codici ASCII o EBCDIC. 

Al contrario per la ricerca del cognome e del nome, sarà necessario porre in B$ la 
concatenazione C$ + N$, ma in questo modo si corrono dei rischi; si considerino i 
seguenti nomi: 


ROSSI UGO 
ROSSINI ALDO 


83 


per concatenazione semplice si avrà: 


ROSSIUGO 
ROSSINIALDO 


e il confronto darà ROSSINIALDO < ROSSIUGO. 
Per evitare tale rischio si inserirà tra cognome e nome un carattere “spazio” che 
in ASCII possiede un valore inferiore alla lettera A. Si avrà quindi: 


ROSSI UGO 
ROSSINI ALDO 


così il confronto fornirà i risultati voluti. Per inserire il carattere ‘‘spazio‘ è 
sufficiente scrivere: 


470 B$ (I) = C$ (I) + “A “ + N$ (1) 


spazio tra le virgolette 


100 REM FROGRAMMA DI RUBRICA TELEFONICA 

110 REM 

120 DIM S$C100),C$(100),N$(100),T(100),E$(100) 
130 GOSUE 800 

140 FRINT "INDICAR 
1500 FRINT "1 = KR 


E IL TIFO DI RICERCA" 
2A DEL COGNOME" 


160 FRINT ERCA DEL COGNOME E NOME" 
170 FRINT ZA DEL NOME" 
180 FPEINT 2A DEL NUMERO TELEFONICO" 


190 PRINT 
200 INFUT 
210 REM FROVA DI VALIDITA” 


CAZIONE DEL NUMERO DI CODICE" 


220 IX = INT CI) 
2200 IF I = 0 THEN STOF 
240 IFIES= 6 THEN STOP 


250 REM SELEZIONE DEL LAVORO 

260 ON I GOTO 400,4560,520,580,640 
270 __STOF 

290 REM RICERCA DEL COGNOME DÌ 

in AR LI si domanda solo la ricerca 
410 ECT) = 0$C1) sul cognome 

420 NEXT I 

4230 COSUE 1000 

440 GOTO 1 

4500 REM RICERCA DEL COGNOME E DEL NOME : . 

430 FORT = 1 TO N si domanda solo la ricerca 


470 ESCL) e C$0I) +" " + N$CI) sul co nome i 
4800 NEXT IT 8 e il nome 


420 GOSUE 1000 

UNO COTTO 140 

10 REM RICERCA DEL NOME 
520 FORT = 1 TON 

Fid Cc 


si domanda la ricerca 
= N$CI) sul nome 


T) 


NUMERO TELEFONICO 


TOA ricerca sul numero 


590 STRE <TCI)) telefonico 
400 : 
410 : 1000 
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620 GOTO 140 

4300 REM FSTRAZIONE DEL NUMERO ICE 
6400 PRINT "IMPOSTARE TL CODICE SCELTO "3; 
6500 TNFUT A$ 

6400 PRINT 

4204) Ù) 

#20 PRINT "ELENCO DEI NOMI CON CODICE "3;A$ 
49N *REINT 

700 FORT := 1 TON 

710 TOSSE T) « : AD THEN 72 

PISO FPINT CECI), N60), 86 CDI, TCI) 
717 Ju do 1 

7220 NEXT I 

7309 PRINT 

ZAN PRINT "IN TOTALE "505" NOMI NELL ELENCO" 
250 COTO 140 

220 REMO SOTTOPROGRAMMA DI LETTURA 
snn FESDO N 

S1)0 FOR T=1T0N 

9200 READ D6&C01),N$(I),5$(T), TI) 
9200 TF C8(0) = "ZZZ" THEN 860 

240 NEXT T 

840 N = TIT L 

8200 RETURN 

990 REM SOTTOFROGRAMMA DI RICERCA SECONDO Il 
1000 Ds 1 

1010 p=2xD 

1020 TFNDE = N THEN 1010 

1020 D = INT (D- 13 # 2) 

1040 TF Dos D THEN 1400 

1050 FOR T= 1 TON - D 

1040 FOR Js= I TO 1 STEP - D 

1070 L= d+D 

1080 IF E$(.J) = = B$(L) THEN 135 
1090 X$ = C$(C64) 

1100 CE) s C$CL) 

1110 CEL) = X$ 

1120 X& = N$(J) 

1130 N$C(J) = N$(L) 

1140 N$CL) = X$ 

1150 X$ = S$(J) 

1160 S5$(4) = 56$(L) 

1170 S$(L) = X$ 

1180 X = T(4) 

1190 TH) # TL) 

1200 TL) = X 

1205 X$ = B$(J) 

1206 B$(J) = B$(L) 

1207 B$(L) = X$ 

1210 NEXT 4 

1350 NEXT T 

1360 GOTO 1030 

1400  GOSUE 1500 

1410 RETURN 

1490 REM SEQUENZA DI STAMPA 

1500 FRINT 

1510  FRINT "COGNOME", "NOME", "CODICE", "NUMERO" 
1520 FRINT 

1530 FOR K = 1 TON 

1540 FRINT C$(K),N$(K),5$C(K),TCK) 
1550 NEXT K 

1560 FRINT 

1570 FRINT "IN TOTALE "iN3" NOMI NELL’ELENCO" 
1580 RETURN 

2000 DATA 10 

2010 DATA ROSSI,FIERO,BE,100 

2020 DATA ROSSINI,GIANNI,BE,110 
2030 DATA VERDI, CLAUDIO,LAE:0,310 
2040 DATA ZANI,DARIO,FA,115 

2050 DATA ZINNI,FRANCO,COM,300 
2060 DATA CARLI,FAOLO,SEC,301 
2070 DATA ZZzZ,2Z,2Z,0 


estrazione del numero 


di interno 


METODO DI 


SHELL 


Figura 4.21 
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Esempio di esecuzione: 


INDICARE Il TIFO DI RICERCA 


1 MERCA DEL COGNOME 

2 CA DEL COGNOME E NOME 

3 N CA DEL NOME 

4 = RI 2A DEL NUMERO TELEFONICO 

5 = E ZIONE DEI. NUMERO DI CODICE 

P1 

COGNOME NOME CODICE NUMERO 

CARI.I FADLO SEC 301 

FOSSI FIERO E 100 

CONT GIANNI EE 110 

CLAUDTO LABO G10 
DARIO FA 115 
FRANCO COM 300 


DEL COGNOME E NOMI 

DEL NOME 

DEL NUMERO TELEFONICO 
ESTRAZIONE DEL NUMERO DI CODICE 


COGNOME NOME UODILCL SUmizEt 
SCI SUD! 
(HS 100 
GIANNI Fit 16 
CLAUDTO LALO du 
DARIN Fri LIS 
FRANCO Cis 100 


TNOTOTALE 6 NOMI NELL ‘ELENCO 
TNDICARE IL TTEO DO RICERCA 


1 DEI NOME 

7 A DEL COGNOME E NOMI 

k) (DA DEL NO 

4 CRA DER MUMERO TELLIONICO 

Moe TRAZIONE PELO MUMERO DT CONTRE 

5a 

COGNOME NOME MUMIRDO 
VERDT CLAUDIO LARG 2310 
ZANI DARTO FA Lit 
ZTNNT FRANCO cor ano 
ROSSINI GTANNI Lib 14U 
CARLT FADLO SEL 301 
ROSST FIEEO GE 100 


Figura 4-22 (continua) 


IN TOTALE 6 NOMI NELL ‘EL 


CO 


TNDICARE IL TIFO DI RICERCA 

1 =* RICERCA DEL COGNOME 

nas CA DEL COGNOME E NOME 

3 A DEL HOME 

4 | :DA DEL NUMERO TELEFONICO 

© = FETRAZIONE DEL NUMERO DT CODICE 

DA 

COGNOME NOME CODICI NUMERO 
ROSSI FIERO EE 100 
ROSSINI GIANNI BE 110 
ZANI DARIO FA 135 
ZINNT FRANCO DUM ano 
CARLI FADLO SEI 301 
VERDI CLAUDIO LAEO 310 
IN TOTALE 6 NOMI NELL’ELENCO 

INDICARE IL TIFO DI RICERCA 

1 = RICERCA DEL CUGNOME 

2 RCA DEL COGNOME E NOME 

3 RCA DEL NOME 

4 = RICERCA DEL NUMERO TELEFONICO 

5 = ESTRAZIONE DEL NUMERO DI CODICE 

25 

IMPOSTARE IL CODICE SCELTO ?EE 

ELENCO DEI NOMI CON CODICE BE 

ROSSI FIERO EE 100 
ROSSINI GIANNI EE 110 


IN TOTALE 2 NOMI NELL’ELENCO 
INDICARE IL TIFO DI RICERCA 


1= RCA DEL COGNOME 

2= CA DEL COGNOME E NOME 

3 = RICERCA DEL NOME 

4 = RICERCA DEL NUMERO TELEFONICO 
5 = ESTRAZIONE DEL NUMERO DI CODICE 
?0 


EREAK IN 230 


Figura 4.22 (fine) 


A partire dal momento in cui la ricerca è alfabetica, la ricerca del numero di 
telefono si può ottenere trasformando questi numeri in stringhe di caratteri da cui la 
linea 590 


590 B$(1) = STR$ (T(1)) 


La sequenza della ricerca assomiglia molto a quella dell’esercizio 4.1 ma le 
permutazioni devono essere effettuate su C$, N$, S$, T e anche su B$. 


Critiche al programma: Il grosso difetto di questo programma è la dimensione 
limitata dell’elenco dovuta al fatto che i dati sono inclusi nel programma e risiedo- 
no durante l’esecuzione nella memoria del calcolatore. Tuttavia con la stessa 
struttura generale è possibile costruire un nuovo programma che lavori con accesso 
diretto al file su disco, in questo caso è necessario porre attenzione a non aumentare 
inutilmente il numero di accessi al disco. 
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4.6 CONCLUSIONI 


Gli esercizi orientati alla gestione non presentano particolari difficoltà dato che 
non richiedono l’accesso a file. E’ tuttavia necessario lavorare con metodo. 


Gli esercizi che fanno appello ai file non sono più difficili a livello di flow-chart 


ma la loro programmazione dipende dall’interprete BASIC impiegato, che è rara- 
mente compatibile per quanto concerne le istruzioni di trattamento dei file. 
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CAPITOLO 5 


CALCOLI MATEMATICI 


5.0 INTRODUZIONE 


Il BASIC si presta bene all’esecuzione di semplici calcoli matematici. General- 
mente i flow-chart sono molto semplici da realizzare e la programmazione presenta 
poche difficoltà. 


Per contro, in certi casi, il propagarsi degli errori di arrotondamento provoca il 
raggiungimento di risultati poco precisi. Per premunirsi contro l’influenza di questi 
errori è possibile l’attuazione di diverse tecniche: 

— prima scegliere degli algoritmi poco sensibili a questo tipo di errore, ciò non è 
sempre facile, 

— in seguito programmare in modo che la perdita di precisione risulti più piccola 
possibile. 

Nel contesto di un libro di esercizi, non è possibile dilungarsi su problemi di 
questa importanza, il lettore interessato potrà consultare delle pubblicazioni di 
analisi numerica. Nel frattempo, per attirare l’attenzione del lettore sulle possibili 


conseguenze, lo invitiamo a studiare attentamente i risultati ottenuti dal calcolo di 
“pi greco‘ con il metodo dei poligoni regolari inscritti. 


5.1 DIVISIONE DI UN POLINOMIO PER (x - 5) 
Si consideri un polinomio p(x) di grado n di cui si conoscono i coefficienti: 
P(x) = ao x" + arx" + ax? +... + ani + an 
Bisogna trovare il polinomio Q(x) di grado n-1 così che: 
P(x)=(x—s5Qx)+R 
il resto R è una costante. Quindi: 
Q(x) = box" + p. +... + baz + dei 
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si otterà: 

bo= ao 

bi= a1+ sbo 

b-= a-+ sbi eR=an+0B", 
B-= ai+ sbi 


bni = an! + SDn-2 


Problema: Scrivere un programma che determini i coefficienti di Q(x) partendo dai 
coefficienti di P(x) e di (x—A1), essendo A1 ugualmente dato. 


Soluzione: La parte relativa all’algoritmo è particolarmente semplice poichè per I 
che varia tra 1 e N—I si può scrivere: 


bi =, 41 + Sbit «_______—__ Coefficiente di Q(x) 
Coefficiente di P(x) 


Si può comunque calcolare bn con questa formula e fare: 
R= bn 
Tutte le difficoltà consistono dunque nello scrivere bene le istruzioni di I/O. 


Il programma di Fig.5.1 costituisce un esempio di ciò che si può fare. E’ possibile 
evidentemente avere altre presentazioni dei risultati. 


20 REM DIVISIONE DI UN FOLINOMIO FER X-Al 

30. REM N = GRADO DEL FOLINIMIO 

40 REM IL VETTORE ‘A’ CONTIENE I COEFFICIENTI DI F(X) 
50 REM IL VETTORE “E” TIENE IL CALCOLO DEI COEFFICIENTI DI Q(X) 
60 REM “R” E” IL RESTO DELLA DIVISIONE TRA F(X) E X-Al 
100 DIM AC20),E(21) 

105 REM LEGGE UN INFUT 

110 READ N,Al 

115 PRINT "DIVISIONE DI F(X) FER (X-"3A13") 

116 PRINT 

120 FORI = 1 TON 

130 READ ACI) 

140 NEXT I 

145 REM CALCOLO DEI COEFFICIENTI DI Q(X) 

150 E(0) = AC0) 

160 FORI =1TON- 1 

170 BECI) = ACI) + Al x ECI —- 1) 

180 NEXT I 

190 R = ACN) + Al x ECN — 1) 

195 REM STAMFA DEI RISULTATI 

200 FRINT "COEFFICIENTI DI F(X) "; 

210 FORI = 1 TON 

220 PRINT ACI); 

230 NEXT I 

240 PRINT 

250 PRINT 
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260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 


Esempio di esecuzione: 


PRINT "COEFFICIENTI DI Q(X) "3; 
FORI = 0 TON - 1 
PRINT ECI); 

NEXT I 

PRINT 

PRINT 

PRINT "RESTO? "3R 
STOF 

DATA 6,1 

DATA 3,2,-1,5,6,4,1 
END 


Figura 5.1 


DIVISIONE DI F(X) FER (X-1) 
COEFFICIENTI DI F(X) 32-1564 
COEFFICIENTI DI Q(X) 0354915 
RESTO? 19 

BREAK IN 330 


Figura 5.2 


Si verifica facilmente che: 


3x+2x°-x'+5x° +6x°+4x+1= 


(x-1) (3x° +5x*+4x°+9x°+15x+19) +2 


Commenti: 


1. In parecchi sistemi l’istruzione BASE 0 è destinata a mostrare che le variabili 


indicate hanno l’indice che parte da 0. 


Secondo le norme standard ANSI l’istruzione BASE 0 dovrebbe essere scritta 


OPTION BASE 0. 


2. Con un BASIC dove che gli indici devono obbligatoriamente partire da 1, è 
sufficiente sottrarre 1 dall’indice nel programma di figura 5.1. 

3. La sequenza di stampa utilizzata può evidentemente essere modificata in 
funzione di una concreta applicazione. Si può notare tuttavia la necessità di inserire 
delle istruzioni PRINT per rendere più chiara la presentazione. 


5.2 CALCOLO DI UN INTEGRALE DEFINITO 


Tra i numerosi metodi che permettono di calcolare l’integrale definito di una 
funzione limitata su un intervallo ugualmente limitato, si possono citare: 


— il metodo di SIMPSON 
— il metodo di WEDDLE 


che sono particolarmete semplici e saranno oggetto di due esercizi. 


Metodo di Simpson: Con n pari, per calcolare: 


s=S fm dx 


Si opera nel modo seguente: 

— dividere l’intervallo a, b, in n intervalli con: 
_b-a 

n 


h 


— calcolare poi: 


S= STICA) +4 f(x) +2 f(x) +4 f(x) +2 f(x) + +4 fn) + fn] 


con Xo = 4 € Xn = d. 


Metodo di Weddle: Con n multiplo di 6, si calcola 4 come in precedenza, poi si valuta 
S con: 


3 h 
S= gg H(0) +5 f(a + h) + f(a +2h)+6f(a +3 h)+f(0 +4 h)+ 


+5 f(a +5 h)+ f(b)] 


sen=60osen= 12,18, etc. 


3} 
= UA) +5 fa +h)+f(a+2h)+6f(a +3 h)+ f(a +4 h)+ 


+Sf(a+5h)+2 f(a +6 h)+ + f(b)] 


Esercizio 1: Scrivere una subroutine che calcoli un integrale definito con il metodo 
di Simpson; scrivere un programma principale che utilizzi tale metodo per calcola- 
re: 


1 


s- | di le * dx 
VI 


1:95 


Fare un calcolo con n = 6, 12, 18, 24, 30 senza tener conto del fatto che questa 
funzione è “pari‘*. 


Esercizio 2: Mantenendo lo stesso enunciato precedente si sostituisce al metodo di 
Simpson quello di Weddle. 


Effettuare il calcolo con n = 6, 12, 18, 24. 
Confrontare i risultati ottenuti. 


Soluzione esercizio 1: 


Per costruire un solo loop di calcolo definiamo due variabili S1 e S2 costituite 
rispettivamente dai termini di S con coefficiente 4 e 2. 
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S1 = f(x) + f(x) + +f(x,-1) 
S2 = f(x) + f(x.) + +f(x,-3) 


sostituendo S si avrà 
H 
S= 3 S1+252+ f(a)+ f(b)) 


Il flow-chart relativo alla parte di calcolo si costruisce ora facilmente (Fig.5.3) a 
condizione di ricordare che S1 deve normalmente contenere un termine superiore a 
quello di S2 e che si deve effettuare un loop comune, sarà necessario alla fine 
aggiungere f (xn-1) che sarà stato dimenticato. 


Soluzione esercizio 2: 


Si farà riferimento al metodo dell’esercizio 1 con la differenza che in questo 
esercizio si faranno 6 somme parziali nel loop e, come precedentemente, una 
somma definita all’uscita del loop. Questo metodo allunga notevolmente il pro- 


gramma ma lo rende più semplice. Per questa ragione qui non viene raffigurato il 
flow-chart relativo. 


S1=S1+F(X) 
X=X+H 
S2=S2+ F(X) 
I=I+1 


S=(4S1+4F[X +H]+2 S2+ 
F[A]+F[B]) 5 


CORRISPONDEA f(x,_,) bd 


Figura 5.3 
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10 REM CALCOLO DI UN INTEGRALE CON IL METODO DI SIMFSON 

50 PRINT "CALCOLO DI UN INTEGRALE DEFINITO CON IL METODO DI SIMFSON" 
60 FRINT 

70 PRINT " N INTEGRALE" 

80 PRINT 

100 DEF FNECX)=C*EXF (-XXX%X0,5) 

110 C=1/80QR(2X3,14159267) 

112 A=-1,95 

113 E=1,95 

115 FOR N= é TO 30 STEF 6 

120 GOSUE 3500 

130 PRINT USING 200,N,5 

140 NEXT N 

150 STOF 

200 3484 IRE AGGGzzziLi 

3400 REM SOTTOFROGRAMMA DI CALCOLO DELL'INTEGRALE DEFINITO 
3410 REM METODO DI SIMFESON 

3420 REM ‘“H” RAFFRESENTA IL FASSO DI INTEGRAZIONE 

3500 XF INTC(N/2) /2 THEN 3640 

3510 N2=N/2-1 


2640 PRINT " IMPOSSIBILE PERCHE” N E” DISFARI " 
3650 STOP 
3660 END 


Esempio d’esecuzione 


CALCOLO DI UN INTEGRALE DEFINITO CON IL METODO DI SIMFSON 


N INTEGRALE 

(e) +94895226178 
12 +19488106717 
18 +9A882 14551 
24 +9488231313 
30 + 9488295761 


Figura 5.4 
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10 REM CALCOLO DI UN INTEGRALE CON IL METODO DI WEDDLE 

80 PRINT "CALCOLO DI UN INTEGRALE DEFINITO CON IL METODO DI WEDDLE" 
90 PRINT 

100 


FNE CX) =OMEXF(MAXAO LS) 


110 1/8QR(2%*3, 141592653) 
120 1,95 

190 +95 

140 PRINT " N INTEGRALE " 
145 FRINT 


150 FOR N=6 TO 30 STEF 6 
150 GOSUE 4000 

170 FRINT USING 200,N,8 
175 NEXT 
180 ST 

200 34EE IRE AGGGzza Li 
4000 IF N-6XINTC(N/6) 
4010 
4020 
4030 
4040 £ 
4050 $ 
4060 £ 
4090 
4100 
4110 S1=S1+FNEC(X+H) 


0 THEN 4200 


4120 S2=S2+FNE(X+2%H) 
4130 S3+FNE (X+3XH) 
4140 S2=S2+FNE (X+4XH) 
4150 S1=S1+FNE(X+5%H) 
4160 S6=56+FNE (X+6XH) 


4165 XsX+6XH 

4170 NEXT I 

4180 S=,3XHX(FNE(A)-FNE(B)+5xS1+S2+6x53+2xS6) 

4190 RETURN 

4209 FRINT "IMFOSSIBILE FERCHE” "3N3" NON E’ MULTIFLO DI é" 
4210 STOF 

4220 END 


CALCOLO DI UN INTEGRALE DEFINITO CON IL METODO DI WEDDLE 


N INTEGRALE 

é +9501323783 
12 +9488272534 
18 +9488241485 
24 +9488239269 
30 +9488238929 


Figura 5.5 


Confronto dei risultati: Il risultato teorico esatto è 0.95. La tabella seguente mostra i 
risultati ottenuti con ciascuno dei metodi utilizzati. 


‘19485226178 +9501323783 


+9488106717 +9488272534 
+ 9488214551 +9488241485 
+9488231313 +9488239269 
+9488235761 +9488238929 
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Questa tabella mostra che il metodo di SIMPSON, meno preciso di quello di 
WEDDLE, fornisce risultati migliori quando si aumenta almeno per un certo 
tempo il numero dei passi per la durata dei calcoli. 


Invece, i risultati migliori vengono ottenuti con l’altro metodo perN= 6eN= 12 
ed in seguito il risultato si allontana lentamente dal valore esatto. Questo può essere 
dovuto a cause diverse, per esempio per l’influenza degli errori di arrotondamento. 


5.3 CALCOLO DI rx CON IL METODO DEI POLIGONI REGOLARI 


Si possono calcolare dei valori approssimativi di x assimilando il perimetro di un 
poligono regolare alla circonferenza di un cerchio inscritto o circoscritto, questo 
permetterà così di avvicinare con dei valori arrotondati per eccesso o per difetto. 


Il poligono di partenza è quadrato. Partendo dal quadrato si possono calcolare 
due valori di x. Ad ogni iterazione, si darà il numero dei lati del poligono, e si 
effettuerà un nuovo calcolo: 


@ calcolo del lato, 
@ calcolo di n. 


Quando l’intervallo che contiene x non diminuisce più per gli errori di arrotonda- 
mento, si arresteranno le iterazioni. 


Analisi del problema: Si studierà ora il caso del poligono inscritto e di quello 
circoscritto: 


Poligono inscritto: Si considera un cerchio di raggio |. 
il lato del quadrato è dato da: 


AB° = OA? + OB? 
da cui AB= C= 72 


Figura 5.6 
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Si calcola ora la lunghezza di AC del nuovo lato: 


AC? = IA" +IC° con IA = 


IC= OC -0I 
OI è dato da OI = OA° - IA° 
sa 
= 1 7 


Portandolo in AC° si ottiene: 


0° e E 
AC°==—+1+ Za 
gate A 


INIZIO 


STAMPARE 
N, C, P 


Figura 5.7 
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Partendo dalla lunghezza di AC si otterrà n con: 
8-AC=2r 
in modo generale © = cs X. lato del poligono 
Flow-chart provvisorio: Comporta una parte inizializzazione che corrisponde al 
quadrato e in seguito al calcolo iterattivo (Fig.5.$). 


Poligono circoscritto: come prima si parte da un cerchio di raggio 1. Inizialmente il 
lato del quadrato è dato da: 


AB = 2AJ. «con AJ=1 


da cui 


La lunghezza del nuovo lato HK dev'essere determinata. 
Per questo si calcoli: 
IK° = AK° — IA° 


con AK = - — IK poichè IK= KJ 


» 
x 
[5 
(ee) 


Figura 5.8 
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da 4 bdwigie dA vi+S- 1 


riportando in IK° otterremo: 


4 


C° C° 
Lace 


IK°=- c.1K +1k°-2-£ +2 


Figura 5.9 
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da cui la lunghezza del nuovo lato 


Si avvicinerà x connHK = 2 p. 


Flow-chart provvisorio: Con questo metodo di approccio dall’alto, si otterrà un 
nuovo flow-chart che assomiglierà molto al precedente. 


Flow-chart definitivo: Si possono combinare i due flow-chart al fine di ottenere un 
flow-chart che permetta ad ogni iterazione di ridurre l’intervallo di scarto. 


INIZIO 


dove C=V2-V4ZC 


| ara 
dove D spl +D'° — 4] 


Figura 5.10 
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Partendo da questo flow-chart la programmazione non presenta particolari 
difficoltà. L'esperienza mostra che l’influenza degli errori di arrotondamento è 
notevole e che i risultati divergono assai velocemente. Si calcola quindi la differenza 
tra i limiti inferiori e superiori ottenuti da: 


E=Q-—P 


e si arresta il programma quando E diviene negativo. 


Inoltre si può calcolare x mezzi corrispondente alla media aritmetica di P e Q. Si 
constaterà che questa media si avvicina molto velocemente ar. Nonostante i calcoli 
vengano eseguiti con grande precisione, si constata che il valore esatto che comincia 
per 


n: 3,14 1592653589793 


non può essere avvicinato con questo metodo con precisione, poichè la nona cifra 
del valore medio d’intervallo è errata. Questo è dovuto al processo di calcolo che 
determina errori di arrotondamento che procedendo nel calcolo condurranno a 
risultati inesatti. Con un calcolatore avente una precisione inferiore, la precisione di 
tt sarà ancora peggiore. 


CON IL METODO DEL FOLIGONO 
TTTO E DEL POLIGONO CIRCOSCRITTO 
De? 
1 Ped RARE 
{I MD 
Di Me, fog (Fe) 
1 EseP 
PRINT USING 225 
PRINT 
NT USING 280,N,C,D,F,Q;î 
ISF (LD ZERO) ) 
( 150,25) 1) /0 
Fim 0, ERANO 
NxD 
{F4Q) 
OUSING 2890,N,C,D,F,Q,M 
0 THEN 290 
THEN 290 
Q THEN 180 
RO e D FICINF FI-SUF FI-MEZZI 


280 LEERERE ia di zzzi ze e iitii. gua iiiicicizidaztiiiiiizizi: jesi icizizzzizLi 
290 STOF 
300 END 


Figura 5.11 
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Esempio d’esecuzione: 


NUMERO c D PIZINF PI-SUP FT-MEZZI 
4 1,4142136 2,0000000 2,22842712475 4.00000000000  ‘3.41421356237 
8 .7653669 8224271 3.06146745292 3.31370949899. 3,12758797895 
16 3901206 39782247 3.12144515 3.12259727807 3.152024151517 
32 .1960343 1969822 3.1365 355 3.15172490743 3.14413669299 
64 0921359 09825397 3.14033115695 3,14411839524 3.14222477110 
129.0490825 .0490972 3,14127725093 23,14222362994 3.14175044043 
256.0245431 ,0245449 3,14151320115 3.14175036914 3,14163202514 
512.0122718 0122720 3.14157294028 3.14163202047 3,14160251037 
1024 .0061359 .0061359 3.14152772485 3,14160250955 3,14159511721) 
2042 0030680 0030680 3.14159142158 3,141595113937 3.14159326747 
4096 .0015340 0015340 3.14159234176 2,14159326923 3,14159280579 
8192. .0007670 0007670 3,14159256943 2,14159277972 2.14159267458 
16324 .0003235 0003835 3,14159260732 2,14159273259 3,14159267298 
2762 0001917 .0001917 3,14159720921 3,14159186903 2,141592/)9642 

Figura 5.12 


5.4 RISOLUZIONE DI UN'EQUAZIONE CON DICOTOMIA 


Data una equazione f(x) = 0 che ammetta almeno una radice nell’intervallo (4,b). 
Si supponga che la funzione sia continua e limitata in questo intervallo. Si suppon- 
ga inoltre che f/a) e f(b) siano di segno contrario. 


Soluzione: Il sottoprogramma deve poter funzionare correttamente qualunque 
siano le equazioni e i dati forniti. In particolare fornendo due punti A e B in modo 
che f(A)e f(B) siano dello stesso segno, il sottoprogramma dovrà accorgersene e 
fornire la segnalazione. Questo porta dunque a un primo controllo. Nell’esempio 
dato, si utilizza una variabile L alla quale si assegnano tre possibili valori: 


— 1 se f(A)e f(B) sono dello stesso segno, 

— 0 se, dopo una o più iterazioni, f(x) <= EI, 

— 1 se, dopo unao più iterazioni, vi è un intervallo di lunghezza inferiore o uguale a 
E. 


Nell’esempio di soluzione si è supposto che l’utilizzatore abbia provveduto alla 
corretta impostazione dei valori di El e E2 che dovevano essere obbligatoriamente 
di segno positivo. 


Per non modificare nel sottoprogramma i valori di A e B, si utilizzano due 
variabili ausiliarie A1 e B1 che rappresentano i limiti dell’intervallo. Si giunge così 
al flow-chart di Fig.5.13. 


Nel programma di Fig.5.14 la divisione per 2 è stata sostituita da una moltiplica- 
zione per 0.5 che è più rapida (linea 1050). 


L’esempio trattato permette di ottenere il valore X= 7.9999945 quando la radice 
esatta è 8. 
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Algoritmo: A ciascuna iterazione si divide per due la lunghezza dell’intervallo. 


1. Calcolare x = e = b 
y= f(x) 


2. Se f(x) e y sono dello stesso segno: 


porre a = x 
passare in 3 


altrimenti 


porre b = x 
passare in 3 


3. Provare se una delle seguenti condizioni è soddisfatta: 


Iy|ze 
lb-a|<zn 


con € e n dati 


— se nessuna condizione è soddisfatta tornare in l e proseguire le iterazioni, 
— se una condizione è soddisfatta si arrestano le iterazioni. 


Problema: Scrivere una subroutine che può ridurre una equazione secondo il 
metodo esposto. 
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100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
300 
310 
1000 
1020 
1030 
1040 
1050 
1060 


RETURN 


RETURN 


Figura 5.13 


DEF FN FCZ) = (Z + 1) x (Z - 8) 
Ei = 0,0001:E2 = 0,001 

A = 1.6:E = 100 

GOSUE 1000 

IF L= - 1 THEN 300 

PRINT "L= "iL 

FRINT “SOLUZIONI TROVATE X= "3X}" 
FRINT 


= "55 


PRINT "INTERVALLO A = "3A1};" E = "El 
STOF 
FRINT "NESSUNA SOLUZIONE} FC(A) E F(E) STESSO SEGNO" 
STOF 
Ai = A:El = E 
U = FN FC(AI) 
V= FN FC(E1) 
IF U x V = 0 THEN 1190 
X = 0,5 x (Al + El) 
W = FN FX) 


Figura 5.14 (continua) 


1070 IF ABS (W) < = El THEN 1170 
1080 IFUx W « 0 THEN 1120 

1090 fl = X:U= W 

1100 GOTO 1140 

1120 Bi = XIV = W 

1140 IF AES (E1 - A1) > E2 THEN 1050 


1150 L = 1? RETURN 
1170 L = 0? RETURN 
1190 L = - 1? RETURN 
1200 END 


Figura 5.14 fine 


Esempio d’esecuzione 
RUN 
L= 0 
SOLUZIONI TROVATE X= 7,9999939 Y = -5,49014286E-05 
INTERVALLO A = 7,99924317 E = 8,00074464 
BREAK IN 190 


Figura 5.15 


5.5 VALORE NUMERICO DI UN POLINOMIO 


Si cerca di calcolare il valore numerico di un polinomio P(x) di cui si conoscono i 
coefficienti e i valori della variabile. Per questo si utilizza il metodo di Héerner che 
minimizza il numero delle operazioni. Per la valutazione del valore di P(x) dato da: 

P(x) = a x"+ ax + +0,_1X + 4, 
si calcola: 
P:(...((a9x+ a) x + a))x+a;)x+- +0, )X+A, 

Nel programma di Fig.5.16, si è utilizzato il secondo metodo che è più semplice. Il 
sottoprogramma comporta 5 istruzioni (linee 1000-1040) compresa l’istruzione 
RETURN. 


Problema: Costruire un programma che per mezzo di un sottoprogramma valuti 
P(x) con i valori di x dati dal programma principale. 


Soluzione: La formula data significa che sarà necessario fare i seguenti calcoli: 


P=a,x+a, P= a 
Poi : Poi 
P=(P+a,))x fino a P=Px+a, 
Poi : Poi 
P=(P+a,)x P= Px+a, 
e così di seguito fino a: e così di seguito fino a: 
Pie (PE QX P=Px+a, 
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infine: 
P=P+a, 


Da qui si può vedere come può essere utilizzata la seguente iterazione: 


100 REM VALORE NUMERICO DI UN FOLINOMIO SECONDO LO SCHEMA DI HOERNER 
110 DIM AC100) 

115 FRINT "IMFOSTARE IL GRADO DEL FOLINOMIO "; 

120 INFUT N 

130 PRINT "IMPOSTARE I "iN + 1}" COEFFICIENTI IN ORDINE DECRESCENTE" 
140 FORI = 0 TON 

150 INFUT ACI) 

160 NEXT I 

170 FRINT "IMPOSTARE IL VALORE DI X FER IL VALORE DEL FOLINOMIO"} 
190 INFUT X 

200 IF X = 0 THEN STOF 

210 GOSUE 1000 

220 FRINT 

230 FRINT "VALORE DEL FOLINOMIO = "}F 

240 FRINT 

250 GOTO 170 

990 REM VALUTAZIONE DEL FOLINOMIO CON IL METODO DI HOERNER 

1000 F = ACO) 

1010 FORI=1T0N 

1020 F_= Fx X + ACI) 

1030 NEXT I 

1040 RETURN 


1050» (END Figura 5.16 


IMPOSTARE IL GRADO DEL FOLINOMIO ?2 

IMFOSTARE I 3 COEFFICIENTI IN ORDINE DECRESCENTE 
?1 

P1 

BI 


IMFOSTARE IL VALORE DI X FER IL VALORE DEL FOLINOMIO?2 
VALORE DEL FOLINOMIO = 7 

IMFOSTARE IL VALORE DI X FER IL VALORE DEL FOLINOMIO?3 
VALORE DEL FOLINOMIO = 13 

IMFOSTARE IL VALORE DI X FER IL VALORE DEL FOLINOMIO?10 
VALORE DEL FOLINOMIO = ill 

IMFOSTARE IL VALORE DI X FER IL VALORE DEL FOLINOMIO?0 
EREAK IN 200 


Figura 5.17 
5.6 CONCLUSIONI 


Questi esercizi mostrano che la soluzione dei problemi matematici iterativi 
(calcolo di integrali definiti, risoluzione di equazioni, etc.) presenta pochi problemi 
dal punto di vista della programmazione. Del resto il flow-chart è più semplice da 
impostare rispetto agli esercizi di aritmetica presentati nel Capitolo 2. Inoltre, per 
questo tipo di problemi, è necessario poter effettuare un gran numero di operazioni 
di calcolo con precisione e il microcalcolatore si rivela adatto a questo genere di 
calcolo. 


Qualche volta si pone per contro il problema della ‘validità dei risultati‘, poichè 
gli errori di metodo e gli errori di arrotondamento possono avere un’influenza seria 
soprattutto quando i calcoli diventano più complessi. Il lettore interessato a questi 
problemi è invitato a leggere dei testi ‘di analisi numerica“, che trattano metodi di 
calcolo per microcalcolatori e dell’influenza degli errori. 
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CAPITOLO 6 


CALCOLI FINANZIARI 


6.1 PREVISIONE DI CIFRE DI AFFARI 


In questi esercizi si suppone noto il tasso di crescita e si cerca di prevedere 
l’evoluzione della cifra di affari. I due programmi, facili da realizzare, sono rara- 
mente utilizzabili nella pratica, ma possono servire da base alla redazione di 
programmi specifici. 


Primo problema: Un’impresa realizza una data cifra di affari C e prevede un tasso di 
crescita T per gli N anni futuri. Determinare la cifra di affari prevista. I dati sono: 


A = anno d'inizio, 

C = cifra d’affari per l’anno A, 

T = tasso di crescita espresso in percentuale (%), 

N = numero degli anni per i quali si desidera la previsione. 


Esempio: 
A = 1978 
C = 20.000.000 
T=20% 
N=5 


Soluzione: Le sole difficoltà che s’incontrano sono quelle relative alla presentazio- 
ne. Si dovrà tener conto del fatto che il tasso si esprime in percentuale (%) e dà luogo 
a un coefficiente moltiplicativo TI dato da: 


pedi 


100 
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100 PRINT "PREVISIONE DELLE CIFRE D'AFFARI" 
110 PRINT 

120 PRINT "IMPOSTARE ANNO E CIFRA INIZIALE "} 
130 INPUT A,C 

140 PRINT 

150 PRINT "IMPOSTARE TASSO DI CRESCITA "3; 

160 INFUT T 

170 PRINT 

180 PRINT "IMFOSTARE NUMERO ANNI DI FREVISIONE "3; 
190 INFUT N 

200 PRINT 

210 FRINT "ANNI", "CIFRE D'AFFARI" 

220 FRINT 

230 PRINT A,C 

240 Ti s 1 + 0,01 x T 


260 ASA4+1 
270 C=Cx TI 
280 PRINT A,C 
290 NEXT I 
300 END 


Figura 6.1 


CIFRE D'AFFARI 
IMFOSTARE ANNO E CIFRA INIZIALE ?1978,1220 
IMPOSTARE TASSO DI CRESCITA 213 


IMPOSTARE NUMERO ANNI DI FREVISIONE 246 


ANNI CIFRE D'AFFARI 
1978 1220 
1979 1378.6 
1980 1557,818 
1981 1760,33434 
1982 1989,1778 
1983 224777092 
1984 2539 ,98114 
Figura 6.2 


Secondo problema: Si parte questa volta da due dati di base: un volume di vendita e 
una cifra d’affari. Si prevede un aumento in volume di Q% e un rialzo annuale dei 
prezzi di T%. 


Sarà necessario ora prevedere per gli N anni futuri l’evoluzione in volume di 
vendita e l’evoluzione della cifra d’affari. 


I dati sono: 


Anno d'inizio 
Volume 

Cifre d’affari 
Aumento del volume QQ (in percentuale) 
Rialzo dei prezzi T (in percentuale) 
Numero degli anni 


O<d> 
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Soluzione: Si potrà partire dal programma precedente. Si dovrà questa volta, 
tuttavia, tener conto dei due rialzi. 


100 FRINT "IMPOSTARE ANNO, VOLUME E CIFRE D'AFFARI "} 
110 INFUT A,V,C 
120 PRINT "IMFOSTARE L'AUMENTO IN VOLUME E RIALZO DEL PREZZO "; 
130 INFUT Q,T 
140 PRINT "IMPOSTARE NUMERO DEGLI ANNI DI PREVISIONE "} 
150 INFUT N 
140 Qi 14 ,01.xQ 
170 Ti: QLUX (1 + .01x T) 
180 PRINT 
190 PRINT "ANNO", "VOLUME", "CIFRE D'AFFARI" 
200 PRINT 
210 PRINT A,V,C 
220 FORT = 1 TON 
230 A=SA4+1 
240 Va Ux 01 
260 PRINT A,V,C 
270 NEXT I 
280 END 
Figura 6.3 


Esempio di esecuzione: 


POSTARE ANNO, VOLUME E CIFRE 
TARE L/ AUMENTO TN VOLUME 
IMPOSTARE 


D'AFFARI 21978,100,15000 
TALZO DEL FREZZO 25,10 
NUMERO DEGLI ANNI DI FREVISIONE 26 


ANNE VOLUME CIFKE D'AFFARI 
1979 100 15000 

1979 105 17325 

192N 110,25 20010 ,375 

1991 7 23111,9831 
1982 26694,3405 
1982 30831,9633 


1009 


35610,9176 


Figura 6.4 


Nota: In pratica sarà interessante utilizzare le istruzioni PRINT USING e IMAGE 
al fine di ottenere una migliore presentazione dei risultati. 
6.2 RIMBORSO DI UN CAPITALE 

Il rimborso di un capitale può essere fatto in molti modi. Qui vengono presentati 
due metodi relativamente semplici e vengono date delle indicazioni concernenti 
possibili variazioni. 
6.2.1 Primo metodo di rimborso 


Un capitale C è rimborsato in N anni. Alla fine di ogni anno viene rimborsata la 
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stessa frazione del capitale C/N e gli interessi corrispondenti alla frazione del 
capitale che non è stato ancora rimborsato. 


Esempio: capitale C 
numero degli anni N 
tasso T 


alla fine del primo anno si rimborsa: CVN+C XT 
alla fine del secondo anno si rimborsa: C/N + C — C/NXT 
e così di seguito... 


Scrivere un programma che effettui questi calcoli e che stampi anche la somma 
dei rimborsi. 


Soluzione: Si utilizzeranno prima le seguenti variabili: 


R1 = C/N rappresenta la parte del capitale rimborsato ogni anno, 
I importo dell’interesse versato ogni anno, 

R2 importo del rimborso annuale, 

(R2= RI+1) 


Per calcolare I, è necessario conoscere il capitale sul quale vanno conteggiati gli 
interessi. Sia Y questo capitale; alla partenza Y è uguale a C e ogni anno, è 
necessario togliere da Y la quantità R1 rimborsata. Questo conduce allo schema 
seguente: 


T=T/ 100 poichè T è dato in % 


Calcoli preliminari R1= C/N 
Y=C 
I=YXT 

poi per ciascun anno R2=RI+I 
Y=Y_-RI 


Per calcolare la somma totale rimborsata, vi sono due possibili metodi: 


— effettuare il cumulo dei rimborsi: all’inizio Q = 0, 
effettuare il cumulo a ogni rimborso Q = Q + R2 
— effettuare il cumulo degli interessi: all’inizio Q = C, 
a ogni rimborso effettuare il cumulo Q=Q+ I, 


110 


dal punto di vista teorico i due metodi sono identici. Ma in pratica, il primo metodo 
permette di tener conto con maggior precisione degli arrotondamenti. 


LEGGERE C, N,T 


NO 
STAMPARE Q 


Figura 6.5 
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100 REM PROGRAMMA DI CALCOLO ANNUALITÀ’ +. IN 

110 REM QUESTO PROGRAMMA, IN CIASCUN ANNO 

115 REM SI RIMEORSA LA STESSA FRAZIONE DEL CAFITALE, 

120 INPUT "IMFOSTARE CAFITALE, TASSO E NUMERO ANNI ";C,T,N 
130 T= Tx 0.01 

140 Q= CîiY= C 

150 R1 =C/Z/N 

160 PRINT 

170 PRINT "NUMERO INTERESSE  MONTANTE" 

180 FORJ = 1 TON 


190 I=Y=XT 

200 R2 = Ri + I 

210 q=Q+I 

220 Ya Y- RI 

230 PRINT Ji TABC 9);}I; TABC 20);R2 


240 NEXT 4 
250 PRINT ? PRINT "SOMMA TOTALE RIMEORSATA = ";Q 
260 END 


Figura 6.6 


Esempio di esecuzione: 


TMFOSTARE CAFITALE, TASSO E NUMERO ANNI 10000,10,10 
NUMERO INTERESSE MONTANTE 
1 1000 2000 

2 900 1900 

2 800 1800 

4 700 1700 

5 500 1600 

6 500 1500 

7 400 1400 

8 300 1300 

9 200 1200 

10 100 1100 


SOMMA TOTALE RIMEORSATA = 15500 


IMFOSTARE CAFITALE, TASSO E NUMERO ANNI 10000,12,7 
NUMERO INTERESSE MONTANTE 


1 1200 2628,57143 

2 1028,57143 2457 +14286 
3 857,142857 2285,71429 
4 685,714286 2114,28571 
5 514,285714 1942,85714 
È) 342,857143 1771.,42857 
7 171,428571 1600 


SOMMA TOTALE RIMEORSATA = 14800 


Figura 6.7 


6.2.2 Rimborso con mensilità costanti 


Un capitale C è stato prestato a un tasso annuale T1. Ogni mese si rimborsa la 
stessa somma. Sapendo che con N mensilità si deve rimborsare tutto, calcolare 
l’importo della mensilità e su domanda dell’utilizzatore la ripartizione del capitale 
rimborsato, e di quanto pagato per gli interessi. 
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Per questo si opererà nel seguente modo: 


— calcolo del tasso equivalente T1 (corrisponde al tasso mensile). Questo tasso TI 
è definito con (1+T1) =1+T. 


da cui, estrapolando, si ricava: 
TI=(1+T)-1 


se T è dato in percentuale, sarà necessario a priori dividerlo per cento. 
E° da notare che sovente le banche applicano una formula che è loro più 
favorevole: 


T1 mensile = T annuale 
12 N 
Xi ARRE = THKI+TI) 
— calcolo dell’importo della mensilità per M = C d+TD°=1 


— calcolo eventuale del dettaglio. 
La parte del capitale rimborsata si ottiene con un ragionamento molto semplice: 


per la prima mensilità, la quota interesse è: C X TI, 

dunque il capitale rimborsato è: M — C X TI, 

il capitale restante C — (M — C X TI) servirà per il calcolo degli interessi della 
seconda mensilità e così di seguito. 


Enunciazione del problema: Scrivere un programma che: 


— legga i dati: 
— capitale, 
— tasso annuale espresso in percentuale 
— numero delle mensilità; 


— faccia i calcoli e stampi: 
— tasso equivalente, 
— importo delle mensilità, 
— somma totale rimborsata; 


— domandi all’utilizzatore se desidera il dettaglio dei rimborsi e se si, inumeri della 
prima e dell’ultima mensilità per le quali si desidera il dettaglio. 


Soluzione: L’inizio del programma deriva dall’enunciazione e dalle formule date a 
condizione di: 


— leggere T in percentuale, 
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— calcolare T = T/100, 
— verificare il valore T1] ottenuto, ma stampare 100 X T1 per darlo in 
percentuale. 
Per la seconda parte verrà disegnato un flow-chart nel quale A e B rappresentano 
i numeri rispettivamente della prima e dell’ultima mensilità per le quali si desidera- 
no i dettagli. 


Ri=T1xC1 
R2=M_- RI 


VISUALIZ= 
ZAZIONE 
DI K R1, R2 


Figura 6.8 


CI rappresenta la parte del capitale non ancora rimborsato. In questo flow-chart 
(Fig.6.8), è stato previsto un loop per K da 1 a N per eventuali estensioni. Se ci si 
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limita al problema esposto è sufficiente prevedere un loop per K che varia da 1 a B, il 
che permette di sopprimere il test K > B. 


100 
110 
12 

130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
230 
340 
350 


REM CALCOLO DELLE MENSILITA’ FER RIMBORSARE UN CAFITALE 
INFUT "MONTANTE DEL CAPITALE PRESTATO? ";C 
INFUT "TASSO ANNUALE IN % 13) 
INFUT "NUMERO DELLE MENSILITA” 1 "iN 
#31 TZ IO0)A CI ZI) = 1 
* TL / (1- (1+ T1) 4 ( - N)) 


"TASSO EQUIVALENTE = "3100 x Ti;" %" 
"MONTANTE DI OGNI MENSILITÀ” = "M 
"SOMMA TOTALE RIMEORSATA = "MAN 


TINFUT "VOLETE IL DETTAGLIO DI ALCUNE MENSILITÀ‘ ? "3R$ 
IFOR$ «2» "SI" THEN STOF 
INFUT "NUMERI DELLE MENSILITÀ? ? "3W,E 
PRINT "NUMERO", "INTERESSE", "AMMORTAMENTO" 
PRINT 3îC1 = € 
FOR K_= 
IFK a 

Ri = TL * C13R2 M- RI 
IFK < W THEN 330 

PRINT K,R1,R2 

Cis C1 - R2 

NEXT K 
END 


Figura 6.9 


Esempio di esecuzione: 


MONTANTE DEL CAFITALE FRESTATO: 
TASSO ANNUALE IN % : 
NUMERO DELLE MENSILITA’ 3 
TASSO EQUIVALENTE = .486755045 % 
MONTANTE DI OGNI MENSILITA’ = 322,438624 
SOMMA TOTALE RIMEORSATA = 46431,1619 
VOLETE IL DETTAGLIO DI ALCUNE MENSILITA’ ? SI 

NUMERI DELLE MENSILITA’” ? 50,55 


3322 


NUMERO INTERESSE AMMORTAMENTO 
50 119,151941 203,286683 

Si 118,162433 204,276191 

52 117.168108 205,270516 

53 116,168943 206,269681 

54 115,164915 207,273709 

55 114,156 208,282624 
MONTANTE DEL CAFITALE PRESTATO! 60000 

TASSO ANNUALE IN % * 8.5 

NUMERO DELLE MENSILITA” * 180 


TASSO EQUIVALENTE = .682149362 % 

MONTANTE DI OGNI MENSILITA’ = 579845213 
SOMMA TOTALE RIMBORSATA = 104372,138 
VOLETE IL DETTAGLIO DI ALCUNE MENSILITÀ’ ? NO 

BREAK IN 240 


MONTANTE DEL CAPITALE FRESTATO: 30000 
TASSO ANNUALE IN % 
NUMERO DELLE MENSILITÀ” 
TASSO EQUIVALENTE = .643403037 % 

MONTANTE DI OGNI MENSILITA’ = 359,57257 
SOMMA TOTALE RIMEORSATA = 43148,7084 
VOLETE IL DETTAGLIO DI ALCUNE MENSILITÀ‘ ? NO 

BREAK IN 240 


8 
120 


+ 00 


Figura 6.10 
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6.3 CALCOLO DEL TASSO DI CRESCITA 


Questo esercizio non dev'essere affrontato se non prima dell’esercizio sul calcolo 
della regressione lineare (9.3). Si conoscono le cifre d’affari degli anni precedenti di 
un’impresa. Si pensa che la cifra d’affari segua una legge matematica del tipo: 


C(1 + 7) 


t rappresenta il tasso di crescita e i l’anno corrente. 


Il problema è determinare C e 1 di prevedere le cifre d’affari per gli anni futuri. Si 
potrà limitare la previsione a 5 anni. 


Analisi matematica: Per risolvere questo problema si adotterà un artificio sul piano 
matematico al fine di renderlo più facilmente risolvibile. Si cercherà di determinare 
C e tin modo da minimizzare non: 


DA 


Q = Y (Log[c(1 + 7)]- Log y,Y 
= Y (Log c + i Log(1+7)- Log y,Y 


poniamo: 


Log C = B Log y; = 2, 
Log (1+7)=A 


E° necessario minimizzare la quantità Q data con: 


Q= D(B+iA-z,) 


Questo ci riporta al problema della regressione lineare esposta nel Capitolo 9.3. 


Problema: Partendo dal programma dato nel Capitolo 9.3, costruire un programma 

che calcoli il tasso di crescita t ed effettui una previsione per i cinque anni futuri. 
In questo programma 1 sarà rappresentato dalla variabile TO. Si supporrà che gli 

anni siano immagazzinati in una tabella T e le cifre d’affari in una tabella X. 


Soluzione: Questo programma comporta tre parti distinte: 


— lettura dei dati N, tabelle T(I) e X(I) e costruzione della tabella Y come: 
Y(1) = Log (X(1)) 
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— richiamo della subroutine che effettua la regressione lineare e calcola dei 
coefficienti C e TO. Partendo da A e B dati dalla subroutine, si ottengono 
questi coefficienti: 

C=eB et TO= 4-1 


— stampa di TO e dei risultati: 
— Per ogni anno conosciuto, cifra d’affari vera e cifra d’affari stimata. 
— Per ognuno dei cinque anni futuri, unicamente la cifra d’affari stimata. 


Per evitare la stampa dei decimali inutili si rimpiazzerà la stampa delle cifre 
d’affari stimate Z con: 


INT (100 X Z) / 100 


Il flow-chart generale riportato in Fig.6.11, è quindi molto semplice. 


LEGGERE I DATI PASSAGGIO 
IN LOGARITMI 
DI RIGRESSIONE 
LINEARE 
VISUALIZZAZIONE 
DI C.A. REALI e C.A. STIMATI 


VISUALIZZAZIONE 
DELLA PREVISIONE 


LINEE 110-150 


LINEE 155:+180 e 1.000-+1.140 


LINEE 210-280 


LINEE 300-340 


Figura 6.11 


Il programma di Fig.6.12 costituisce un esempio di quanto si può scrivere, sarà 
possibile migliorarlo stampando: 


— coefficiente di correlazione 


— un intervallo stima per la previsione che sarà utile ma complicherà il 
programma. 
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Attenzione: Una tale previsione non può essere utilizzata in pratica senza precauzio- 
ni. In effetti la cifra d’affari di un’impresa dipende da numerosi fattori in modo 
particolare, dalla congiuntura economica e dalla concorrenza, e questi fattori 
possono influenzare notevolmente le realizzazioni future. 

100 DIM TC(15),X(15),Y(15) 

110 READ N 


120 FORI = 1 T0N 
130 READ TCI) XI) 


140 YCI) = LOG (XCI)) 
150 NEXT I 

155 NO = T(1) 

1600 GOSUE 1000 

170 (XP (E) 

180 EXF CA) — 1 


190 PRINT "TASSO STIMATO = "3100 x TO 
200 FRINT 
210 T2 = 1 + TO 

220 PRINT "ANNO", "ATTUALE", "PREVISTO" 

2230 PRINT 

240 Z= € 

260 FOR I = 1 TON 

280 PRINT TCI),XCI), INT (100 x Z) / 100 
2985 Z= Z% T2 
290 NEXT I 

300 FOR I = 1T05 
310 T3= + I 

320 Z Zx T2 

330 PRINT T3, INT (100 x 2) / 100 

340 NEXT I 

400 DATA 6 

410 DATA 1975,99,2,1976,110,1977,121.3 
420 DATA 1978,133,1,1979,146,3 

430 DATA 1980,160 

500 STOF 

1000 Ul = 0 

1010 42 = 0 

VI = 03V2 = u 


50 FOR I = 1 TON 

3. TA = TI) - NO 

U1 = Ul + T4 

1070 VI = VI + Y(I) 

1080 U2 = UZ + T4 x T4 
1090 V2 = V2 + YI) * YI) 
1100 Wes W+ TA * YL) 
1110 NEXT I 

1120 A so (We UL x VI Z N) / CU? UL Xx UL / N) 
1130 BE = (VI - A x Ul) / N 
1140 RETURN 

1150 END 


Figura 6.12 
Esempio di esecuzione: 


TASSO STIMATO = 10,0084650 


ANNO ATTUALE PREVISTO 
1975 99,2 

1976 110 

1977 121,3 

1978 123.1 

1979 146.3 

1980 160 

1981 

1982 

1983 35 

1984 258.96 

1985 284,98 : 

BREAK IN 500 Figura 6.13 
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6.4 CALCOLO DELL'AMMONTARE DELLE IMPOSTE 


Supponiamo di voler calcolare a quanto ammontano le tasse dovute da un 
semplice salariato e supponiamo che la legge fiscale preveda quanto segue: 


— Determinazione del reddito imponibile poi calcolo del quoziente familia- 
re. Questo reddito imponibile è definito da R = 0,72 X reddito netto + 
pensione — deduzioni. 
Il quoziente familiare è dato da: Q = R/N, con N il numero dei componenti la 
famiglia. 


— Calcolo dell’importo d’imposta; sono possibili due metodi: 
— l’applicazione di una formula diretta, 
— il calcolo per passi successivi utilizzando la tecnica degli scaglioni tariffari, 
— calcolo dell’imposta relativa ad ogni scaglione 
— cumulo 
— poi moltiplicazione per il numero dei componenti 


— Si deducono in seguito dall’importo lordo alcune eventuali voci per ottenere 
l’importo netto. 


Per i redditi incassati nel corso dell’anno la tabella di Fig.6.14 indica gli ipotetici 
scaglioni tariffari. 


REDDITO IN MILIONI ALIQUOTE 


7.925 
8.300 
9.925 
15.700 
20.625 
25.925 


31.350 
36.175 
62.600 
86.125 
105.950 
125.050 
oltre 


Figura 6.14 


Esempio: Supponiamo che una famiglia composta da tre persone abbia incassato 
redditi per 130.000.000 di lire. Questa famiglia può dedurre 5.000.000 di lire 
rappresentati dagli interessi di un prestito contratto per l’acquisto di una casa. 
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Reddito imponibile: 0,72 X 130.000.000 — 5.000.000 = 88.600.000 
Quoziente familiare: 88.600.000/3 = 29.533.000 
Importo d’imposta: 


( 8.300.000 - 7.925.000) X 5/100= 18.750 
( 9.925.000 - 8.300.000) Xx 10 / 100 = 162.500 
(15.700.000 - 9.925.000) Xx 15 / 100 = 866.250 
(20.625.000 - 15.700.000) X 20 / 100 = 985.000 
(25.925.000 - 20.625.000) X 25 / 100 = 1.325.000 
(29.533.000 - 25.925.000) X 30 / 100 = 1.082.400 


TORALE niro ata 4.439.900 
Questo importo moltiplicato per i componenti dà: 4.439.900 X 3 = 13.319.700 


Da questa somma si potrà dedurre un eventuale “avere fiscale‘. Si constati che 
questa famiglia raggiunge lo scaglione tariffario del 30%. 


Problema: Il fine di questo esercizio è evidentemente la scrittura di un programma 
che, partendo da alcuni dati, calcoli l’importo dell’imposta sul reddito. Si propone 
di seguire i passi sotto indicati: 


— precisare il metodo in modo da ridurre il numero degli scaglioni (scaglione 
con tasso 0) al fine di ottenere 12 scaglioni utili, 


— costruire in seguito il flow-chart di dettaglio unicamente per la parte 
calcolo dell’imposta, 


— scrivere un programma che chiama un sottoprogramma per la lettura 
degli scaglioni e delle percentuali e un altro sottoprogramma per il calcolo 
del montante d’imposta. 


Soluzione: Siano F e T due tabelle contenenti rispettivamente i valori “più bassi‘ 
degli scaglioni tariffari e le relative aliquote. Non si è interessati allo scaglione 0% 
quindi si avrà: 


F(1) = 7 925 F(12) = 125 050 
T(1) = 0,05 T(12) = 0,6 


Se Q < F(1) nessuna imposta, 
se F(1) < Q < F(2) si calcolerà S= (Q — F(1)) X T(1) 
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— in generale, se F(K—1) < Q < F(K) si calcolerà: 


S=(Q-F(K-1)xT(K-1)+ (F)- F(a-1)xT1- 1) 


Questo porta alla parte di flow-chart in Fig.6.15. 


S=S+(F(1) FI 1)x TI - 1) 
I=1+1 


S=S+(Q- F(1-1))x TI - 1) 


NO CON M “PIATTAFORMA” 


S=S+(Q- T(M)) x T(M) 


S = INTERO (SxN+0.49) 


RETURN 


E’ necessario ora scrivere un programma che sia facilmente adattabile alla 
modifica degli scaglioni: 


PER EFFETTUARE 
L'ARROTONDAMENTO 


Figura 6.15 


— il sottoprogramma di lettura (linee 1000 - 1170) legge il numero degli 
scaglioni M poi le coppie F(I) e T(1) 
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— il sottoprogramma di calcolo (linee 1200 - 1290) corrisponde al flow-chart 


di 


Fig.6.15. 


Il programma principale è quindi molto corto: pone le domande necessarie e 
verifica che i dati impostati siano corretti. Calcola in seguito il reddito imponibile e 
il quoziente familiare poi chiama il sottoprogramma per il calcolo delle tasse. 


80 D 
90 G 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
205 
210 
220 
230 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 


IM F(15),T(15) 

OSUE 1000 

INFUT "IMFOSTARE IL REDDITO NETTO ";R1 

IF Ri « 0 THEN 100 

INFUT "IMFOSTARE IL MONTANTE DELLA FENSIONE "}F 
IF FP € 0 THEN 120 

INFUT "IMFOSTARE IL MONTANTE DELLE DEDUZIONI ";D 
IF D £ 0 THEN 140 

R = .72 x RL +F - D 

INFUT "IMFOSTARE IL N. DEI COMFONENTI "}N 


IFUN < 0 ORN< > INT (N) THEN 170 
Q=RZN 

FRINT ! FRINT "REDDITO IMPONIEILE = "3R 

PRINT "QUOZIENTE FAMILIARE = "5Q 

GOSUE 1200 

FRINT "MONTANTE DELLE IMPOSTE = "38 

STOF 

REM SOTTOFROGRAMMA DI LETTURA DEGLI SCAGLIONI E DEI TASSI 
READ M 


FORI =1 TOM 

READ FCI), ASTCI) = .01 x A 

NEXT I 

RETURN 
DATA 12 
DATA 7925,5 
DATA 8300,10 
DATA 9925,15 
DATA 15700,20 
DATA 20625,25 
DATA 25925,30 
DATA  31350,35 
DATA 36175,40 
DATA 62600,45 

DATA 86125,50 

DATA  105950,55 
DATA 25050, 60 

S = 0} IFQ = F(1) THEN RETURN 
FORI=2TO0M-1 

IF Q « FCI) THEN 1270 

S..= 8 + (FI) - F(I — 1)) x TI — 1) 
NEXT I 

S = S + (Q — F(M)) x TCM) 

GOTO 1280 

S = St+ (Q- FCI — 1)) * TI - 1) 
S = INT (N x S + ,49) 

RETURN 


Figura 6.16 


Esempio di esecuzione: 
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IMFO 
IMFO 
IMFO 
IMFO 
REDD 
QUOZ 
MONT 
EREA 


STARE IL REDDITO NETTO 1000000 
STARE IL MONTANTE DELLA FENSIONE 0 
STARE IL MONTANTE DELLE DEDUZIONI 0 
STARE IL N. DEI COMFONENTI 1 

ITO IMPONIBILE = 720000 

TENTE FAMILIARE = 720000 

ANTE DELLE IMFOSTE = 394712 

K IN 230 


Figura 6.17 


6.5 INCIDENZA DI UNA REDDITO SUL POTERE D’ACQUISTO 


Il potere di acquisto di una famiglia può essere definito dal reddito netto meno le 
spese fisse imposte (spese professionali, rimborsi di prestito per l’acquisto di una 
casa, imposte da pagare). Poichè in Italia le imposte sono pagate con un ritardo 
medio di un anno, si può considerare la loro incidenza sul potere d’acquisto 
leggermente più bassa poichè in questo periodo è possibile accumulare degli 
interessi, è sufficiente porre: 

C 


—— con U è il tasso d’interesse 
1+U 


per poter pagare con un anno di ritardo la somma C. Il tasso U dipende dall’utiliz- 
zatore. 


Problema: Completare il programma precedente in modo che venga stampato il 
tasso dello scaglione più elevato per il quale il contribuente deve pagare delle 
imposte e poi che, a partire da un reddito addizionale netto (supposto senza spese) e 
da un tasso di interesse, calcoli l'aumento reale del potere di acquisto. 

Per questo esercizio, si supporrà che il supplemento d’imposta è pagato in media 
un anno dopo aver incassato il reddito addizionale. 

Si deve riflettere sulle condizioni di validità del programma e su altre modifiche 
per migliorare le condizioni di validità. 


Soluzione: Il flow-chart generale del programma è molto “lineare‘*, poichè non vi 
sono loop di calcolo e “calcoli iterativi‘* etc. 


Poichè si farà due volte il calcolo dell’importo delle imposte, è consigliabile 
chiamare in questo caso una subroutine di calcolo. 


Partendo dal programma precedente è sufficiente aggiungere una parte che legga 
il reddito addizionale, lo incorpori nel reddito precedente al fine di permettere un 
nuovo calcolo delle imposte dovute.In seguito, si calcoli l'interesse rappresentato 
dalla differenza tra i due importi che si suppone investiti durante l’anno, al tasso 
dato U dall’utilizzatore. 


Si può infine calcolare il potere d’acquisto addizionale che è uguale a: 
reddito supplementare — supplemento di imposta + interesse 


Critica: In questo programma si suppone che il reddito addizionale non dia luogo a 
delle spese supplementari non deducibili (esempio: asilo d’infazia, riduzione del- 
l'ammontare degli assegni familiari, etc.). Se è così, sarà sufficiente dedurre dal 
potere d’acquisto addizionale, dato dal programma, l’ammontare delle spese sup- 
plementari o la riduzione degli assegni, per ottenere l’aumento reale del potere 
d’acquisto. 
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10 REM PROGRAMMA DI CALCOLO DEL MONTANTE DELL’IMFOSTA SUL REDDITO 
20 REM DELLE FERSONO FISICHE 

30 REM IL SOTTOFROGRAMMA LINEE 990-1170 LEGGE GLI SCAGLIONI E I TASSI 
40 REM E’ NECESSARIO MODIFICARLO QUANDO CAMEIA L’ANNO 

80 DIM F(15),T(15) 

90 GOSUE 1000? FRINT "IMFOSTA SUL REDDITO DELL'ANNO "}:A9: FRINT 
100 INFUT "IMPOSTARE IL REDDITO NETTO ";R1 

110 IF Ri « 0 THEN 100 

120 INFUT "IMFOSTARE IL MONTANTE DELLA FENSIONE ";F 

130 IF FP < 0 THEN 120 

140 INFUT "IMFOSTARE IL MONTANTE DELLE DEDUZIONI "3}D 

150 IF DD 0 THEN 140 

160 R = +72 X R1 + F - DîR2 = R 

170 INFUT "IMFOSTARE IL N, DEI COMPONENTI "3N 


180 IFN< 0 0RN< © INT (N) THEN 170 
190Q0=R/ZN 

200 PRINT : PRINT "REDDITO IMPONIBILE = ";R 
205 PRINT "QUOZIENTE FAMILIARE = "30Q 


210 GOSUE: 1200 

220 FRINT "ENTRATE NELLO SCAGLIONE AL "3T1;" %" 
230 PRINT "IL MONTANTE DELLE VS, IMPOSTE E’ "365 
240 L = Ri + P_- D- S°SI1 = S 

250 FRINT "DISFONETE DI "}L3" FER VIVERE" 


260 FRINT 

270 INFUT "QUAL/E/ IL MONTANTE DI UN REDDITO ADDIZIONALE "A 
275 IF A £ = 0 THEN 500 

280 R = R2 + .72 *X A:Q=R/ZN 

290 FRINT "NUOVO REDDITO IMPONIBILE = "3R 

300 FRINT "NUOVO QUOZIENTE FAMILIARE = "3Q 

310 GOSUE 1200 

320 FRINT "NUOVO MONTANTE DELLE IMFOSTE = "35 


325 FRINT "A QUALE TASSO FOTETE FIAZZARE LA SOMMA CHE RAFFRESENTA IL" 
330 INFUT "SUFFLEMENTO DI IMFOSTA CHE FAGHERETE L'ANNO SUCCESS, "3U 
335 Li = A + S1 - S + 0.01 x U x (S + S1) 

340 PRINT "DISFONETTE INFATTI DI UN FOTERE DI ACQUISTO ADD. DI "; 

345 PRINT Li 

500 INFUT "ALTRO CASO DA TRATTARE "3R$ 


510 IF R$ < 2 "SI" THEN STOF 
520 FRINT ? GOTO 100 
980 STOF 


990 REM SOTTOFROGRAMMA DI LETTURA DEGLI SCAGLIONI E DEI TASSI 
1000 READ M,A9 

1010 FOR I = 1 TOM 

1020 READ FCI), AI:TCI) = .01 x A 
1030 NEXT I 

1040 RETURN 

1050 DATA 12,1978 

1060 DATA 7925,5,8300,10 

1080 DATA 9925,15,15700,20 
1100 DATA 20625,25,25925,30 
1120 DATA 31350,35,36175,40 
1140 DATA 62600,45,86125,50 
1160 DATA 105950,55,125050,60 


1200 S = 0:T1 = 0? IFQ< = F(1) THEN RETURN 
1210 FORI =2 TOM- 1 

1220 IF Q « F(I) THEN 1270 

1230 S = S + (FCI) — FCI — 1)) x TCI - 1) 

1240 NEXT I 

1250 S = S + (Q - F(M)) x TCM) 

1260 GOTO 1280 

1270 S= S+ (Q- FCI — 1)) x TCI - 1) 

1280 S = INT (N x S + +49):T1 = 100 x TC(I - i) 
1290 RETURN 


Figura 6.18 


Nota: In questo programma, la barra obliqua / separa più istruzioni sulla stessa 
linea. 
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Esempio di esecuzione: 


IMFOSTA SUL REDDITO DELL'ANNO 1978 

IMPOSTARE IL REDDITO NETTO 10000000 

IMFOSTARE IL MONTANTE DELLA FENSIONE 500000 

IMFOSTARE IL MONTANTE DELLE DEDUZIONI 180000 

IMFOSTARE IL N. DEI COMPONENTI 4 

REDDITO IMFONIEILE = 7520000 

QUOZIENTE FAMILIARE = 1880000 

ENTRATE NELLO SCAGLIONE AL 55% 

IL MONTANTE DELLE VS, IMFOSTE E’ 4362850 

DISFONETE DI 5957150 FER VIVERE 

QUALE” IL MONTANTE DI UN REDDITO ADDIZIONALE 3000000 
NUOVO REDDITO IMFONIEILE = 9680000 

NUOVO QUOZIENTE FAMILIARE = 2420000 

NUOVO MONTANTE DELLE IMFOSTE = 5658850 

A QUALE TASSO FOTETE FIAZZARE LA SOMMA CHE RAFFRESENTA IL 
SUFFLEMENTO DI IMFOSTA CHE FAGHERETE L'ANNO SUCCESS, SI 
PREENTER 

SUFFLEMENTO DI IMFOSTA CHE FAGHERETE L'ANNO SUCCESS, 1650000 
DISFONETTE INFATTI DI UN FOTERE DI ACQUISTO ADD. DI 2,1385704E+10 
ALTRO CASO DA TRATTARE 1210000 

BREAK IN 510 


Figura 6.19 
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CAPITOLO 7 


GIOCHI 


7.0 INTRODUZIONE 


I giochi intellettuali che hanno luogo tra l’essere umano e il calcolatore interessa- 
no molte persone: i giovani per l’aspetto ‘“gioco‘, i meno giovani per l’aspetto 
competizione contro questa macchina così rapida e gli specialisti che fanno ricerche 
sulle strategie. 


L’esperienza mostra che la scrittura dei programmi di giochi costituisce, in 
generale, un lavoro lungo e difficile che non è, salvo casi eccezionali, alla portata del 
principiante. Questo è evidentemente il caso del gioco degli scacchi, ma per dei 
giochi più semplici si constata che il programma è piuttosto lungo e che le prestazio- 
ni del calcolatore non sono sempre buone. 


Esistono dei giochi la cui programmazione è molto facile, o perchè non vi sono 
strategie da implementare (esempio il gioco alto-basso), o perchè la strategia è 
molto semplice (esempio il gioco di Marienbad). 


Nel momento in cui le strategie si complicano il programma inevitabilmente si 
complica. 


I tre giochi qui presentati sono particolarmente facili dal punto di vista della 
programmazione. 


7.1 UN GIOCO: ALTO-BASSO 


Prima parte: Bisogna trovare un numero intero N posto tra 0 e A, scelto a caso dal 
calcolatore. Quanto s’imposta un numero X, il calcolatore lo confronta con il 
numero da trovare N. 


se X= N. il calcolatore stampa INDOVINATO DOPO I TENTATIVI, 
dove I rappresenta il numero dei tentativi effettuati dal giocato- 
re, 

se X<N. il calcolatore stampa E’ ALTO 

se X>N. il calcolatore stampa E’ BASSO 
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Analisi: Per questo problema sarà necessario utilizzare 4 variabili: 


A limite superiore dell’intervallo 
N numero da trovare 

X numero impostato dal giocatore 
I numero di tentativi effettuati 


La variabile A non è indispensabile, ma il suo utilizzo riduce il numero delle linee 
da modificare se si desidera cambiare il limite superiore dell’intervallo. 

La variabile I è in effetti un ““contatore‘ che si incrementa di una unità ad ogni 
nuovo tentativo. 


Flow-chart: Tra tutti i metodi possibili, il flow-chart seguente sembra essere uno dei 
più semplici. 


DEBUT 


LEGGERE A 


VISUALIZZARE “DOVETE 


TRAOE A” 


ESTRARRE N CASUALE 
IN (0,A)I = O 


FINE 


Figura 7.1 
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Per la programmazione, bisogna scegliere N intero casuale nell’intervallo (0, A). 
Per questo bisogna scrivere N = INT (A X RND(X). Questo metodo di scrittura 
può cambiare in minima parte da un sistema all’altro. 


100 REM CALCOLO DELLE MENSILITÀ‘ FER RIMEORSARE UN CAFITALE 
110 INFUT "MONTANTE DEL CAFITALE FRESTATO? ";C 

120 INFUT "TASSO ANNUALE IN % + "3T 

130 INFUT "NUMERO DELLE MENSILITÀ‘ 1 "3N 

140 T1 = (1 + T/ 100) 4 (1 / 12) - 1 

150 M= C x T1 / (1 - C1 + T1) 4 ( - N)) 


160 FRINT 

170 FRINT "TASSO EQUIVALENTE = "3100 x T1lj" %" 

180 PRINT 

190 FRINT "MONTANTE DI OGNI MENSILITA” = "} TAEC 35);M 

200 PRINT 

210 FRINT "SOMMA TOTALE RIMEORSATA = "3 TAEC 33);M X N 
220 PRINT 


230 INFUT "VOLETE IL DETTAGLIO DI ALCUNE MENSILITA” ? "R$ 

240 IF R$ < > "SI" THEN STOF 

250 INFUT "NUMERI DELLE MENSILITA/ ? "3;W,E 

260 PRINT ? PRINT "NUMERO "3 TABEC 8);" INTERESSE"} TABC 23); "AMMORTAMENTO 


Figura 7.2 


Seconda parte: Giocando ci si accorge che bisogna fare un piccolo sforzo per 
ricordare in quale intervallo si trova il numero, poichè, normalmente, questo 
intervallo diventa più piccolo. Per aiutare il giocatore, il calcolatore può ricordar- 
gli, dopo ogni tentativo errato, quali sono i limiti del nuovo intervallo. 


Analisi: Per questo miglioramento saranno necessarie due nuove variabili C e D che 
rappresentano i limiti dell’intervallo. 


Inizialmente C=0e D= A, 
dopo un tentativo errato: 


se X< N porre C = MAX (C,X) 
se X>N porre D= MIN (D,X) 


rigorosamente, se si ha la certezza che il giocatore non imposterà mai un numero 
superiore all’intervallo, sarà sufficiente porre C = X e D= X. 


L’ordine di stampa diventa lo stesso in ogni caso errato poichè è sufficiente dare i 
limiti di intervallo. 


Flow-chart: E’ di immediata realizzazione sulla base del precedente. 
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INIZIO 
LEGGERE A 


VISUALIZZARE 
“TROVARE UN 
NUMERO TRA 0EA 


ESTRARRE N CASUALE 
NELL'INTERVALLO (0,A) 
I=0 C=0 D=A 


Terza parte: Un buon metodo per raggiungere rapidamente il numero cercato è 
quello di provare, a ogni tentativo, un numero che si trovi nel mezzo dell’intervallo. 


Piuttosto che fare questo calcolo a mente è possibile modificare il programma per 
rendere automatica questa operazione. 


Analisi: E’ sufficiente nell’istruzione di stampa aggiungere il valore di (C + D)/2. 
Non vi è dunque alcuna modifica da apportare al flow-chart. 
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Programma: 


10 INFUT "IMFOSTARE IL NUMERO MASSIMO "3A 

20 PRINT 

300 FRINT "DOVETE TROVARE UN NUMERO COMPRESO TRA 0 E "3A 
40 PRINT 

500 PRINT "IMPOSTARE IL TENTATIVO " 

60 N=INTCAXRNDCOX) ) 

70 Ist? C=0 : D=sA 


800 Ist+1 

90 INFUT X 

100 XF XeN THEN GOTO 190 

1109 IF XEN THEN GOTO 112 ELSE GOTO 116 


112 IF XSC THEN Cex 3 GOTO 115 
116 IF XD THEN D=x 


118 PRINT "TRA "3C5" E "3Di" META” = "5 (0+D)/2 
120 GOTO 80 
190 PRINT "INDOVINATO DOFO ";1}" TENTATIVI" 
140 END 

Figura 7.4 


Quarta parte: Partendo dal momento in cui il giocatore introduce il numero 
suggerito dal programma, si può far giocare la macchina contro se stessa, in questo 
caso il giocatore diviene spettatore. 


Analisi: Il flow-chart sarà di poco modificato, l’istruzione di lettura di X sarà 
rimpiazzata da X = INT((C + D ) / 2) o da 


X= INT(C+D4+1)/2) 


e l’istruzione che stampa i limiti dell’intervallo potrà essere soppressa. Invece se si 
desiderano conoscere i tentativi effettuati dalla macchina, sarà necessario inserire 
un’istruzione di stampa X. 


Commento: In questo programma utilizzando il video la visualizzazione si effettua 
con una cadenza molto rapida. Per dare all’utilizzatore il tempo di leggere, è stata 
aggiunta l’istruzione SLEEP 5, che sospende l’esecuzione per 5 secondi dopo ogni 
gioco. Questa istruzione, SLEEP, non è disponibile per tutti i sistemi. Si può 
ottenere lo stesso risultato in altri modi: 


— con l’istruzione WAIT quando disponibile, 


— inserendo un loop di calcolo che impedisce all’unità centrale di passare da 
un tentativo all’altro troppo rapidamente. 
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10 INFUT "IMFOSTARE Il. NUMERO MASSIMO "3A 


20 PRINT 
30 PRINT "DOVETE TROVARE UN NUMERO COMPRESO TRA 0 E ";A 
40 PRINT 


45 SLEEP 5 
50 PRINT "IMPOSTARE IL TENTATIVO " 

60 N=INTCAXRND(X)) 

70 I=0 3 C=0 3: D=A 

80 I=t+1 

90 X=INT((C4D)/2) 

100 IF X=N THEN GOTO 130 

110 IF X€N THEN GOTO 112 ELSE GOTO 116 
112 IF X5C THEN CeX 3 GOTO 118 

116 IF XD THEN 


118 PRINT "TRA "Cp" E o "3D}" META? = "3 (C+D)/2 
120 GOTO 80 
130 PRINT "INDOVINATO DOFO "313" TENTATIVI" 
135 GOTO 45 
140 END 

Figura 7.5 


7.2 TROVA IL NUMERO 


Questo gioco, che deriva da quello precedente, consiste nel trovare un numero 
scelto a caso dal computer dando una serie di intervalli (due numeri) che devono 
contenere il numero stesso. Il calcolatore dice: 


— se il numero è stato inquadrato, 
— se è più piccolo, 
— se è più grande. 
Per esempio se il numero da trovare fosse 55 e fossero proposti i numeri 18,24, il 
computer risponderebbe “PIU” GRANDE“. 


Primo problema: Studiare un semplice programma che permetta di giocare a questo 
gioco e che conti il numero dei tentativi effettuati dal giocatore. 


Secondo problema: Proporre un programma più sofisticato che sia capace di capire 
se il giocatore gioca bene. 


Soluzione al primo problema: Saranno necessarie le seguenti variabili: 


A = numero massimo dell’intervallo (1,A), 

N = numero da trovare, 

X e Y limiti dell’intervallo proposti dal giocatore, 
I contatore del numero dei tentativi 
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INIZIO 


OA = 1000 


ESTRARRE N CASUALE 


IN (0,A)I=0 


VISUALIZZAZIONE 

INIZIALE (1) 

LEGGERE XE Y 
SI DIAGNOSTICA 
ea 00 
One» 


I PROVE 


Figura 7.6 


Il flow-chart di Fig.7.6 conduce facilmente al programma di Fig.7.7, come prima è 
possibile migliorare il programma in modo da poter suggerire al giocatore il 
migliore intervallo da provare e nello stesso tempo far giocare il calcolatore contro 
se stesso. 


NOTA: Qui si ha l’interesse di dividere l’intervallo in tre intervalli uguali e di tentare 
l’intervallo centrale. 
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Esempio: per (0, 1000) tentare (333, 666) 


Se la risposta è “Più grande‘ (667, 1000) provare allora l’intervallo (778, 889). 


Questa strategia è la migliore poichè apporta il massimo dell’*informazione*‘, la 
parola informazione è qui utilizzata in senso matematico e non informatico. 


100 REM TROVA IL NUMERO 

110 A = 1000 

115 N = INT CA * RND (1)) 

120I=0 

130 FRINT "TROVARE IL NUMERO COMPRESO TRA 0 E "3A}" FER INQUADRAMENTO" 
140 INFUT X,Y 

145I=I+1 

150 IF X «= Y THEN 180 

160 FRINT "IL FRIMO NUMERO NON DEV'ESSERE FIU’ GRANDE DEL SECONDO" 
170 GOTO 140 

180 IF N « X THEN 210 

190 IF N > Y THEN 220 

200 IF X = Y THEN 230 

205 FRINT "CENTRATO "? GOTO 140 

210 FRINT "FIU’ FICCOLO "? GOTO 140 

220 FRINT "FIU” GRANDE "i? GOTO 140 

230 PRINT "TROVATO DOPO "3I}" TENTATIVI "? FRINT 

240 INFUT "ANCORA UNA FARTITA "R$ 

250 IF R$ = "SI" THEN 115 

260 END 


Figura 7.7 


Esempio di esecuzione: 


TROVARE IL NUMERO COMFRESO TRA 0 E 1000 FER INQUADRAMENTO 
2400,600 

FIU’' FICCOLO 

200,400 

CENTRATO 

2300,350 

PIU’ FICCOLO 

2250,300 

CENTRATO 

2270,300 

CENTRATO 

2290,300 

PIU’ FICCOLO 

2280,290 

CENTRATO 

2285,1 

CENTRATO 

2285, 285 

PIU‘ GRANDE 

9287 ,287 

FIU’ GRANDE 

92878,288 

IL FRIMO NUMERO NON DEV'ESSERE FIU’ GRANDE DEL SECONDO 
?288,288 

PIU’ GRANDE 

?289,289 

TROVATO DOFO 13 TENTATIVI 
ANCORA UNA FARTITA 


Figura 7.8 
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7.3 GIOCO DI MARIENBAD 

Questo gioco molto semplice ha la particolarità di offrire una ‘‘strategia‘* vincen- 
te per il giocatore che muove per secondo. Le regole del gioco sono le seguenti: 
Si dispongono all’inizio del gioco 21 fiammiferi che potranno essere rimossi da 
ciascun giocatore, al proprio turno, nella misura da 1 a 4. Il giocatore costretto a 
prendere l’ultimo fiammifero perde. 

La strategia del giocatore che gioca per secondo è di completare a S il numero dei 


fiammiferi presi dal primo giocatore. Così il primo giocatore si troverà di fronte a 


gruppi composti rispettivamente da 21, 16, 11, 6 e 1 fiammiferi e sarà quindi 


obbligato a prendere l’ultimo. 
3 2 
3 
l 
4 


Problema: Costruire un programma che dia al computer il ruolo del secondo 
giocatore applicando quindi la strategia vincente. 
Il programma dovrà cercare eventuali ‘“imbrogli‘** del primo giocatore. 


Esempio: 


LEGGERE LA PRESA 
DEL PRIMO GIOCATORE 
SIA | QUESTO VALORE 


STAMPARE: “NE PRENDO” 5-1 
“NE RESTANO” N 


Figura 7.9 
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Soluzione: V’algoritmo può essere rappresentato con il flow chart di principio 7.9. 
Per conteggiare gli eventuali “imbrogli” del primo giocatore, si aggiungeranno 
dei test. 


— test che I è intero 
— test che I è compreso tra 1 e 4 


Se un test fornisce una risposta negativa, sarà necessario visualizzare un nuovo 
messaggio di errore e ritornare alla lettura del valore di I per poter reimpostare la 
giocata. 


100 REM GIOCO DI MARIENEAD 

110 REM 

160 FRINT "CI SONO 21 FIAMMIFERI. CHI FRENDE" 
170 FRINT "L'ULTIMO FERDE. SI GIOCA A TURNO " 
175 FRINT "E SI DEVONO FRENDERE 132,3 0 4" 
180 FRINT "FIAMMIFERI ALLA VOLTA." 

185 FRINT iN = 21 

190 INFUT "QUANTI NE FRENDI ? "I 

200 IFIS= INT CI) THEN 400 

210 IFI << 1 0KI > 4 THEN 430 

220 N=N- 5 

230 FRINT "NE PRENDO "35 — I}" NE RESTANO ";N 
240 IF N 1 THEN 190 

250 FRINT "DEVI PRENDERE L'ULTIMO " 

260 FRINT 

270 FRINT "HO VINTO !!"? PRINT 

280 INFUT "ANCORA UNA FARTITA ? "3}R$ 

290 IF R$ = "SI" THEN 185 

300 STOF 

400 FRINT "IL NUMERO DEV'ESSERE INTERO " 

410 GOTO 190 

430 PRINT "NON IMEROGLIARE !!, DEVI FRENDERE 1,2,3 0 4 FIAMMIFERI" 
440. GOTO 190 


Figura 7.10 


CI SONO 21 FIAMMIFERI. CHI FRENDE 
L'ULTIMO FERDE. SI GIOCA A TURNO 
E SI DEVONO FRENDERE 1,2,3 0 4 
FIAMMIFERI ALLA VOLTA. 

QUANTI NE FRENDI ? 1 

NE FRENDO 4 NE RESTANO 16 
QUANTI NE FRENDI ? 1 

NE PRENDO 4 NE RESTANO 11 
QUANTI NE FRENDI ? 2 

NE FRENDO 3 NE RESTANO 6 
QUANTI NE PRENDI ? 4 

NE FRENDO 1 NE RESTANO 1 

DEVI FRENDERE L'ULTIMO 

HO VINTO !! 

ANCORA UNA FARTITA ? SI 

QUANTI NE FRENDI ? 3 

NE FRENDO 2 NE RESTANO 16 
QUANTI NE PRENDI ? 2 

NE PRENDO 3 NE RESTANO 11 
QUANTI NE FRENDI ? 3 

NE FRENDO 2 NE RESTANO 6 
QUANTI NE PRENDI ? 1 

NE PRENDO 4 NE RESTANO 1 

DEVI PRENDERE L'ULTIMO 

HO VINTO !! 

i ANCORA UNA FARTITA ? NO 

Figura 7.11 EREAK IN 300 


Esempio di esecuzione: 
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Il programma di Fig.7.10 scritto partendo dal flow-chart 7.9 comporta i due test 
alle linee 200 e 210, i messaggi di errore alle linee 400 e 430 che poi rinviano alla 
lettura dell’input. 


7.4 CONCLUSIONI 
I tre giochi qui presentati sono molto semplici da programmare per due ragioni: 


— assenza di strategia o strategia molto elementare, 
— nessuna “posizione strategica da valutare‘ come nei giochi tipo dama o 
scacchi. 


Per tutti i giochi che si svolgono su una “scacchiera‘* (OTELLO, DAMA, 
SCACCHI,, si devono manipolare degli indici che aggiungono una complicazione 
supplementare alla difficoltà della strategia. 


Si consiglia ai lettori interessati a questo tipo di problema, di affrontare i giochi 
gradatamente più complessi, per esempio seguendo la seguente progressione: 


Gioco del NIM (simile al gioco di Marienbad ma con più gruppi) 

Gioco del MASTER MIND. 

Gioco di OTELLO implementando prima una strategia semplice poi provare 
a raffinarla progressivamente. 
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CAPITOLO 8 


RICERCA DEL METODO OPERATIVO 


8.0 INTRODUZIONE 


Una parte importante dei problemi di ricerca operativa porta al trattamento di 
grafici (Esempio: Problema del rappresentante, PERT, ricerca topologica etc.). Per 
questa categoria di problemi, la manipolazione degli indici è molto delicata, anche 
per un programmatore esperto. 


Inoltre, gli indici hanno, in generale, dei valori interi e sono molto più adatti per 
questo tipo di problemi gli interpreti BASIC che accettano la rappresentazione 
intera e la rappresentazione in ‘virgola flottante“. 


Il lettore è quindi invitato a non affrontare questa parte fino a quando non abbia 
acquisito una sufficiente pratica con gli esercizi precedenti. 


8.1 RICERCA TOPOLOGICA 


Siano TI, T2, ..., Ty lavori da realizzare rispettando dei vincoli di precedenza: un 
certo numero di coppie (i, j) indicano che il lavoro j non può essere intrapreso fino a 
quando non è stato terminato il lavoro i; un’ ultima coppia (0, 0) indica che la lista è 
terminata. 


Problema: Partendo dai seguenti dati - numero dei lavori e liste di coppie (i, }) - 
trovare un ordine di esecuzione dei differenti lavori soddisfacendo ai vincoli di 
precedenza. 


Per risolvere questo problema si potrà utilizzare il seguente metodo: 
— una tabella T inizializzata a 0 viene modificata dalla lettura delle coppie (I, 
J), quando si legge una costrizione (i, j) che significa che il lavoro J non 
può essere intrapreso fino quando il lavoro i non è terminato, si pone la 
costante 1 in T (I, J), 


— in seguito si cerca un lavoro eseguibile che non abbia costrizioni, o le cui 
costrizioni siano state soddisfatte. Un lavoro K è caratterizzato da: 


T(I, K) = 0 qualunque sia I 
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L’esecuzione di questo lavoro K libera alcuni vincoli. Si esegue dunque: 
T(K, J) = 0 qualunque sia J 


Questo lavoro può essere scelto e il suo numero può essere stampato. Per non 
doverlo selezionare un’altra volta si eporrà: 


T(K, K)= 1 
Si continua cercando un nuovo lavoro fino a che non si è terminato. 
Due casi si possono allora presentare: 


Primo caso: Si contano i lavori selezionati ed è stato possibile stampare gli N lavori. 
Si è dunque trovata la soluzione. 


Secondo caso: Dopo aver fatto tutte le prove possibili, si constata che il numero di 
lavori stampati è inferiore a N. Il problema non ammette quindi soluzione, è 
necessario dunque stampare un messaggio. 


Si applicherà il programma all’esempio dato dallo schema di Fig.8.1 nel quale, 


una freccia che và per esempio dal nodo 8 verso il nodo 5, significa che il lavoro 5 
non può essere intrapreso prima che il lavoro 8 sia stato ultimato. 


4 


Figura 8.1 
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Questo esempio si concretizzerà con i dati della Fig.8.2 che segue: 


1000 DATA 8 

1010 DATA 1,3 
1020 DATA 2,3 
1030 DATA 8,3 
1040 DATA 3,5 
1050 DATA 4,6 
1060 DATA 6,7 
1070 DATA 6,8 
1080 DATA 2,4 
1090 DATA 8,5 
1100 DATA 7,8 
1110 DATA 0,0 

Figura 8.2 


Soluzione: Si scomporrà il problema in tre parti come suggerisce il seguente 
metodo: 


— inizializzazione della tabella a 0, 
— lettura dei dati e preparazione della tabella T, 
— esecuzione dell’algoritmo. 


Nel programma dato in Fig.8.3, queste 3 parti sono l’oggetto di tre sottopro- 
grammi. 


Paragrafo inizializzazione a 0: per alcuni interpreti, tutte le variabili sono, 
all’inizio, inizializzate a 0, ma poichè questo non è prerogativa di tutti i sistemi, è 
preferibile farlo comunque come alle linee 500-540. 


Lettura di dati e preparazione della tabella T e stampa dei dati: durante questa 
lettura, si verifica che i numeri dei lavori sono verosimili e che un vincolo è 
caratterizzato da numeri differenti (altrimenti un lavoro dovrà essere preceduto da 
se stesso). Se si trova un errore, viene stampato un messaggio. Se un vincolo è 
accettato, si pone: T (K, L) = 1. 


Esecuzione di algoritmo (linee 800-960): 


questa parte è più corta di quella precedente, poichè è semplice e vi è una sola 
istruzione di uscita. 


Si noti che per un grafico che comporta un totale di n lavori, vi sono al 


massimo Mach vincoli. Questo appare nell’istruzione FOR della linea 605. 
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100 REM QUESTO FROGRAMMA DI RICERCA TOPOLOGICA DETERMINA IN QUALE ORDINE 
105 REM SI FOSSONO ESEGUIRE DEI LAVORI CHE DEVONO SODDISFARE DELLE 
110 REM COSTRIZIONI DI FRECEDENZA. 

120 FRINT "RICERCA TOFOLOGICA "3? PRINT 

150 DIM T(20,20):N9 = 20 

170 REM INIZIALIZZAZIONE DELLA TABELLA T 

180 GOSUE 500 

190 REM LETTURA, CONTROLLO E STAMPA DEI DATI 

200 GOSUE 600 

210 REM CHIAMATA DELL’ALGORITMO 

220 GOSUE 800 


230 STOF 
500 FOR I = 1 TO N9 
510 FOR J = TO N9 


He 


520 TI, J) 
530 NEXT 4 
540 NEXT I: RETURN 

600 READ N? FRINT "NUMERO DEI LAVORI = "iN? FRINT 
603 PRINT "LISTA DELLE COSTRIZIONI "!? FRINT 

605 FOR I = 1 TO (N — 1) x (N - 1) 

610 READ K,L 
620 IF K = 0 AND L = 0 THEN 720 

630 IFK > L THEN 650 

640 FRINT "ERRORE FOICHE' 2 LAVORI HANNO LO STESSO NUMERO "3kKî STOF 
650. IF K = 0 AND K « N9 THEN 670 

660 FRINT "ERRORE DI NUMERAZIONE DEL LAVORO "3Kî STOF 

670 IF L > 0 ANDL « N9 THEN 690 

680 FRINT "ERRORE DI NUMERAZIONE DEL LAVORO "iL? STOF 

690 TCK,L.) = 18 PRINT Ki TABC 6);L 

700 NEXT I 

710 FRINT "ERRORE NEI DATI "î? STOP 

720 C = I - 1° FRINT 

730 FRINT "NUMERO DELLE COSTRIZIONI = "} 
800 FRINT ? FRINT "ORDINE DEI LAVORI ":I 
810K= 1 

820 FOR J = 1 TO N 

830 IF TC4,}K) = 1 THEN 920 

840 NEXT 4 


C: RETURN 
= 0 


850 I = I + 1? PRINT TAEC 5 x I - 4)}K; 
860 FOR J = 1 TO N 
870 TUK,J) = 0 


880 NEXT YJ 

890 TC(K,K) = 1? GOTO 810 

920 K = K + 13 IFK = N THEN 820 

930 IF I = N THEN RETURN 

940 PRINT "NESSUNA SOLUZIONE FOICHE/ NON SI FUO/ REALIZZARE "3N - I}" LAV 
ORI" 

960 STOF 


Figura 8.3 
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RICERCA TOFOLOGICA 


NUMERO DEI LAVORI = 8 


LISTA DELLE COSTRIZIONI 


w 


NONOS Ss SWUONE- 
DONRSONSUAILDWwW 


NUMERO DELLE COSTRIZIONI = 10 


ORDINE DEI LAVORI 
LR: 4 6 7° ‘8-:3..5 


Figura 8.4 


8.2 IL PERCORSO PIU’ BREVE IN UN GRAFICO 


Il programma qui proposto tratta il caso dell'ordinamento in sequenza di lavori 
di durata conosciuta. 


La numerazione dei lavori è stata stabilita per mezzo di numeri crescenti al fine di 
semplificare la programmazione, ciò permette così di ottenere dei buoni risultati 
anche con un personal computer. 


Tenuto conto della lunghezza del programma, non verrà posto un problema da 
risolvere ma verrà presentata direttamente una soluzione. 
Presentazione dei dati: 


L’insieme dei lavori può essere rappresentato da un grafico orientato avente 
un’entrata ed in principio una sola uscita. Questo grafico orientato non deve 
comportare dei “‘cicli*. Per esempio si potrà avere il grafico della Fig.8.5. 


2 4 


5 


Figura 8.5 
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Nel grafico ogni freccia corrisponde ad un lavoro. Per esempio tra i nodi 2 e 5 
figura un lavoro la cui durata è 4. 


Infatti ogni lavoro è caratterizzato da: 


— un numero di nodo di partenza, 

— un numero di nodo di arrivo, 

— una durata (l’unità può essere una qualsiasi), 
— una designazione. 


Così il grafico della Fig.8.5 corrisponde al listato dei lavori riportato in Fig.8.6. 


2010 DATA 1,2,5, "SOLLEVAMENTO" 
2020 DATA 1,3,9, "SMONTAGGIO" 
2030 DATA 2,3,5, "CAMEIO RUOTE" 
2040 DATA 2,5,4, "REVISIONE" 
2050 DATA 3,4,6, "DISCESA" 

2060 DATA 4,5,1, "SERRAGGIO" 
2070 DATA 0,0,0," 


Figura 8.6 


Per poter leggere questi dati, si utilizzeranno tre tabelle: 


— tabella S che conterrà i nodi di partenza di ogni lavoro 
— tabella F che conterrà i nodi di arrivo di ogni lavoro 
— tabella D che conterrà la durata di ogni lavoro. 


Una stringa di caratteri C$ conterrà le designazioni. Il programma è stato 
limitato al trattamento di problemi con non più di 20 lavori. Per semplificare la 
lettura una variabile, stringa di caratteri, D$ limitata a 20 caratteri riceverà la 


designazione di ogni lavoro. 


Tutte le tabelle saranno quindi dimensionate a 20. La stringa C$ sarà una 
lunghezza di 20 x 20 = 400 caratteri. 


Un sottoprogramma di lettura effettua il seguente lavoro: 


— lettura dei dati, 

— individuazione della fine dei dati caratterizzati da NI = N2 = 0, 

— controllo N1 < N2, 

— inizializzazione di alcune tabelle, 

— conteggio del numero dei lavori. Il risultato è posto nella variabile N. 
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950 

980 

990 

1000 
1010 
1020 
1025 
1030 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 


REM SOTTOFROGRAMMA DI LETTURA E STAMPA DEI DATI 

PRINT "DA A DURATA DESIGNAZIONE" 

FRINT 

FOR I = 1 TO N9 
READ N1,N2,DCI),D$ 
IF Ni = 0 AND N2 = 0 THEN 1100 
IF NI « N2 THEN 1050 
FRINT "NUMERO NON INCROCIANTE"? STOF 

ECN1) = 0?EC(N2) = 0:LCN1) = 0?LCN2) = 0 

SCI) = NI:FCI) = N2!C$(20 x I - 19) = D$ 

FRINT SCI)} TAEC 5);F(I)} TABC 13)}D(I)} TABC 18);D$ 

NEXT I 

N = N9!? GOTO 1110 


N=I- 1 

PRINT 
PRINT "NUMERO DEI LAVORI = "}N 
RETURN 


Figura 8.7 


In questo sottoprogramma dato dal listato di Fig.8.7, le variabili avranno il 
seguente significato: 


N9 


variabile che è stata inizializzata dal valore 20 al lancio del 
programma principale. 


NleN2 variabili che ricevono al momento della lettura rispettivamente il 


numero del nodo di partenza ed il numero del nodo di arrivo, i 
valori di N1 e N2 vengono trasferiti nelle tabelle S e F. 


Tabella E: tempo minimo di inizio di ogni lavoro. 


E è inizializzato a 0. 


Tabella L: tempo massimo di esecuzione di ogni lavoro volendo mantenere 


la tempistica stabilita dal progetto. 
L è inizializzato a 0. 


Il programma principale, dopo aver fatto leggere i dati dal sottoprogramma già 
esaminato, procede con il calcolo del tempo minimo di inizio di ogni lavoro. 


Si prenda l’esempio del grafico in Fig.8.5. Il lavoro che va dal nodo 3 al nodo 4, 
non può cominciare fino a che i lavori che terminano in 3 non siano stati ultimati. 
Questo tempo E(3) è dunque caratterizzato da: 


9 < E(3) 


5+5<E(3) 


} questo dà E(3)=10 


In generale si potrà partire da E(1)=0 e fare il seguente calcolo: 


E(F(1)) = MAX [E(F(1), E(S(A) + DOM] 
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questo si trova alle linee 270-300. 


Per calcolare il tempo massimo di realizzazione di ogni lavoro si partirà in senso 
inverso, per l’ultimo nodo si avrà: 


L(F(N))=E(F(N)) 
L rappresenta la tabella dei tempi massimi di esecuzione, così: 


L(S(1)) = MIN [L(S(1)), L(F(I)) — DM] 


questo calcolo è realizzato alle linee 320 - 360. 
Partendo dal momento in cui si conosce per ogni lavoro: 


— il tempo minimo di inizio E(S(I)), 
— il tempo massimo di esecuzione L(F(I)), 
— la durata normale di esecuzione D(I). 


Si ottiene la durata del ritardo massimo che si può ammettere su questo lavoro 
con: 


FI(I) = L(F(1)) — E(S(1) — D(1I)) 


se F1(I) = 0 significa che ogni ritardo nell’esecuzione di un lavoro ritarderà l’intero 
progetto. La variabile C1 è destinata a contare questi lavori che mostrano tempi 
critici. Questo calcolo è realizzato nelle linee 410-440. 


E’ possibile ora stampare, per ogni lavoro, le seguenti informazioni: 


— numero del nodo d’arrivo 


— numero del nodo di partenza 
dati 
— durata 


— tempo minimo di inizio 
— tempo massimo di esecuzione risultati del calcolo 
— durata del possibile ritardo 


Questo è realizzato dalla sequenza di istruzioni 495-550. 


In seguito vi è il calcolo della durata totale del cammino critico per mezzo della 
variabile C3. 


C:= MAX (C3,L(F(I))) 
linee 570 - 590. 


Poichè il cammino critico comporta dei lavori la cui durata di ritardo ammissibi- 
le è zero, lo si può stampare partendo dal nodo d’inizio (linee 595-720) 


— ricerca del lavoro di partenza (linee 640-660), 
— stampa di questo lavoro, 
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— ricerca del lavoro successivo (linee 700 - 720), stampa e proseguimento 
fino alla fine. 


Il flow-chart di questa parte è mostrato in Fig.8.8. 


RICERCA 
DEL LAVORO 
DI PARTENZA 


LAVORO 


NO 


VISUALIZZAIONE “PIU' DI 
UN PERCORSO CRITICO” 


FINE 


SI 


Figura 8.8 
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100 REM IL FERCORSO FIU’ BREVE CON UN GRAFICO 

110 DIM SC20),F(20),D(20),E(20),L(20),F1(20),C$(400),D$(20) 
115 N9 = 20 

120 REM 

130 REM LETTURA E STAMPA DEI DATI 

150 GOSUE 980 

250 REM INIZ: E CALCOLO DELLA DATA D'INIZIO LAVORI 
260 Cis 0302 = 03C3 = 0 

270 FOR I 1 TON 

2800 Mis EC(5C1)) + DCI) 

2900 IF ECFCOI)) « ss MI THEN EC(FCI)) = MI 

300 NEXT I 


310 REM 

320 LFON)) «= ECFEON)) 

330 FOR IT = N TO 1 STEP - 1 

340 LI GCI)tM2 = LC6FCI)) — DCI) 

9500 IF LIL) è: # M2 OR LCL1) = 0 THEN LILL) = M2 
350 NEXT T 

400 


410 FORI = 1 TON 

4200 F1CT) = LFC0) — ECSCI)) - DUI) 
4230 IF FICI) «= 0) THEN Cis C1 #1 
440 NEXT T 

495 PEXNT 

oo NT "ANAL 
510 PRINT 

520 PRINT "DA —"p"A", "MIN, INIZ.","MAX.,ESEC,","NORM.DUR," 

525 PRINT 

SQ0 FORI = 1 TON 

59350 PRINT SCI)" "3F (I) ESCL), L(FC(I)),F1(I),C$(20 x I - 19) 
550 NEXT T 

540 REM 

570 FORT = 1 T0N 

5900 IF LUPO) è D)3 THEN C3 = LIFCI)) 

590 NEXT I 


SI DEL CAMMINO CRITICO" 


Gislai 

#00 LUNGHEZZA DEL CAMMINO CRITICO "}C3 
6510 

è20 UVA DA". "A" 

520 


FORT = 1 TON 
XF F1CI) = 0 THEN 670 


660 NEXT I 
570 PRINT SCI),FCI) 


680 C2 = C2 + 1 

690 IF I = N THEN 730 

700 FOR ds 1 TO N 

710 IF GC < 3 FCI) OR F1C4) «3 0 THEN 72 

715 Ts J? GOTO 670 

720 NEXT YJ 

7930 IF CL © 8 C2 THEN PRINT "FIU” DI UN CAMMINO CRITICO" 
740 PRI 
2005 
950 FEM O SOTTOFROGRAMMA DI LETTURA E STAMFA DEI DATI 
980 PRINT "DA", "A", "DURATA", "DESIONAZ," 

990 PRINT 

10009 FORI = 1 TO N9 

1010 READ NI,N2,DCI),D$% 

1020 TF NI O AND N2 O THEN 1100 

1025 IF NI N2 THEN 1050 

1030) PRINT "NUMERO NON INCROCIANTE"? STOF 

1050 ENI) = 03FECN2) = OSLCONI) = O:LCONZ) = 0 

1060 ST) s NIIFCOI) = N21C$020 X I - 19) = D$E 

1070 PRINT SCI),FCI),DCUI),D$ 

1080 NEXT I 

1090 N = N92! GOTO 1110 

1100 N=1- 1 


Figura 8.9 
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L’esecuzione porta alla figura seguente costituita da tre parti: 


DA à 
1 2 
1 2 
2 3 
Z o] 
6; 4 
4 Ò 


NUMERO DET LAVORI = 6 


ANALIST DEI. CAMMINO CRITICO 
DA É& MIN.INIZ. 

È. SL Ù) 

1 3 Di) 

23 5 

205 aj 

3 4 10 

4005 16 
LUNGHEZZA DEL CAMMINO CRITICO 
VA DA A 

1 2 

2 3 

3 4 

4 Òò 


BREAK IN 800 


DURATA 


= sO 


MAX.ESEC. 
b) 
10 
10 
17 
16 
17 


17 


Figura 8.10 


DESIGNAZ. 


SOLLEVAMENTO 
SMONTAGGIO 
CAMEIO RUOTE 
REVISIONE 
DISCESA 
SERRAGGIO 


NORM.DUR. 


D) SOLLEVAMENTO 
1 SMONTAGGIO 

0 CAMEIO RUOTE 
8 REVISIONE 

D) DISCESA 

(i) SERRAGGIO 


— stampa dei dati e numero totale dei lavori, 
— analisi del cammino critico, per ogni lavoro: 


— numero del nodo di partenza, 
— numero del nodo di arrivo, 


— data di inizio lavoro, 
— data di fine per non ritardare la realizzazione del progetto 
— durata delle soste, 


— itinerario del cammino critico. 


8.3 IL PROBLEMA DEL COMMESSO VIAGGIATORE 


Un rappresentante deve visitare i suoi clienti disseminati in N città. Si conoscono 
i costi d; per andare dalla città i alla città j. I costi possono essere espressi in distanza 
chilometrica o con altre unità. 


In quale ordine il rappresentante deve visitare i suoi clienti per minimizzare il 
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costo dello spostamento ? In questo problema il rappresentante dovrà, dopo il suo 
giro tornare al punto di partenza. 


Esempio: 
Trento 


Bergamo 


Milano 


Pavia 


Piacenza 


Figura 8.11 


Metodo suggerito: Non si utilizzerà l'algoritmo esatto che è abbastanza complesso e 
lento nell’esecuzione. Si utilizzerà il metodo euristico seguente: 


1. Si sceglie una città di partenza. 

2. Si decide di andare poi nella città più vicina. 

3. Quindi si andrà nell’altra città più vicina non ancora visitata, e così di seguito 
fino a che nonsisiano visitate tutte le città. Si ritorna poi alla città di partenza. Si 
noti il costo del pércorso così determinato e si ricominci scegliendo tutte le città 
come punto di partenza. 


Problema: 


1. Analizzare il problema in vista di una scomposizione in parti relativamente 
corte. 

2. Stabilire un flow-chart di principio molto coinciso facendo apparire dei sotto- 

programmi. 

Stabilire in seguito dei flow-chart di dettaglio per ciascun sottoprogramma. 

4. Scrivere il programma. 


dro 


Si utilizzeranno le seguenti variabili: 


V$ Matrice di stringhe di caratteri che conterrà il nome delle città. 
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D Tabella a due dimensioni che contiene la matrice dei costi 
d; = costo per andare dalla città i alla città j 
di = 0 
T Tabella che contiene una soluzione in corso di esecuzione. 
T2 Tabella che contiene la migliore soluzione trovata durante l’esecuzione del- 
l’algoritmo. 
S Costo della migliore soluzione trovata. 
C Costo della soluzione in corso di studio. 


Soluzione: Questo problema non è difficile a condizione di lavorare con metodo. 


1. Analisi: Il programma totale comporterà differenti parti: 


— lettura dei dati, 

— stampa dei dati, 

— ricerca del migliore itinerario (parte algoritmo), 
— stampa di quell’itinerario. 


Per facilitare la messa a punto della parte relativa all’algoritmo è possibile la 
stampa di itinerari provvisori. E’ quindi facile realizzare la visualizzazione per 
mezzo di un sottoprogramma. Sarà sufficiente quindi inserire un’istruzione GO- 
SUB per questa messa a punto. 


La matrice dei costi può essere simmetrica o asimmetrica. Si esamineranno nella 
lettura i due casi. Così, il lavoro dell’utilizzatore sarà ridotto se c’è una simmetria. 


Si giunge così al flow-chart di principio di Fig.8.12. 


La visualizzazione della matrice dei costi potrà essere la stessa se è simmetrica 
oppure no. Per contro si dovrà tener conto della lunghezza delle linee. 


Per ottenere una corretta stampa si potrà utilizzare una tabella V$ di stringhe di 
caratteri che conterrà il nome delle città. La matrice dei costi sarà rappresentata da 
una tabella quadrata D a due dimensioni. 


Nel flow-chart di Fig.8.12, si è evitato di descrivere l’algoritmo. Progressivamen- 
te entreremo nel dettaglio; occorrerà costruire una descrizione, memorizzarla e 
comparare il suo costo con un’ altra soluzione. Si utilizzeranno le variabili che 
seguono. 


Una tabella T che conterrà la sequenza dei numeri delle città nell’ordine in cui 
appaiono nel circuito. Si avrà dunque: 


C= Y D(T(1, T(1 + 1)) + D(T(N), T(1)) 
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Si potrà ora costruire il flow-chart di Fig.8.13, nel quale la ricerca della città 
successiva non è stata dettagliata. 


INIZIO 


LETTURA E CHIAMARE S/P 1000 PER LETTURA 
VISUALIZZAZIONE 
DEI DATI CHIAMARE S/P 1500 PER VISUALIZZAZIONE 


ESCLUSIONE CHIAMARE S/P 2000 
DELL'ALGORITMO 


VISUALIZZAZIONE 
DELLA CHIAMARE S/P 3000 
MIGLIORE SOLUZIONE PER LA VISUALIZZAZIONE 
TROVATA 


FINE 


Figura 8.12 


Per questo si esamina la situazione in corso di definizione del percorso: si 
supponga che L—I città siano state trovate, i loro numeri sono in T(1)a T(L—I). 
Successivamente si cercherà tra tutte quelle città il cui numero J soddisfa: 


J# T(K) perK = 1,2,..,L—1 


quelle di cui il costo D(T(L—1),J) è minimo al fine di selezionare e immagazzinare 
in T(L) i loro numeri. Questo porta al flow-chart di Fig.8.14 che permette di 
intraprendere la programmazione. 


Il programma è stato scritto per un microcalcolatore che ammette le tabelle di 
stringhe di caratteri. Per un microcalcolatore che non ammette questo tipo di 
tabella, occorrerà apportare delle modifiche simili a quelle che figurano net pro- 
gramma PERT (vedere problema 9.2). 


Per facilitare la comprensione del programma, questo è stato diviso in sottopro- 
grammi, il programma principale provvede alla chiamata di quattro subroutine. 


995 lettura della tabella dei costi 
1500 stampa della tabella dei costi 
2000 algoritmo 
3000 stampa della soluzione trovata 
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RICERCA DELLA 
CITTA' SEGUENTE 
SIA L1 N° CITTA' 
TROVATA, 
CORSO Ci 


MEMORIZZAZIONE 
T(L) = LI 
CUMOLO C = C + Ci 
L= L+1 


creazione di un percorso partendo dal nodo | 


NO 
C=C + D(T(N), T(1)) 
PER RITORNARE AL PUNTO 


naz DI PARTENZA 


MEMORIZZAZIONE 
DI QUESTA 
SOLUZIONE 

SE È LA MIGLIORE 


NO 


RETURN 


Figura 8.13 
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SI 
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D(T(L - 1) J=>C1 


INIZIO 


K=1 | LACITTA JÈ GIA 
STATA SCELTA 


C1= D(T(L — 1), J) 
L1=J 


NO 


RETURN 


Figura 8.14 


La subroutine ALGORITMO fa, anche lei, appello a due altri sottoprogrammi: 


2500 ricerca della successiva città 
2700 test se la soluzione trovata è migliore e se si immagazzinare questa 
soluzione. 


20 REM V$ RAFFRESENTA IL NOME DELLE CITTA” 
25 REM T = TABELLA DI LAVORO I CUI ELEMENTI CONTENGONO 


30 REM I NUMERI DELLE CITTA‘ DEL FERCORSO 

35 REM T1 CONTIENE I NUMERI DELLE CITTA” 

40 REM DELLA MIGLIORE SOLUZIONE TROVATA 

45 REM = E’ LA MATRICE DELLE DISTANZE 0 DEI COSTI 


90 DIM V$(20),T020),T1(20),D(20,20) 

100 FRINT "FROGRAMMA DEL COMMESSO VIAGGIATORE " 
110 PRINT 

120 READ S$ 

125 IF S$ = "SYM" THEN GOSUE 995 

126 IF S$ < = "SYM" THEN GOSUE 800 

130 GOSUE 1500 

140 GOSUE 2000 

150 GOSUE 3000 

780 STOF 

790 REM LETTURA DELLA TABELLA DEI COSTI QUANDO LA MARICE NON E” SIMMETRICA 


800 READ N 

810 FOR I = 1 TON 

820 READ V$C(I) 

825 NEXT I 

827 FORI = 1 TON 

830 FOR J = 1 TON 
840 READ DUI, J) 

850 NEXT 4 

860 NEXT I 

870 RETURN 

990 REM LETTURA DELLA TABELLA DEI COSTI DI MATRICE SIMMETRICA 
995 READ N 

1000 FORI = 1 TON 
1010 READ V$C(I) 

1020 NEXT I 

1030 FORI = 1 TON 
1040 DCI,I) = 0 

1050 FORJ =I+ 1 TON 
1060 READ DCI,J) 
1070 D(J,I) = DI, J) 
1080 NEXT 4 

1090 NEXT I 

1100 RETURN 


1480 REM 

1490 REM SUEROUTINE FER STAMFARE LA MATRICE DEI COSTI 
1500 FRINT "TABELLA DEI COSTI" 
1510 FRINT 

1520 FOR I = 1 TON 

1530 FRINTO TAEC 5 * T)}V$(I)} 
1540 NEXT I 

1550 FRINT 

1555 PRINT 

1560 FOR I = 1 TO N 

1570 FRINT V$(I); 

1580 FOR 4 = 1 TO N 

1590 PRINT TAEC 5 * J);D(I,J); 
1600 NEXT J 

1605 FRINT 

1608 FRINT 

1610 NEXT I 

1620 RETURN 


12700REN Figura 8.15 
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1970 
1980 
1990 
2000 
2002 
2005 
2010 
2020 
2030 
2040 
2050 
2060 
2065 
2070 
2090 
2110 
2470 
2480 
2490 
2500 
2510 
2515 
2520 


2525 


2530 
2540 
2550 
2560 
2570 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
3000 
3010 
3015 
3020 
3030 
3040 
3050 
3055 
3060 
3070 
3080 
5000 


REM 
REM INIZIO DELL’ALGORITMO 
REM 
S = 1E38 
FORI =1 TON 
C=0 
TI) = I 


FOR L = 2 TO N 
GOSUE 2500 


T(L) = L1 
C= C+ ci , A 0a 
NEXT L aggiungere il costo di ritorno 
C = C+ DAUTCON),TC1)) <,,% F 

GOSUE 2700 della città di partenza 

NEXT I 

RETURN 

REM 

REM RICERCA DELLA CITTA’ SUCCESSIVA 

REM 

Ci = 1E38 


FOR J = 1 TON 
FOR K = 1 TOL - 1 


la città seguente non deve essere 
IF TC0K) = 4 THEN 2560 


NEXT K già stata scelta 

IF D(TC(L — 1),4) > = C1 THEN 2560 ; n NA 
C1 = DITE — 1),4) } si considera la più vicina 
li = 4 

NEXT 4 

RETURN 

REM 


REM LA SOLUZIONE TROVATA E’ LA MIGLIORE ? 

REM SE SI, TRASFERIRE T IN T1 E C INS 

IF S< = C THEN 2750 

S=C 

FOR K = 1 TON 

T1C(K) = TCK) 

NEXT K 
RETURN 
FRINT 
PRINT TABC 9); "SOLUZIONE" 
PRINT 
FOR L = 1 TON - 1 

PRINT V$(T1(L));" AU "3V$(T1C(L + 1)),D(T1(L),T1(L + 1)) 
PRINT 

NEXT L 

PRINT V$(T1(N));" AU "3V$(T1C(1)),D(T1C(N),T1C1)) 
PRINT 

PRINT "COSTO TOTALE î? "} TABC 32);5S 

RETURN 

END 


Figura 8.16 


PROGRAMMA DEL COMMESSO VIAGGIATORE 


TABELLA DEI COSTI 


VR ES FU TR 
VE D) 45 67 13 
ES 47 o 29 37 
FU 68 30 Di) 73 
TR 19 26 74 0 
CR 20 24 22 43 
EG 57 293 25 60 
FO n9 40 13 98 
MI 81 36 37 75 

SOILUZIONE 

TRA UI 12 
UR A CR 40 
CECA Pu 22 
PU A FI 12 
FCC A EG 36 
EG A MI 13 
MI A Be 36 
ES. A TR 37 
COSTO TOTALE è 209 


Dati corrispondenti: 


155 
160 
170 
180 
190 
200 
210 
22 


230 


DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 


"SYM",8 


e 
ts] 


46 


REM LETTURA DELLA TABELLA DEI COSTI QUANDO LA MARICE NON E’ SIMMETRI 


45,67,13,40,68,89,81 
29,37 ,22,23,41,36 
73,21,24,12,37 
42,60,95,73 
36,33,49 

36,13 

47 


Figura 8.17 
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FROGRAMMA DEL COMMESSO VIAGGIATORE 


TABELLA DEI COSTI 


VR 


FC 


VR ES FU TR CR 
D) 5 67 13 40 
45 0 29 37 22 
67 29 Di) 73 21 
13 37 73 D) 42 
40 22 21 42 0 
68 23 24 60 26 
89 41 12 95 33 
81 26 a7 73 49 


SOLUZIONE 


VR 13 
CR 40 
FU 21 
FC 12 
EG 36 
MI 13 
ES 36 
TR 37 


O TOTALE ! 208 


Dati corrispondenti: 


155 


160 
170 
180 
190 
200 
210 
220 
230 
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DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 


"SYM",8 
VR,ES,FU,TR,CR,EG,FC,MI 
45,67, 13,40,68,89,81 
29,37 ,22,23,41,36 
73,21,24,12,37 
42,60,95,73 

36,33,49 

36,13 

47 


Figura 8.18 


47 


Percorso corrispondente al risultato precedente: 


Trento 


Bergamo 


Cremona 
Piacenza 


Figura 8.19 


Questo risultato non corrisponde al circuito più corto che può essere trovato e il 
cui costo è di 206. Questo percorso di costo minimo è il seguente: 


Trento 


Brescia 


Piacenza 


Figura 8.20 


NOTE: 1. Se una città è ‘‘equidistante‘ da due altre città e se i costi sono minimi, 
l’algoritmo programmato non esegue test successivi con ciascuna delle città, ma 
sceglie sistematicamente la prima in ordine di indice. 


Per obbligare l’algoritmo a effettuare il secondo test, occorre aumentare artifi- 
cialmente di poco il costo del primo tragitto e domandare una seconda esecuzione. 
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2. Per aumentare la velocità dell’algoritmo è possibile modificarlo, imponendogli 
dei tratti che fanno parte del percorso di costo minimo. 


8.4 CONCLUSIONI 


Si sono visti tre semplici programmi di ricerca operativa. Il lettore ha potuto 
rendersi conto della lunghezza di questi programmi nonostante la semplicità 
apparente del problema. In generale ogni volta che è necessario ‘‘percorrere un 
grafico‘ è necessario ricorrere a manipolazioni piuttosto delicate di indici che 
generano difficoltà di programmazione. 


Si consiglia al lettore interessato a questo problema di documentarsi sui seguenti 
argomenti: 


— l’algoritmo di Kruskal, 

— il problema dell’ottimizzazione del trasporto (algoritmo dello 
“‘steppingstone‘‘), 

— l’ottimizzazione del flusso in un grafico (algoritmo di Ford-Fulkerson), 

— programmazione lineare (algoritmo del simflesso). 
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CAPITOLO 9 


STATISTICA 


9.0 INTRODUZIONE 


Per l’elaborazione di diverse statistiche si fa spesso riferimento a microcalcolato- 
ri le cui qualità (velocità di calcolo, accesso ad un gran numero di dati) li rendono 
particolarmente adatti a questo tipo di applicazioni. 


In questo capitolo si studieranno dei programmi di statistica semplici, ma molto 
utili. Per esempio per il calcolo del tasso di crescita studiato al Capitolo 6, si fa 
riferimento a un sottoprogramma di regressione lineare spiegato in questo capitolo. 


Il numero di esercizi presentati è limitato al fine di non annoiare il lettore che non 
sia interessato a questo tipo di programmi. E° necessario sapere, tuttavia, che esiste 
un gran numero di programmi per questo campo. 


9.1 MEDIA DI UNA SERIE DI MISURE 


Calcolare la media aritmetica M di una serie di misure. In questo esercizio si 
suppone che tutti i dati siano incorporati nel programma. Un valore numerico 
uguale a —999 indica la fine dei dati. 


Problema: Analizzare il problema, disegnare il flow-chart e scrivere quindi il 
programma. 


Soluzione: Ogni misura è utilizzata una sola volta per il calcolo nel corso della 
somma. Non è dunque necessario utilizzare una tabella. Il numero di misure sarà 
registrato in una variabile N al fine di poter effettuare la divisione. 


Alfine di non sommare il valore —999, è necessario effettuare un test A = —999 
se A# —999 effettuare la somma:M = M+A 
N=N+1 


se A = —999 tutti i dati sono stati posizionati ed è necessario effettuare la 
divisione M = M/N. 
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Questo condurrà al flow-chart di Fig.9.1. 


INIZIO 


IL VALORE 
DIM 


FINE 


Figura 9.1 


La programmazione di questo flow-chart è facile. La sola difficoltà risiede nel 


curare la presentazione dei risultati. 


100 M = 0:N= 0 
110 READ A 
120 IF A= - 999 THEN 170 
130 N=N+ 1 
140 M=M+ A 
150 GOTO 110 
170 M=M/ N 
180 FRINT "NUMERO DELLE MISURE= "iN 
190 FRINT 
200 FRINT "MEDIA "} TABC 20);"= "iM 
210 DATA 12,25,15,0,-999 
220 END 
Figura 9.2 
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9.2 MEDIA, VARIANZA, DEVIAZIONE STANDARD 


Per stimare la media, la varianza e la deviazione standard di una serie di N misure 
si possono utilizzare le seguenti formule: 


; 1 
Media a 
v x A(I) 
DI N ni 
Varianza V= Reî D (Aq) —- MY? 


Deviazione standard E = VV 


Come precedentemente supposto i dati sono incorporati nel programma e il 
valore numerico —999 indica la fine dei dati (ricerca della fine dei dati). 


9.2.1 Primo problema: Si tratta di partire da una serie di misure supposte incorpora- 
te nel programma per stimarne la media, la varianza e la deviazione standard per 
mezzo delle formule precedenti. Per questo viene proposto di trattare questo 
problema in tre parti. 


Parte a: Disegnare il flow-chart che permette di ottenere il valore numerico 
delle tre stime. 


Parte b: Modificare la formula che fornisce il valore di V affinchè il flow- 
chart abbia un solo loop di calcolo. 


Parte c: Scrivere un programma corrispondente al secondo flow-chart. 
Soluzione: 
Parte a: Si è condotti a costruire il flow-chart di Fig.9.3 che consta di due parti: 
— una prima parte destinata a calcolare la media, 
— una seconda parte che permette di ottenere la varianza e la deviazione 


standard. 


Questo dà il flow-chart di Fig.9.3 che comporta due loop di calcolo e due letture 
dei dati. 
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Nel caso in cui i dati siano poco numerosi e incorporati nel programma, questa 
doppia lettura di dati non costituisce una particolare difficoltà. Per contro nella 
pratica, quando si tratta di archivi di dati di una certa consistenza, una doppia 


INIZIO 


Figura 9.3 
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lettura raddoppia i tempi di calcolo se si è in ‘*multiprogrammazione‘‘* e aumenta 
notevolmente il ritardo di restituzione se il microcalcolatore utilizzato lavora in 
“multiprogrammazione‘. Per questa ragione si prova a minimizzare il numero 
d’accessi all’archivio. 


Parte b: Sviluppando la formula che fornisce il valore di V si ottiene: 


DE 0a 1 
Veg Db A(i) - 2M Y AG)+nM?| 


i=1 


con 


1 > er 
NSSTZi b A(i) - NM | 


Questa formula permette di calcolare M e V per mezzo di un solo loop di calcolo 
come si vede nel flow-chart che segue. 


LEGGERE A 


N=N+1 
M=M+A 
A2=A2+AxA 


Figura.9.4 
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Parte c: Questo flow-chart è molto semplice e la programmazione è facile. La sola 
difficoltà consiste nel curare la presentazione dei risultati (Figg.9.5, 9.6). 


100M= 0 

110 N = 0 

120 A2 = 0 

130 READ A 

140 IF ÀA= - 999 THEN 190 

ISO N=N+ 1 

160 M=M+ A 

170 A2 = AZ+AX A 

180 GOTO 130 

190 M=M/ZN 

200 V = (AZ - N X Mx.M) / (N -— 1) 
210 E = SQR (4) 

220 FRINT "NUMERO DELLE MISURE= "3N 
230 PRINT " MEDIA= "3M 
240 FRINT " VARIAZIONE= "}V 
250 FRINT "DEVIAZIONE STANDARD= "3E 
260. STOF 


300 DATA 9,9,9,10,8.5,9;10.1 
310 DATA 10,9,8,10,2 

320 DATA -999 

330 END 


Figura 9.5 


Esempio di esecuzione: 


NUMERO DELLE MISURE= 9 
MEDIA= 9.,61111112 
VARIAZIONE= .373611014 
DEVIAZIONE STANDARD= .611237281 


Figura 9.6 


9.2.2 Calcolo dei coefficienti d’asimetria e appiattimento 


Modificare il programma precedente per calcolare il valore numerico delle stime 
seguenti: 


: = Des =; 
Coefficiente di asimmetria NE Y (A) - MY 
Isl 


Y (AQ) - MY 


i i iattimento ——; 
Coefficiente di appiattime NE 


Soluzione: Si noti che il momento d’ordine 2 si scrive: 
] da . sr vr) 
M, = N di (A(i) — MY 
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Corrisponde a una stima per vie traverse della varianza. 


Si pone , 
Vis > (A(1) — M)° 


Si avrà 
let 
CON-1°' 
1 
M.=% i 


Si possono ora scomporre le due formule che danno S e P 


| = 1 = 
S= Gi [® A(i)' - 3M X A(i) +3 M° £ A(i) — NM'] 
N(R) 
| da si 
i 
N) 
! gra = = 
P= — 7 [E Ali)" — 4M ® ACI) +6 MÈ X Ai) — 4M'Z ACI) + NM'] 
N(S) 


N nà Re S: 
= 73 [D Aci) - 4M ® Ai)' + 6M° © ACi)È — 3 NMI] 


E' sufficiente ora inserire nel loop i calcoli di X A(i) e di Z A(i)‘ che saranno posti 
nelle variabili A3 e A4, per ottenere in seguito i coefficienti cercati. 


1 > 
S= [A3-3M A2+2NM7] 


so) 


P= 33: [A4-4M A3+6M°A2-3NM'] 


La costruzione del programma di Fig.9.7 non presenta a questo punto più alcuna 
difficoltà. 
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100 AL = 0 

110 N = 0 

120 A2 = 0 

125 A3 = 0 

127 A4.= 0 

130 READ A 

140 IF A= - 999 THEN 190 

150 N=Nt+ 1 

155 Al = Al + A 

160 X= Ax A 

162 A2 = A2 + X 

165 A3 = AB + Xx A 

167 A4 = AA + XxX 

180 GOTO 130 

190 M = AL /N 

200 V = CA2 - N * Mx M) / (N - 1) 
210 E = SQR (UV) 

220 FRINT "NUMERO DELLE MISURE= "} 
230 PRINT" MEDIA= "};M 
240 PRINT" VARIAZIONE= "}V 


250 FRINT "DEVIAZIONE STANDARD= "E 

253 M2 = MX M 

255 S = (A3 — 3 X Mx AZ + 2 Xx M2 x A1) / (N x V x E) 
F_= CA4 - 4 Xx Mx A3 + 6 X MZ x AZ — 3 X N X M2 x M2) / (N x V x V) 

270 FRINT " COEFFICIENTE DI SIMMETRIA= "85 

280 FRINT "COEFFICIENTE DI AFFIATTIMENTO= "}F 

285 STOF 

300 DATA 2,2.5,3,3,5,4 

310 DATA -999 

320 DATA -999 

330 END 


Figura 9.7 


Esempio di esecuzione: 


NUMERO DELLE MISURE= 5 
MEDIA= 3 
VARIAZIONE= ,625 
DEVIAZIONE STANDARD= .790569415 
COEFFICIENTE DI SIMMETRIA= 0 
COEFFICIENTE DI AFFIATTIMENTO= 1,088 


Figura 9.8 


Altra serie di dati ad esecuzione corrispondente 


NUMERO DFI LE MISURE= 5 
MEDIA= 3 
VARIAZIONE= 2,5 
DEVIAZIONE STANDARD= 1.58113883 
COEFFICIENTE DI SIMMETRIA= 0 
COEFFICIENTE DI APFIATTIMENTO= 1,088 
BREAK IN 285 


Figura 9.8 
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NOTE: 


1. Questi ultimi due sistemi di stima non sono valevoli con tutte le leggi. Devono 
quindi essere utilizzate con prudenza. 

2. Coefficiente di simmetria vale 0 se il campione è simmetrico. 

3. Coefficiente di appiattimento ingrandito nel valore assoluto con l’appiattimento 
della funzione di densità. 


9.3 REGRESSIONE LINEARE 


Cercare una retta che approssima nel migliore dei modi possibile l’insieme dei 
punti (x;y) sperimentali. Per questo si utilizza generalmente il criterio dei ‘minimi 
quadrati‘ che consiste nel determinare i coefficienti a, © in modo che: 


Y (ax, + b — y,) 


im 
Per minimizzare questa quantità è necessario calcolare a e b come 
nZxy=Ex)Ey) 
nXx-(Xx)Y 


nZx-A%Zx, 


IENA 


Per assicurarsi della validità ‘‘statistica‘* del calcolo, si può calcolare il coefficien- 
te r dato da 
pot Y 
x(Y,-Y) 
Se r è vicino a l il calcolo è valido dal punto di vista statistico; in caso contrario, la 


regressione lineare si applica piuttosto male alla serie dei dati. Si può allora 
calcolare una varianza di a e d e determinare l’intervallo limite di a e db. 


r = segno di db 


Esercizio: Scrivere un sottoprogramma che partendo dalle tabelle T(100) e Y(100) 
determini una retta di regressione e calcoli il coefficiente r. 


Il calcolo dei coefficienti A , B costituisce un sottoprogramma che inizia in 1000, 
e il calcolo del coefficiente R dev'essere realizzato in un sottoprogramma che inizia 
in 600. 


Soluzione: La parte calcolo di A e B deriva dalle formule precedenti con l’ausilio di 
un loop di calcolo, si ottiene: 


x, X Vox, 2 x Vis 


poi partendo da queste quantità determinare i valori di A e B di cui al flow-chart di 
Fig.9.9 
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Ul=U2=W=0(1) 


UI=UI+T(1) 
VI=VI+Y(l) 
U2=U2+T(1)*T(1) 
W=W+T(1)*Y(1) 
1=1+1 


UlxVI 


N 
dI 
Ul: dia 

N 


MIAXUI 


B= N 


Figura 9.9 


Per il calcolo di R si effettuerà un altro LOOP di calcolo 


Ul=U2=0 


U1=U1+(Y(1)-AxT()-8)? 


U2=U2+(A-T(1)+ BIO 


I=I1+1 


R=SIGN (B) WII, 


FINE i 
Figura 9.10 
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Partendo da questi due flow-chart, scrivere il programma con l’esempio di 
esecuzione di Fig.9.11 è cosa facile. 


DIM TC(100),YC100) 


READ N 
FORI = 1 TON 
READ TCL) YI) 
NEXT I 
GOSUE 1000 
FRIN' TAEC 9); "INCLINAZIONE = ";A 
NORDINATA ALL'ORIGINE = "3E 


ONT", "MISURA Y". "CALCOLO Y" 


LTON 
3 AMO OTCI) + E 
PRINT TCI).,Y(I),Yil 
NEXT I 
STOF 
DATA 5 
DATA 0,1,1,1,5,2,2,4,3,6,4 


1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 


+ YI) 

+ TCI) * TCI) 

1090 : + YCI) X YCI) 

1100 MO + TCI) * YCI) 

1110 NEXT I 

1120 Asso (W — UL X V1 / N) / (U2 — U1 * Ul / N) 
1130 E = (V1 - Ax Ul) / N 

1140 RETURN 

1200 END 


INCLINAZIONE = ,5 
ORDINATA ALL'ORIGINE s: 1 


Ti MISURA Y CALCOLO Y 
0 1 1 
1 1,5 1,5 
2 2 2 
4 3 3 
6 4 4 
Figura 9.11 


Regressione lineare senza calcolo del coefficiente r. 
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DIM TC(100),Y(100) 
READ N 
FORI = 1 TON 

READ T(I),YCI) 

NEXT I 

GOSUE 1000 

GOSUE 600 

PRINT TAEC 9); "INCLINAZIONE = "; 
PRINT "ORDINATA ALL'ORIGINE = "E 
PRINT TAEC 7)}"COEFFICIENTE Rs ";R 
PRINT 

PRINT "T", "MISURA Y", "CALCOLO Y" 
PRINT 

FORI = 1 TO N 

YI Gs A XTC) + E 

PRINT TOD,Y(I),Y1 


NEXT I 
STOP 
DATA S 
DATA 0,1,1,1,5,2,2,4,3,6,4 
Ul se 0 
= 0 


1 TON 
UL + (CYCT) — fi * TOI) — Ei) 42 
UZ + CA X TAI) + E VI ZON) 4 2 
L 
SEN (E) * SQR (1 U1 7 U2) 


1000 Ul = 0 


1010 0) 

1020 n 

1040 = 0 

1050 \I=1T0N 

1060 Ul + TAI) 

1070 = V1 + YCI) 

1080 = U2 + TCI) * TCI) 
1100 Wo+ TL) x YeIo 


1110 NEXT I 

1120 A Gs (4 > UL XxX VI Z N) 7 (42 — Ul * Ul / N) 
1130 Eos (VI — A GX UL) / N 

1140 RETURN 

1200 END 


Figura 9.12 


Il programma della figura 9.12 corrisponde al calcolo dei coefficienti a e d e del 
coefficiente r. 


Esempi di esecuzione con diversi set di dati. 


INCLINAZIONE © +5 

ORDINATA ALLORIGINE 1 

COEFFICIENTE KR 1 
Ti MISURA Y CALCOLO Y 
Ù) 1 1 
1 1,5 1,5 
2 è 2 
4 3 3 
é 4 4 
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INCLINAZIONE # 503125 
ORDINATA ALL'ORIGINE = 1.00312 
COEFFICIENTE R = 992166 


T MISURA Y CALCOLO Y 

0 .95 1.00312 

1 1.55 1.50625 

2 2.05 2,0937 

4 2.95 3.01562 

4 3.05 3.01562 
INCLINAZIONE © .480603 


ORDINATA ALL'ORIGINE 1.25043 
COEFFICIENTE R = ,932274 


T MISURA Y CALCOLO Y 
(0) +95 1.225043 
1 1.55 1.73103 
2 2.95 2.21164 
4 3.05 3,17224 
6 4 4. 13405 


Questo esempio mostra la sensibilità del coefficiente di correlazione utilizzato. 


9.4 CONCLUSIONE 


Questi esercizi mostrano che la programmazione di calcoli elementari del tipo 
media, varianza, etc. non presenta particolari difficoltà. L’esercizio di calcolo di 
una regressione lineare è molto utile poichè viene utilizzato nel Capitolo 6 per il 
calcolo del tasso di crescita. 


Invece quando è necessario addentrarsi in calcoli più sofisticati (prove statist 
che, regressione multipla, regressione polinomiale, etc.), i programmi si allungano 
diventano sensibili agli errori di arrotondamento. 


Il lettore troverà in libri specializzati degli esercizi supplementari: 


— istogrammi, 
— regressione lineare multipla. 
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CAPITOLO 10 


VARI 


10.0 INTRODUZIONE 


Questo capitolo contiene un insieme di esercizi che non fanno parte di nessun 
argomento trattato nei capitoli precedenti, ma che presentano un interesse sul 
piano informatico sia perchè fanno appello a delle astuzie di programmazione, sia 
perchè l’elaborazione del flow-chart non è evidente. 


10.1 I SEGNI ZODIACALI 
Partendo dal giorno e dal mese di nascita, determinare il ‘segno zodiacale‘ al 


quale si appartiene. Questo è facilmente realizzabile in quanto è sufficiente consul- 
tare la tabella 10.1 per ottenere la risposta. 


SEGNO PERIODO 


CAPRICORNO 23 dic. al 19 gen. 
AQUARIO 20 gen. al 19 feb. 
PESCI 20 feb. al 20 mar. 
ARIETE 21 mar. al 19 apr. 
TORO 20 apr. al 20 mag. 
GEMELLI 21 mag. al 20 giu. 
CANCRO 21 giu. al 21 lug. 
LEONE 22 lug. al 22 ago. 
VERGINE 23 ago. al 23 set. 
BILANCIA 23 set. al 22 ott. 
SCORPIONE 23 ott. al 21 nov. 


SAGITTARIO 22 nov. al 22 dic. 


Figura 10.1 
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Esercizio 1: 


Scrivere un programma che partendo dal giorno e mese di nascita determini il 
segno zodiacale corrispondente. 


Si supporrà di disporre di un BASIC che accetti le tabelle costituite da stringhe di 
caratteri. 


Esercizio 2: 


Come il precedente ma disponendo questa volta di un BASIC che non accetta le 
tabelle costituite da stringhe di caratteri. 


Soluzione esercizio 1: 


Si devono confrontare il giorno del mese J con un limite L che dipende dal mese 
(L da 20 a 23). 


— Se J< L l’indice da utilizzare è I = M. 

— Se J2L l’indice da utilizzare è I= M + 1 salvo che questo diaI= 13in tal 
caso è necessario attribuire a I il valore 1 (caso di una persona nata tra il 23 
e il 31 dicembre). 


Perchè L assuma il valore corretto si dovrà prima attribuirgli il valore 20 poi 
utilizzare un’istruzione ON...GOTO per effettuare il salto a una serie di istruzioni L 
=L+1 


In questo modo si eviteranno parecchie prove e istruzioni GOTO. 
Questo metodo può essere rappresentato dal flow-chart di Fig.10.2 che si presenta 
più difficile del programma di Fig.10.3 corrispondente. 


Soluzione esercizio 2: 


Il programma precedente determina un indice J che indica il numero del segno, 
ma questo indice non può essere utilizzato in questa forma. 

I segni zodiacali più lunghi da scrivere sono CAPRICORNO e SAGITTARIO 
con 10 lettere. Si potrà quindi utilizzare una variabile, stringa di caratteri, A$ con 
lunghezza 120 che potrà contenere i dodici segni di 10 caratteri (quelli di lunghezza 
inferiore saranno completati con l’inserimento di spazi). 


Partendo dall’indice I ottenuto precedentemente si scriverà: 


MID$ (A$, 10 X I — 9,10) 
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Occorre ora poter inizializzare correttamente A$ per questo è possibile utilizzare 
diversi metodi, per esempio: 


115 DIM A$(120) 
120 A$="CAFRICORNOACQUARIO miFESCI 


ann ARIETE aan 1 TORO nana" 
125 AS=A$+ "GEMELLI anna CANCRO4 a a ALEONE ann n VERGINE anna E TLANCI Anna" 
150 A$=A$+"SCORFIONEMmSAGITTARIO" 


Ogni segno completandolo con un numero di caratteri ‘“spazio‘, inizia nella 
giusta posizione (1,11,21,31 etc.) e potrà essere facilmente selezionato. 


INIZIO 


LEGGERE 
GIORNO, MESE 


STAMPARE 
AS(I) 


FINE Figura 10.2 


177 


115 DIM A$(C12) 

120 FOR I = 1 TO 12 

125 READ A$CI) 

130 NEXT I 

140 PRINT "GIORNO, MESE DI NASCITA "} 
145 INFUT J,M 


180I=M 

190 L = 20 

200 ON M GOTO 600,600,500,600,500,500,400,300,300,300,400,300 
300L=L4+ 1 

400L=L4+ 1 

500L=L4+1 

600 IF J € L THEN 610 
605I=14+1 

610 IF I = 12 THEN 620 
615 I= 1 

620 FRINT "VOI SIETE "3A$(I)} 
630 PRINT 


650 GOTO 140 

900 DATA "CAFRICORNO","AQUARIO","FESCI", "ARIETE" 

910 DATA "TORO", "GEMELLI", "CANCRO", "LEONE", "VERGINE" 
920 DATA "EILANCIA","SCORFIONE", "SAGITTARIO" 

990 END 


Figura 10.3 


GIORNO, MESE DI NASCITA 212,10 
VOI SIETE EILANCIA 


GIORNO, MESE DI NASCITA ?19,1 
VOI SIETE CAFRICORNO 


GIORNO, MESE DI NASCITA ?20,1 
VOI SIETE AQUARIO 


GIORNO, MESE DI NASCITA ?21,5 
VOI SIETE GEMELLI 


GIORNO, MESE DI NASCITA 223,9 
VOI SIETE EILANCIA 


GIORNO, MESE DI NASCITA ? 
PREENTER 
Ls 


Figura 10.4 


Attenzione questo programma non si ferma da solo! 


10.2 REGINE SU UNA SCACCHIERA 
Un problema classico per il gioco degli scacchi è il seguente: trovare tutte le 


soluzioni che permettono di piazzare 8 regine su una scacchiera senza che vi siano 
due regine sulla stessa linea orizzontale, verticale e diagonale. 


Problema: 


Generalizzando una scacchiera N X N, impostare tutte le soluzioni che permetto- 
no di posizionare N regine senza che ci siano due regine ‘in presa“ tra loro. Si farà 
variare N da 2 a 8. 


Si elimineranno le soluzioni che si deducono per simmetria. 
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Metodo proposto: 


Figura 10.5 


Una tabella D conterrà le regine posizionate nella soluzione in corso di elabora- 
zione. 


Per esempio la soluzione proposta in Fig.10.4 corrisponde a: 


D(1)=1 D(5) = 3 
D(2) = 5 D(6) = 7 
D(3) = 8 D(7)=2 
D(4) = 6 D(8) = 4 


Con il modo di rappresentazione dei dati scelti, affinchè due regine non siano “in 
presa‘ occorre: 


— che non siano sulla stessa colonna, e qui è realizzato poichè c’è solo una 
regina per colonna, 


— che non siano sulla stessa linea, ciò significa che è necessario porre: 
D(I) # D(J) 
— che non siano sulla stessa diagonale. Per questo è necessario porre: 


D(J) — D(I) # JHI diagonale a 45 gradi 
D(J) — D(I) # I—J diagonale a —45 gradi 


179 


Questi due test possono essere riuniti nel seguente test: 


ABS( D(I) — D(J))#I—J 


Criteri di ricerca delle soluzioni: 
Si parte da D(1)=1. 


Si cerca la prima posizione ammissibile per D(2), cioè una posizione tale per cui 
D(2) non sia ‘in presa‘ con D(1). 


Si cerca allora una posizione ammissibile per D(3) e si continua fino a trovare una 
posizione per D(N). 


Se non si trova la posizione ammissibile per D(I), si proverà a spostare D(I-1) 
soddisfacendo alle specifiche: D(I—1) non dovrà essere “in presa‘ con le regine 
precedenti (D(1),D(2)...D(I—2)) poi si cercherà una nuova posizione ammissibile 
per D(I). 


Più precisamente l’algoritmo proposto è il.seguente: 
Per I che varia dala N 


1. Porre D(I1)= I 
2. Verificare che la nuova regina D(I) nonsia ‘in presa‘‘ con una delle regine 
già posizionate. 
Se è “in presa‘ andare in 3. 
Altrimenti, se I = N si ha una soluzione che si può stampare procedendo 
oltre (andare in 3). 
Se I< N porre I=I+ 1 e andare in I. 
3. Cercare una nuova posizione per D(I): 
Porre D(I)= D(1) +1 
Se D(I) < N andare in 2 
Se D(I) > N non ci sono posizioni convenienti: porre I=I—1, andare in 3. 


Per eliminare le soluzioni che si deducono per simmetria da una soluzione già 
trovata: 


a) Si farà variare D(1) da 1 a N72 soltanto (si eliminano così tutte le soluzioni 
simmetriche rispetto all’asse orizzontale). 


b) Non si stamperà alcuna soluzione D(1) > D(N) poichè questa soluzione corri- 
sponde a una soluzione simmetrica rispetto all’asse verticale già stampata. 
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Flow-chart 


Il sottoprogramma di Fig.10.6 rappresenta l’algoritmo descritto precedentemen- 
te. 
Il flow-chart di Fig.10.7 rappresenta il programma principale. 


RETURN 


tml-1 0(1)=D(1)+1 


( INIZIO ) 


VISUALIZZ, 


Figura 10.6 


FINE 


Figura 10.7 
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100 REM RICERCA DEL NUMERO DELLE SOLUZIONI CHE FERMETTANO 
110 REM DI FORRE N REGINE SU UNA SCACCHIERA SENZA CHE VE 
120 REM NE SIANO DUE "IN FRESA" 

130 DIM DCIO) 


140 N9 = 8 
150 FOR N = 2 TO N9 
160 S= 0 
170 PRINT 
180 FRINT "N = "iN 
190 PRINT 


200 GOSUE 500 
210 NEXT N 


220 STOF 

500I=1 

510 Ni = INT (N / 2) 

520 DCI) = 1 

530 IF I< > 1 THEN 560 

540 IF D(1) «“ = N1 THEN 600 


550 GOTO 690 

560 FORJ= 1 T0I1I- 1 

570 IF DCI) = D(J) THEN 640 

580 IF ABS (DCI) — D(J)) = I - J THEN 640 
590 NEXT 4 


610 IFIE N THEN 510 
620 IF D(N) > D(1) THEN GOSUE 700 
630 I = N 


640 IF DCI) « N THEN 670 
650 I=I-1 

560 GOTO 640 

670 DCI) = DI) +1 

580 GOTO 530 

690 RETURN 

700 FOR L = 1 TON 

710 FRINT DL); TAEC 4 x L)}; 
720 NEXT L 

730 PRINT 

740 RETURN 

750 END 


Figura 10.8 


Per i sistemi che dispongono di un BASIC che accetta le variabili intere (l'’MBA- 
SIC della Microsoft ad esempio), la velocità di esecuzione sarà aumentata dichia- 
rando tutte le variabili intere e sostituendo l’istruzione della linea 510 con N1=N/2. 


10.3 CONCLUSIONI 


Tutti questi esercizi hanno mostrato che la programmazione propriamente detta, 
non presenta in generale difficoltà, dopo aver acquisito la tecnica di effettuare 
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Esempio di esecuzione: 


+INMMNIMUOIN-t Mt 00 0 Nn Innannennsnineoaomno 


NANO LOMN ON +1 LA NN AINNFLE7 NONE MUMIOMZONAN 


ttt LN LNANNANNaetehnaert Ina MAZ4NTNN 


MNAANOTTTAITMONITITATNNTFNOITNODOO0NTNNNDO 


I10IQITLTNMNILAITNOTTNNNEANDOINDIAMNI IMNNDNMO 


909 +N INNIAIMNINNLNQNNNQIONNTNMNNDQA-AMNTN 


NINNTNINDOOINSNNAITINMNNILIIIAIINNAOTNNNNDONNMDO 


Se INNNNNNNNNANOMIOOOMOOOOMMFPtt Ltfttetfrtreer 


uu 
m (al'ellela) 
es) nam 
+ mario 
cd n 
HUH INN 
z z 


Human co o cnr na 
+ 4h MOMNOMZ ON ht 
NonsnmsnsesgsreM Nye 010 
nosgrnynse re 10 Armor o 
Nrct EMO n ao onde 
menor +nnon gone nh 


AI TENNANANNNAUmmmmo 
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presto riguadagnato durante la scrittura del programma e della sua messa 


a punto. 


procuratevelo. 
— Prima di cominciare riflettere bene: il tempo passato nel fare l’analisi è 


— Cercate se tale programma esiste, in caso affermativo, se è possibile, 


AI termine di questo libro vogliamo dare un ultimo consiglio: prima di lanciarvi 


l’analisi del problema e di ricavarne il relativo flow-chart. Questo è il caso dell’ulti- 
nella scrittura di un programma: 


mo esercizio inerente alle regine sulla scacchiera. 


APPENDICE I 


L’ALFABETO IN BASIC 


L’alfabeto in BASIC comprende: 


— le lettere maiuscole da Aaz 
— i numeri da 029 
— i simboli di operazione 

addizione e sottrazione +tet 

moltiplicazione e divisione X e / 

elevamento a potenza(1) f 
— le parentesi (e) 
— i caratteri di confronto 

uguale = 

minore < 

maggiore > 
— i caratteri di punteggiatura 

virgola e punto e virgola i 

punto e due punti PES 

punto esclamativo e interrogativo le? 
— i caratteri speciali 

“spazi‘* 

apostrofo ’ 

virgolette io) 

dollaro $ 


@ Partendo da questo alfabeto si costruiscono i programmi in BASIC. 


@ Alcuni sistemi ammettono delle differenze o dei caratteri supplementari. 


(1) Su alcuni sistemi questa freccia verticale è sostituita da un accento circonflesso. 
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LE COSTANTI E LE VARIABILI NUMERICHE 


Costanti numeriche: Si può rappresentare un numero per mezzo di: 
— un numero intero con o senza segno: 
Esempio: 11, —162 
— un numero decimale senza esponente: 
Esempio: 3.1415927, —3., 0.12, .12 


La parte intera è separata dalla parte decimale per mezzo di un punto . la virgola 
non viene utilizzata a questo scopo. 


— un numero con esponente: 
Esempio: 1E + 5 — 16E— 19 
esponente 
Il primo segno è il segno del numero, il secondo è quello dell’esponente: 
— 1.6 E —19 rappresenta —1,6 10° 


Poichè si è costretti a scrivere i numeri su una linea, si separa l’esponente dal resto 
con la lettera E. 


Attenzione: Il calcolatore distinque il numero 0 dalla lettera O maiuscola. L’utente 
che utilizza la tastiera deve fare la stessa distinzione. 


LE VARIABILI NUMERICHE 


Si distinguono due categorie di variabili numeriche: 


— le variabili semplici, 
— le variabili con indice dette anche “MATRICI“. 
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Le variabili numeriche semplici sono caratterizzate dal loro ‘‘nome‘ (chiamato 
anche identificatore) che è rappresentato: 


— sia da una lettera da A a Z, 
— sia da una lettera seguita da un numero. 


Esempio: A, B, RI, R9 sono delle variabili. 
In un programma si può dunque avere un totale di 26 variabili a una lettera e 26 X 


10 variabili costituite da una lettera seguita da un numero per complessive 286 
variabili distinte. 


NOTA: La maggior parte dei BASIC per microcalcolatori (soprattutto quando 
l’interprete è stato sviluppato dalla Microsoft) accetta dei nomi composti sia da una 
lettera, che da una lettera seguita da un carattere ‘““alfanumerico“ (lettera o nume- 
ro). 


Le variabili indice sono caratterizzate dal loro nome e da uno o due indici. 
Esempio: A(R,I), B(1), C(I + 10 X K) 


Gli indici possono essere delle costanti, delle variabili o delle espressioni aritmeti- 
che. 


Prima di utilizzare una variabile indice, occorre dichiararla per mezzo di un’i- 
struzione DIM all’inizio del programma. 


Esempio: DIM A(10, 20). 


LE ESPRESSIONI ARITMETICHE 


Le espressioni aritmetiche sono costruite a partire da: 


— variabili e costanti, 

— operatori +, —, *,/, }, 

— funzioni numeriche standard 

— funzioni dell’utilizzatore (numeriche), 
— parentesi. 
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Le parentesi svolgono due funzioni: 


— Inquadrare l’argomento o gli argomenti di una funzione. 
— Impostare l’ordine col quale è necessario effettuare i calcoli per 
“valutare‘ il valore dell’espressione. 


Esempio: A + B:C significa a + (b x c) 


(A + B):C significa (a + b) x c 
A + B-SIN(C + 3) qui le parentesi ‘“inquadrano“ l'argomento che è C + 


3. 

Normalmente le operazioni sono effettuate tenendo conto della “priorità degli 
operatori‘ e a priorità uguale nell’ordine da sinistra a destra. L’ordine di priorità 
decrescente è il seguente: 

parentesi, 

funzione, 

operatore elevamento a potenza, 

moltiplicazione e divisione, 

addizione e sottrazione 

Esempio: l’espressione 

A * (B+ 3.2 « SIN(Y + 3 * Z)) + X 1 4/C 
t- 4 
2.1 


4 9 68 


è eseguita nell’ordine indicato e corrisponde a: 


Ax (B+3.2x SIN(y+32)) +. 
E 


ISTRUZIONI DI ASSEGNAZIONE ARITMETICA 


La forma di questa istruzione è: 
variabile = espressione aritmetica 


| variabile semplice o elemento di tabella 
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Il significato è: calcolare il valore numerico dell’espressione situata a destra 
del segno uguale, e racchiudere il risultato nella variabile situata a sinistra del 
segno uguale. 


Esempio: V = 4 * 3.14159 x R 1 3/3 
X= A 


ISTRUZIONI DI SALTO 


Salto ‘‘incondizionato‘ 
La forma più semplice dell’istruzione GOTO è: 


GO TO E 
oppure 


GOTO E 
con E numero di linea. 


Questa istruzione provoca il salto alla linea n. 


Salto ‘‘calcolato‘‘ 
Questa istruzione ha in generale la seguente forma: 


ON espressione numerica GO TO E; ,E:,E3,Ep 
con E: ,E:,...,.E numeri di linea. 


Nell’esecuzione questa istruzione provoca la valutazione dell’espressione, l’even- 
tuale conversione per arrotondamento in un intero e, quindi il salto 


— alla linea E: se il valore ottenuto è 1 

— alla linea E: se il valore ottenuto è 2 

— alla linea E; se il valore ottenuto è p. 
Esempio: ON I GO TO 100, 200, 600, 200 

se I vale 1 salto in 100 


se I vale 2 o 4 salto in 200 
se I vale 3 salto in 600 
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Se il valore dell’espressione dopo l’eventuale arrotondamento non cade in un’in- 
tervallo (1, p), il risultato è incerto: 


— l’esecuzione prosegue in sequenza, oppure 
— la segnalazione di un codice errore. 


Salto condizionato 


I test sono effettuati per mezzo dell’istruzione IF che può assumere differenti 
forme: 


Prima forma: Si tratta della forma più semplificata: 


espressione numero 
IF di THEN di 
comparazio- 
ne linea 


Se la risposta al confronto è si, vi è un salto alla linea indicata, altrimenti 
l’esecuzione prosegue in sequenza. 


Esempio: IF A < B THEN 600 


Il confronto si può effettuare con i seguenti simboli: 


minore n 
maggiore > 
minore o uguale <= 
maggiore o uguale >= 
diverso da <> 


Questi operatori possono essere utilizzati con delle variabili numeriche e con 
stringhe di caratteri. 


Seconda forma: Si tratta di una forma più efficiente: 


espressione 
IF di THEN Istruzione eseguibile 
comparazione Il 


questa istruzione non può essere 
un’istruzione FOR 


— Se la risposta al test è SI allora l’istruzione è eseguita e se questa istruzione 
non è un salto, l’esecuzione prosegue in sequenza. 
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— Se la risposta al test è NO l’istruzione non è eseguita, vi è il passaggio alla 


linea seguente. 


Esempio: IF A < B THEN X=B 
IF A < B THEN GO TO 600 


Terza forma: Questa forma molto evoluta è di più facile utilizzazione. 


espressione istruzione 
IF di THEN ELSE 
comparazione eseguibile 
Esempi: 


IFA<BTHENC=A:D=BELSEC=B:D=A 
IF A$ = B$ THENI=1I+ 1 ELSE GO TO 100 


NO 
CONFRONTO 


ISTRUZIONE 
ESEGUIBILE 


ISTRUZIONE 
ESEGUIBILE 


Questa istruzione corrisponde al flow-chart sopra riportato. 


istruzione 
eseguibile 


Per alcuni sistemi ogni alternativa di un’istruzione IF può contenere una sola 


istruzione. 
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LOOP DI CALCOLO 


Si realizza un loop di calcolo per mezzo dell’istruzione FOR ... NEXT. 


FOR var = El TO E2 STEP E3 


o ooo 


NEXT var 


nella quale var è il nome di una variabile numerica: 


E1, E2 e E3 rappresentano delle espressioni aritmetiche 
EI dà il valore di partenza 
E2 dà il valore finale 
E3 dà il passo 
Se = El < E2 deve essere E3 > 0 
E] > E2 deve essere E3 < 0 


i valori numerici di El, E2 e E3 sono calcolati prima di ‘“incominciare‘ il loop di 
calcolo. 
Se il passo è I si può scrivere 
FOR var = El TO E2 
Esempio: 100 FORI= I TO 10 
110 FORJ= 1 TO 10 
120 A(I,J) = 0 


130 NEXT J 
140 NEXT I 


STRINGHE DI CARATTERI 


Costanti di caratteri 
Sequenze di caratteri fra virgolette. 


Esempio: ‘“ABCD“ 
“BELLA DONNA“ 
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Il carattere ‘‘spazio‘ è significativo in una stringa di caratteri. La lunghezza 
massima di una stringa di caratteri dipende dal sistema utilizzato. 


Variabili di caratteri 
Il nome di una variabile è definito da una lettera seguita dal carattere dollaro 
Esempio: A$, B$ 
dunque un massimo di 26 variabili distinte in un programma. 
Operazioni possibili sulle stringhe di caratteri 


Confronti: Si dice che A$ minore di B$ se nell’ordine alfabetico la variabile A$ 
precede B$. 


Esempio: A$ = “GIANNI“ 
B$ = “ROSSI 


qui B$ minore di A$ 
Concatenazione: consiste nell’affiancare due stringhe di caratteri. 


Esempio: A$ = “PIETRO” 
B$ = “ VERDI“ 
N$ = A$ + B$ il carattere + è utilizzato in BASIC 


Fircone simbolo di concatenazione 
N$ assume il valore “PIETRO AVERDI"“. 


Funzioni speciali per stringhe di caratteri: L’elenco delle funzioni non è limitativo e 
può variare da un sistema all’altro. La seguente lista rappresenta numerose funzio- 
ni disponibili sui microcalcolatori. 


ASC (X$) fornisce il valore numerico del primo carattere di X$ 
secondo la rappresentazione del codice ASCII (codice 
ISO) per esempio ASC (“A°*) dà 65. 


CHRS$(1) fornisce una stringa di caratteri che rappresenta in 
ASCII il valore numerico di I. 


IN$TR(I,A$,B$) cerca la prima applicazione di B$ in A$ partendo dalla 


I-esema posizione. 
Fornisce la posizione della sottostringa B$ trovata. 
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LEN(A$) 


CHRS$(I) 


HEXS$(X) 


LEFTS$(A8$,1) 


MID$(A$,I,J) 


OCT$(1) 


RIGHTS$(A$,I) 


SPACES(I) 


STRING$(I,J) 


VAL(A$) 


fornisce la lunghezza della stringa A$ 


fornisce una stringa di caratteri che rappresenta, in 
codice ASCII, il valore di I. 


fornisce una stringa di caratteri che rappresenta in 
esadecimale il valore decimale dell’argomento X. 


fornisce una stringa di caratteri contenente gli I caratte- 
ri più a sinistra di A$ 
fornisce una stringa estratta da A$ partendo dal I- 


esimo carattere e di lunghezza J. 


fornisce una stringa che rappresenta in codice ottale il 
valore decimale dell’argomento I. 


fornisce una stringa contenente gli I caratteri situati più 
a destra nella stringa A$ 


fornisce una stringa contenente gli I caratteri ‘‘spazio‘‘. 


fornisce una stringa di lunghezza I e rappresenta il 
valore numerico di J. 


fornisce un numero che contiene il valore numerico 
della stringa ASCII A$. Questa funzione presuppone 
che i caratteri contenuti in A$ rappresentino un nume- 
ro. 


INPUT/OUTPUT 


Lettura di dati ‘‘conversazionali‘‘: 


INPUT lista di variabili 
variabili semplici o elementi di tabelle 


Lettura dei dati inclusi nel programma: 


READ lista delle variabili 
DATA valori numerici separati da una virgola (o da uno spazio in certi 


sistemi) 


RESTORE per ‘riavvolgere‘. 
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Visualizzazione dei risultati: 


Prima forma: 
PRINT lista 


| delle espressioni 
delle costanti 


Esempio: 100 PRINT “X = ©“; X,% Y =“; Y 


I separatori da utilizzare sono la virgola e il punto e virgola. La virgola non limita 
la lunghezza della zona utilizzata per la stampa. Il punto e virgola limita questa 
zona alla lunghezza necessaria per la spampa. 


Un separatore alla fine dell’istruzione PRINT impedisce il passaggio alla linea 


seguente. 


Nella lista, si può inoltre utilizzare, la funzione TAB che permette di inserire delle 


tabulazioni. 
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Non si può imparare la programmazione in BASIC 
senza la pratica. Ecco dunque in questo volume 
una raccolta completa e progressiva di esercizi ri- 
guardanti matematica, gestione, ricerca operati- 
va, gioco e statistica. Ciascun esercizio proposto 
comporta l'enunciazione e l’analisi del problema, 
la risoluzione mediante flow-chart e commenti, 
così come un programma che implementa la solu- 
zione illustrato da semplici esempi rappresentati- 
vi. Questo metodo mette in grado il lettore di veri- 
ficare passo passo le sue conoscenze e il livello di 
apprendimento raggiunto, nonchè di risolvere in 
modo autonomo i problemi secondo un procedi- 
mento ‘inverso’. Cioè, come è possibile scom- 
porre il problema in sottoproblemi ben identificati 
e risolvibili separatamente in modo semplice. Il let- 
tore è così in grado di risolvere ogni difficoltà e 
costruire un programma modulare di facile lettura 
e adattabile alle specifiche applicazioni. Tutti i 
programmi, al fine di renderli di validità generale, 
sono scritti in BASIC Microsoft e girano su AP- 
PLE, PET/CBM, TRS80 e altri computer. 
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